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Mechanicam vero duplicem Veteres constituerunt: Rationalem quae per De- 


monstrationes accurate procedit, & Practicam. Ad practicam spectant Artes omnes 
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operatur, imperfectior est Mechanicus, & si quis accuratissime operart posset, hic 


foret Mechanicus omnium perfectissimus. NEWTON 
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plutét les véritables loix de la Nature dans tout leur éclat, & on y trouvera des raisons 


encore plus fortes, d’en admirer la beauté & la simplicité. EULER 


Ceux qui aiment l’Analyse, verront avec plaisir la Méchanique en deveniy une 


nouvelle branche... LAGRANGE 
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Uber die Potentialfunktionen von Kurvenbelegungen 


CLAUS MULLER & RoLF LEIS 


Das Verhalten der Potentiale einfacher Kurvenbelegungen bei Annaherung 


_ des Aufpunktes an die Kurve wurde zuerst von H. Porncaré [1] behandelt und 


am Beispiel eines homogenen gradlinigen bzw. kreisférmigen Drahtes explizit 


_angegeben. Da diese Frage in der Zwischenzeit nicht weiter behandelt wurde, 
_ erscheint es reizvoll, die Problemstellung wieder aufzugreifen und auf das Ver- 
_halten der Potentiale von Dipolbelegungen auszudehnen. 


Ziel der vorliegenden Arbeit ist es daher, das singulare Verhalten der Potentiale 


von beliebigen Kurvenbelegungen bei Annaiherung des Aufpunktes an die Kurve 


zu beschreiben. Dazu stellen wir die Laplacesche Potentialgleichung AU =0, 


- der die Potentialfunktionen geniigen, an den Ausgangspunkt der Uberlegungen 


und bestimmen aus ihr durch Reihenansatz Lésungen, die ahnliches Verhalten 
zeigen wie die Potentialfunktionen. Durch einen Grenziibergang gelingt es dann 
mit Hilfe der Greenschen Formel zu zeigen, daB die so gewonnenen Lésungen bis 
auf eine tiberall stetige Funktion mit den Potentialen iibereinstimmen. Man 
kann daher an den Potenzreihen das singulare Verhalten der Potentialfunktionen 


~ ablesen. Die wichtigsten Ergebnisse lassen sich dann folgendermaBen formulieren: 


ej 


Mit y werde ein innerer Punkt der Kurve © beschrieben. Es seien und 6 
die Vektoren der Haupt- und Binormalen. Dann verhalt sich bei Annaherung 
an einen inneren Punkt y der Kurve © das Potential 


G ds 
| 0(s) Tr—y(s)) 
€ 


der einfachen Belegung @(s) wie 
—2e(»)Ig|y—z|+0(1) fir zp. 


Gibt der Einheitsvektor 
D(s) = sina h(s) + cosa b(s) 


die Richtung der zur Kurventangente senkrechten Achse des Dipols an, so stellt 


ds 


Jee pines At 


das Potential der ees dar. Diese Funktion verhalt sich bei Annaherung 
an den inneren Punkt y wie 


20(0) (d(0) ra) — eC) *(o) (68) Ig |z — vl + (4) 


wobei x(y) die Kriimmung im Punkte } darstellt. 
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In $1 der Arbeit wird zunachst das singulare Verhalten der Potentialfunk- 
tionen von geraden Linien explizit angegeben. Da das kartesische Koordinaten- 
system zur Darstellung von beliebigen Kurven wenig geeignet ist, werden 
in § 2 krummlinige Koordinaten eingefiihrt, mit deren Hilfe dann in § 3 und § 4 
das singulire Verhalten der Potentialfunktionen von einfachen und doppelten 
Belegungen in inneren Punkten der Kurve behandelt wird. In §5 wird ab- 
schlieBend noch das Verhalten an den Endpunkten nicht geschlossener Kurven 
untersucht. 

Das in der vorliegenden Arbeit behandelte Problem zerfallt in zwei ver- 
schiedene Fragenkreise. Zunachst handelt es sich namlich darum, den expliziten 
Ausdruck zu bestimmen, der das Verhalten in der Nahe der Kurve bestimmt. 
Dann ist zu priifen, unter welchen Differenzierbarkeitsvoraussetzungen die ge- 
wonnenen GesetzmaBigkeiten gelten. 

Die untersuchten Kurven und ihre Belegungen werden dementsprechend zu- 
nachst als analytisch vorausgesetzt, um die formale Seite des Fragenkreises nicht 
mit Differenzierbarkeitsbetrachtungen zu belasten. Nach Erhalt dieser Ergeb- 
nisse ist es dann leicht, die Voraussetzungen auf Differenzierbarkeitsbedingungen 
endlich hoher Ordnung zu reduzieren. 


§ 1. Vorbemerkungen 


Unseren Betrachtungen liege der dreidimensionale euklidische Raum zugrunde. 
In ihm sei in einem kartesischen Koordinatensystem eine endliche doppelpunkt- 
freie Kurve © in der Form 


(1.1) ) = 9(s) 


mit beliebig oft differenzierbarem (s) gegeben. Der Parameter s sei die Bogen- 
lange dieser Kurve. Der Punkt y liege nicht auf €. Dann verstehen wir unter 
der Potentialfunktion U(r) einer einfachen Belegung @,(s) der Kurve € die 
Funktion 


12 Uy () = f sis 
(1.2) 0 (t) J Teo)” 


Die Belegung @9(s) wollen wir ebenfalls als beliebig oft differenzierbar voraus- 
setzen. Es sei §(s) der Hauptnormalen- und b(s) der Binormalenvektor der 
Kurve ©. Es sei ferner 9;(s) eine beliebig oft differenzierbare Funktion, dann 
erhalten wir durch 


(153) Ue) = f als) 9(9) Fy ras 
¢ 

bzw. 

(1.4) Ut) = f 02(9) (3) Ky ds 
¢ t—)| 


die Potentiale der Dipolbelegungen. Dabei ist 


e, 


1.5 EG) 4) 
(1.5) V = ——e,+ i eta 
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der ,,Nabla‘“-Operator. Durch das Symbol VY, deuten wir an, daB beziiglich 
differenziert wird. Es gilt 
1.6 oc eee 
ee "Te ‘lz 9| 

Die Funktionen U(r) (j =0,1, 2) sind fiir alle x, die nicht auf der Kurve 
legen, differenzierbar. Bezeichnen wir mit 


a o2 o2 o2 
(1.7) 4=24+5+4 
den ,,Laplace“-Operator, so folgt fiir alle Punkte r¢€ 
(1.8) AU (t) = 


Unsere Potentialfunktionen geniigen also der Laplaceschen Differentialgleichung 
A= (), 

Um einen ersten Einblick in das Verhalten der Potentialfunktionen bei 
Annaherung von xy an die Kurve © zu erhalten, betrachten wir das Potential 
U(x) fiir die Gerade y=0, z=0, OX xSL mit konstanter Belegung 0). Wir 
fiihren in der (y, 2)-Ebene durch 


(1.9) Ve 7 SIVGu = a7 COS |, 


Polarkoordinaten ein und erhalten 


iL ae 

vt) = == 

) 2 Vé—4)?2 +7 
= @)In|(L—*%+V(L — «)? +P) (« +Vx?+ 7)] — 2e9Inr. 

Fiir 0<*<L folgt aus (1.10), daB U) (x) bei Annaherung von x an die Kurve € wie 


(4.10) 


(1.14) —20,Inr 


singular wird. Fiir U,(r) und U,(r) erhalten wir 


1 L—x iG 
2) U;, (x 1 rsin 9 i aaa ae | 
(1.12) U( af yg Vé—9+" d= = 0, r \V@aonepr ate \err 


und 


¥ COS P -d c. ee | L—%* gz . 
(1.13) G(r) = =a { (gS yep pe E = 0» r \Vtanetr ts} 


Daraus folgt, daB die Funktionen U,(r) bzw. U,(r) fiir 0<*<L wie 


(1.14) Die DAW ene On 


if 


COs P 
ey cia 
singular werden. 

Wir werden erwarten, daB diese Ergebnisse auch bei beliebig geformten 
Kurven und Belegungen die Anfangsglieder einer Reihenentwicklung der Poten- 
tialfunktionen darstellen. Um das zu zeigen und um weitere Glieder dieser Ent- 
wicklung anzugeben, wollen wir die Kurve zunachst in einem geeigneten Koordi- 
natensystem darstellen und formale Beziehungen in diesem Koordinatensystem 
ableiten. 


afi 
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§ 2. Das Koordinatensystem 
Wir fiihren in der Normalenebene der Kurve © Polarkoordinaten 7, p ein 
und kénnen dann durch 


(2.1) r(s, y, 7) =(s) +rsin y h(s) +7 cos  b(s) 


einen Punkt ¢ im Raum festlegen. Um ein Orthogonalsystem zu erhalten, setzen 
wir noch [r(s) sei die Windung von (| 
(2.2) Tis\=—rOat 
So 
und verwenden fortan die Koordinaten 


(2.3) W=s; wW=g—T(s); wr. 


Aus den Frenetschen Gleichungen 


t’ # t 
(2.4) pia | =e ihe 
b’ —T b 
folgt wegen 
Bites Oe ee 
(20>) out ieee Oy’ ou? OQ’ ous er 


unmittelbar mit 6d =1—7~x sin 
0 
Ee fete ae ot, 
(2.6) tase E =7 cos ph —rsin pb, 


é : 
ewer r=singh-+cos pb, 


so daB sich 
07-0).0 
(2.7) Sik = O 7 O 
OROm 


ergibt. In einer Umgebung der Kurve © ist die Darstellung (2.1) umkehrbar 
eindeutig, wenn dort g = det(g,,) nicht verschwindet. Fordern wir daher 


(2.8) r>0;  rxsing<1, 


so ist diese Bedingung erfiillt. Wir wollen unsere weiteren Uberlegungen deshalb 
auf das schlauchartige Gebiet © beschranken, welches die Kurve © und alle 
Punkte, die der Bedingung (2.8) geniigen, enthiilt. 

Wir betrachten nun die in ® enthaltene Koordinatenflache u3= yr = const. 
Der Ma&tensor dieser Flache Z hat die Gestalt 


(2.9) ve ( " 
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und das Flachenelement ergibt sich wegen 


A (ut, u2) Aer 

2.10 Biya 1 
p19) a(s, 9) 0 4 

zu 

(2.11) dF =rddwdw=rddsdpo. 


Den Normalenvektor n der Flache Z erhalten wir dutch 


4 : 
(2.12) Ni Vaer Ce ene et EPP Bia 
ap 


Es sei g(r) eine differenzierbare Funktion der Ortskoordinaten. Dann gilt auf Z 


(2.13) 


SchlieBlich a wir noch den A-Operator im u-System an. Es ist 


AU =—_ 2 git 
ae es a 
1 Geen ra) 1 4) a) 1 re) é 
(2.14) ~ 6 dul 6 Bul Use v26 Ou® Ou? ot 7d our ous 
; ee cos # sin 
= Uji (x pees eh Ui» 5 Us 


Ui Uiee 
7 é2 ih ine Y wr Us. 


Wir wollen aus (2.14) durch Reihenansatz in  periodische Lésungen be- 
stimmen, deren erste Glieder ahnliche Gestalt haben, wie wir es uns am Falle 
der geraden Linie im §1 tiberlegt haben und erhalten 


. 1 Mt C 
Vo (tz) = — 200(s) Igr —xgsingprigr +r lgr(oo — x? 09 sin? y) 


V, (x) == a(s) sin @ + 0, x sin? m — 


—7lg7 {sin yp (or ~ ae 01 i) + cos y(o,t' +27 o1)} + 


01%" 


(2.15) +ir2o, sin? p + — oo 01 + a1 7°) 


Va (x) == ols) cos p + + o»x sin 29 — 


ah ey {eos y (e2 4 oe Q2t”) — sing(Q,t’+ 27 os) oe 


— 2-1729, cos + 7 rigrzu(d.t’ + 27 03) 


mit 
(2.16) AV, (x) =O(4) fiir Kleine 7. 


Im nachsten Paragraphen werden wir zeigen, da unsere Potentialfunk- 
tionen U,(r) sich fiir 70 wie V;(x) verhalten. 
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§ 3. Das singulare Verhalten der Potentialfunktionen von einfachen 
Belegungen 
Im letzten Paragraphen haben wir die Funktion V(x) eingefiihrt, die der 
Abschatzung 


(3.1) AVs(t) =O(4) tordcieine 7 


gentigt. Um zu zeigen, da8 diese Funktion sich in der Nahe der Kurve (fiir 
kleine 7) wie U,(x) verhalt, fiihren wir folgende Bezeichnungen ein: 

Die Oberflache unseres Gebietes © sei Q. 

€  @ enthalte das Gebiet 8 mit der Oberflache Z 

(u? =r= const). Dann wird das Gebiet ¥ = 6 — 3 

von den Oberflachen 2 und Z berandet. In ¥B 


y wahlen wir einen festen Punkt ry, aus, um den 
wir eine geniigend kleine Kugel & mit dem 

tr, Radius @ und der Oberflache K legen. 8 liege 

Fig. 1 ganz in W. yx sei ein beliebiger Punkt in G, 


R(x, Xp) der Abstand zwischen r und Yg. 
Wenden wir nun die Greensche Formel auf das Gebiet 8 —S mit den Funk- 
tionen V(x) und R1(r, x9) an, dann erhalten wir 


[led p— RAM aY = |b de ee Al aP 


RB 
Bee) a. 1 @ ° a1 4 °@ 

pit én RR an Wi) aF — | (% én RR €n K)aF. 

K 
Aus 
go - 3 AO, -_ \Y als 

(3.3) Ges R R an) OF & o : oO ae) Gao 

K K 
folgt 

3 A 4 1 0 
(3.4) lim [ (esa pe oy GF = — 40 Yu). 
Da das Integral 
1 1 

05) f(a h—4anjav 


8 
existiert, erhalten wir aus (3.2) 


— fx atavmsatiod + [OS ae Wear 


6-6) d oo Me ee ey, 
~ an, Ree a) aE. 


Mit r= Yy— y und (vgl. Fig. 4) 


(3.7) 1 =t—r(singyh+cospb)=r—rn 
folgt 


(3.8) R? =|r)— y [2 — 2ren4+72, 
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und es gilt bei festem rz, fiir kleine 7 


1 1 ytn 
ne) — | O 2 : 
ae Spear aires See 
Wegen 
AO eee ne 
8-10) Ito |? r= H|’ 


erhalten wir schlieBlich 


i aa VT aI 


eee 
eg ees aT ae 


OW); 
(3.41) 


Aus der Darstellung (2.15) der Funktion (ry) entnehmen wir 


Vo = — 2e9lg7+ (1) 
(3.12) Ay: ' 


1 
ee = 209 + O(lg7). 


Zusammen mit (3.11) folgt daraus fiir v0 nach (1.2) 


(3.13) Jone ee ifr 


da wegen der vorausgesetzten endlichen Lange der Kurve alle anderen GréBen 
bei festem ry beschrankt bleiben und bei nicht geschlossenen Kurven die Beitrage 
der Endflachen wegen dl =r dyrdq von der Ordnung o(1) sind. 


ae ee )= 42 Up (Xo) +0 (1), 


Wegen 
(3.14) AV, =O (1) 
existiert das Integral 
AN 
(3.15) [ dV 


S 


fiir alle xy), und wir erhalten aus (3.6) durch den Grenziibergang 7y—0 fiir alle 
Yo, die nicht auf € liegen, 


AV, GA I  @ 
(3.16) 42 Uj(%o) = 42 Vo (Xo) 4 ; =p oy ! { (% én R RR @n 1) 4F. 
6 Q 


In (3.16) ist das Integral iiber 2 eine stetige Funktion fiir alle inneren Punkte ry 
von ®&. Da AX iiberall in & beschrankt ist, stellt auch das Integral tiber G fiir 
alle r€® eine stetige Funktion dar [2]. Fir alle inneren Punkte von ©, die 
nicht auf © liegen, stimmt das Potential U,(r) daher bis auf eine tiberall stetige 
Funktion f,(z) mit M(x) iiberein. Durch den Grenztibergang 


(3.17) Lim {Uo(ts) —V(t,)}= cee to(t») = fo (Ge) 


erhalten wir daher 
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Satz 1. Es se 
Vo(t) = — 200(s) lgr — x gosing rlgr + 37 lg7 (oo — %* eosin’ ¢). 
Dann gibt es eine iiberall in & stetige Funktion fy(z), so daB fir alle x in & 


Up (t) = Vo(t) + fol) 
gilt. 
Fiir r>y€ verhalt sich Uj (rx) daher wie 


(3.18) — 209(y) lg |x — 5]. 


§ 4. Das singulare Verhalten der Potentialfunktionen von Dipolbelegungen 

Um das singulare Verhalten der Potentiale von Dipolbelegungen anzugeben, 
gehen wir ahnlich vor, wie bei den einfachen Belegungen. Zur Vereinfachung 
der Rechnungen fassen wir die Funktionen U,(r) und U,(x) in 


(4.4) Om sine Ort ae Oe (ee Oe ereriet 

mit D(s) = sina h(s) + cosa 6(s) zusammen. Entsprechend betrachten wir 
(4.2) V=sinaKW+cosah, 

mit 

(4.3) AV =O(1) 


und wollen zeigen, daB die Funktionen U und V bei Annaherung an die Kurve 
gleiches singulares Verhalten aufweisen. 

Dazu wenden wir analog (3.2ff.) die Greensche Formel mit der Funktion 
V(r) an und erhalten fiir alle r,¢€ 


AN Ae : a1 4 
Re 4x V(e) + | (No en ie oe on) OE 
(4.4) me 
1 IL o@le 
~ | (% on R R om) OF 
Nun gilt wegen (2.12), (3.14) und 
OV 
(4.5) Gn = — 20608(p — a) + O(lgr) 
fiir kleine 7 
QO il i Gs 
{vZ4 R agltF = ‘[feonte p—a)n7 eae aa 
Z 
(4.6) += rd cos (y — yay 
je 2 ea 


2 cos(@ — d 
Ife a) sin px =a pds +0(1) 
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oder 


a] V ; . 
(4.7) [Vv 5 Fag) UF = 47 Ute) — 2msina. Ua(to, % e1) +0(1), 
Zz 


on R R an 


wobei wir unter U)(rz,@,) das Potential mit der einfachen Belegung xo, ver- 
stehen. Durch den Grenziibergang y->0 folgt daher aus (4.4) 


42 U (to) =42V (to) + 2asina Uj (to, x 01) + 
(4.8) AV Of 3 qe 
| i R eg {lv én R R én V) ar. 


2Q 


Die Integrale stellen wieder stetige Funktionen dar, es gilt daher 


Satz 2. Es sez 
V(e) =sina V(x) + cosa Vy (xr) = a3) (sin a + cosa) +--: 


die durch Gl. (2.15) defimierte Funktion. Dann gibt es eine tiberall in & stetige 
Funktion f(x), so daB fiir alle x in & 


U(t) =V(z) + sO, % @) sine + f(z) 
gill. 
Bei Annaherung des Punktes xy an den Punkt y der Kurve € wird die Potential- 
funktion 
1 


(4.9) TU) = f 066) 20) Fo pag A 
daher wie 

O(p) (g—9) =, 
(4.10) o(y) {ae — x (0) (6) le le — ol} 
singular. 


§5. Das Verhalten der Potentialfunktionen an den Kurvenenden 


Wir haben bisher das singulare Verhalten der Potentialfunktionen bei An- 
naherung des Aufpunktes y an einen inneren Punkt der Kurve © diskutiert. 
Wir setzen nun voraus, daB € nicht geschlossen ist. Die Endpunkte seien durch 
s=0 bzw. s=L gegeben. Wir wollen das Verhalten bei Annaherung von ry an 
einen dieser Punkte bestimmen und setzen zur Abkiirzung x (0) =o. 


Fir alle Punkte r€® haben wir die Gleichung 


(5.4) Uy (x) = Vo(t) + fo(t) 

mit 

62) =e] f avs [Cae ele 
105) Q 


bewiesen. Es sei yp ein Punkt mit den Koordinaten (0,7, @). % liegt also auf 
der Grundflache K der schlauchartigen Flache 2. Wir konnen das Flacheninte- 
gral iiber 2 in Integrale iiber K und itiber Q*=Q2—FK zerlegen. Die Flache 


96 Craus MULLER & RotrF LEIs: 


K (s=0; 0SrSc; 0X pySX2za) ist eben. Ihr Flachenelement lautet 


(5.3) aF =rdrag. 
Auf K gilt 

C, gte ie a 
2) CORE 


Zur Diskussion des Verhaltens an den Kurvenenden fiihren wir zunachst in 
(5.1) den Grenziibergang y—>r) durch und gehen anschlieBend mit ry gegen ov. 
Bei beiden Grenziibergangen bleiben das Gebietsintegral und das Flachenintegral 
tiber Q* beschrankt, da sie fiir r€ K stetige Funktionen darstellen. Wir disku- 
tieren daher die Integrale 


i 1 One 
(5.5) [Gata wd [KO Za. 
K K 
Die Integrale in (5.5) sind Potentialfunktionen mit einfacher ee doppelter 


Flachenbelegung. Das Potential mit der einfachen Belegung — ae bleibt 


beim Grenziibergang r—yzy=0 stetig, da die Belegung in einer oe des 
Punktes yg beschrankt ist [3]. Wir erhalten also 


ey pe 1 
a) meee as 0 eek ie as "OR, ah 
K 


und wollen zeigen, daB8 das Integral auf der rechten Seite von (5.6) fiir alle ry 
(auch fiir yy = 0) beschrankt bleibt. Es gilt 


lgvdF 
le eelh 


(5.7) fas gs FoR dF SA ts 


Nuniist mit.) r5 |= % 


lev dF |] ‘leg [lg 7| 
5.8 8 =e. g 
12:8) rte roe ea ee 
d |7—17,|S 47 |r— eae 
un 

lg v| 

Trove] OF SA 70 [18 701 

1—1.|S37% 

(5.9) | Is 


Pani Pee [Psire~ ede e a ro|lg 70]. 


Wir erhalten also 


(5.10) 7 os OR 


0 
oes dF) SB, +B, r|lg 75] <C. 


K 


Fur das Integral mit der Dipolbelegung VY folgt beim Grenziibergang rr, [3] 


: OA 


to 
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Es folgt also 


(5.12) Oy (Lo) = & Mo (Lo) + & (Lo) 
mit 
(5.13) (eit) K < co. furvalle 4. 


Fiir ty verhalt sich die Potentialfunktion U,(x)) daher wie 4% (1). 


al 
2 
Fiir die Potentialfunktionen mit Dipolbelegungen U,(r) und U(r) gelten ent- 
sprechende Uberlegungen. An Stelle von (5.7) erhalten wir jetzt 


(5.14 {GA ge SA, feet 
a i (3 os | Sala r|z—2o| 


K 


sin pdF 


J v|z—tol. 
K 


+ A; 


Die beiden Integrale auf der rechten Seiten von (5.14) verschwindén aus Sym- 
metriegriinden. Es gilt z.B. 


5.15 singdF _ yadndy Sat 
( ) ee lena OF 


und es folgt, daB die Potentialfunktion U(r) fiir ry—>p wie 
(5.16) 3 {V (2) + 3 U, (Xo, % 01) sin a} 
singular wird. Wir erhalten also 


Satz 3. Es set U(r) das Potential einer einfachen oder doppelten Belegung der 
nicht geschlossenen Kurve€. Bet Anndherung an innere Punkte von © verhalte 
sich U(x) wie V(x), wober V(x) eine der in Satz 2 bzw. Satz 3 genannten Funktionen 
darstellt. Liegt der Punkt ¢ in der Normalenebene durch ein Kurvenende, so wird 
das Verhalten des Potentials U(x) ber Anndherung an den Endpunkt durch dte 
Funktion $V(t) beschrieben. 


§6. Abschwachung der Differenzierbarkeitsvoraussetzungen 

Zur einfacheren Formulierung der Ergebnisse haben wir bisher von der 
Kurve © und den Belegungsfunktionen @ beliebig hohe Differenzierbarkeit vor- 
ausgesetzt. Da wir diese Voraussetzung jedoch nicht benutzt haben, wollen wir 
sie nun auf Differenzierbarkeitsbedingungen endlich hoher Ordnung abschwachen. 

Bei der Bildung der Funktion /;(r) sind die GréBen g/Y und 1’ aufgetreten. 
Hohere Ableitungen wurden niemals bendtigt. Es folgt daher, daB unsere Ergeb- 
nisse gelten, wenn wir die Kurve € als sechsmal und die Belegungsfunktionen 9; 
als viermal stetig differenzierbar voraussetzen. Wenn wir auf die Stetigkeit der 
Funktionen /;(x) verzichten und nur noch ihre gleichmaBige Beschraénktheit 
zeigen, lassen sich diese Voraussetzungen weiter verringern. Wir beweisen 


Satz 4. Die Kurve © sei dreimal stetig differenzierbar. Die Belegung 0 (s) sez 
hélderstetig. Dann bleibt in jedem ganz in © enthaltenen abgeschlossenen Telbereich S* 


Up (t) + 2@0(s) Ig |e — 9| 
gleichmafig beschrankt, 
und 
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Satz 5. Die Kurve © sei viermal stetig differenzierbar. Die Belegung @(s) set 
hilderstetig differenzierbar. Dann blevbt im jedem ganz in & enthaltenen ab- 
geschlossenen Tealbereich S* 

Ue) — 2o(u) OY + o(v) x (0) (59) le le — 9| 
gleichmafig beschrankt. 
Wir beweisen zunichst Satz 3 fiir konstante Belegungen @,(s). Dazu setzen wir 


(6.1) Vo" (t) = — 2aolg7. 
Vo*(x) geniigt fiir kleine ry der Abschatzung 

1 
(6.2) AV(*(t) =O (=). 


Es existiert daher das Integral 


(6.3) dee vol [ a ee U/ ar) 49), 


so daB unsere Uberlegungen im § 3 richtig bleiben, wenn wir *(r) statt Y%(r) 
verwenden. In f(r) gehen nur Ableitungen von © bis zur Ordnung von ~ ein. Bei 
der Einfiihrung des u-Systems hatten wir jedoch t benutzt, so daB € als dreimal 
stetig differenzierbar vorausgesetzt werden muB. Satz 3 ist damit fiir konstante 
Belegungen bewiesen. 


Es sei @9(s) nun eine hdlderstetige Funktion. Dann gilt in einem beliebigen 
Punkt x) € © 


(6.4) Uy (aps 00 (s)) =U, (ie 00 (So)) + Uy (to; Qo(S) — Qo (So)) , 

und es ist 

(6.5) Up (x, 20 (So)) + 2.9 (So) lg 7 

fiir alle x gleichmaBig beschrankt. Wir miissen daher zeigen, daB auch 
(6.6) Ug (vo, 0o(S) — Qo (So)) =U, (fo: 00 (s)) 


gleichmaBig beschrankt bleibt. Dazu stellen wir unsere Kurve in der Umgebung 
des Punktes ) in einem kartesischen Tangentensystem dar. Dieses hat folgende 
Eigenschaften: 


1. Yo ist Ursprung des Systems. 
2. Es gibt eine Konstante c>0 so, da der im Innern der Kugel 
2% 4 4% 1 22 < ¢ 


gelegene Teil ©, von © in der Form 


y=f(x);  z=g(x) 


dargestellt werden kann, wobei die Funktionen / und g in x<c zweimal stetig 
differenzierbar sind. Die Konstante c existiert fiir alle Punkte Yo gleichmaBig. 


3. Die Ableitungen /’ und g’ verschwinden im Ursprung. 
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Wir wahlen nun eine feste Zahl 0<1<c und zerlegen Up (tp, of (s)) im Punkte 
Yo=%o —7 Sing hy—7 cosy by in der Form 


(6.7 Oem) = Qa (s) ds 1 e8(s) ds 
) peeinee) | yea tol taco 


Tt 


Die Integrale auf der rechten Seite von (6.7) sind beschrankt und kénnen fol- 
gendermaBen be EAs werden: 


r= f (s) ds 
¢ [do —h+7sin p h)+7 cos @ by| 
_ |a|*d 
—O : x x 
aol SE O(z*), 
6.9 oo (s) ds | _ eas 
| | bi ‘lto— 5] cy Ses 


Es gilt also gleichmaBig beziiglich yy mit 0<a#<1 
(6.10) |% (to, 8(5))| S4a(2 +14) Sy, 


womit Satz 4 bewiesen ist. 
Zum Beweis von Satz 5 gehen wir entsprechend vor. Wir setzen 


—w) d] a 
lc—oP 


(6.11) V(x) =g(o) {77 ) — 2(p) (hd) lg 2-9] + (0) 4 


V*(x) gentigt fiir kleine y der Abschatzung 


(6.12) AV*(z) =O (} 


7 

Es existiert daher das mit V* gebildete Integral (6.3) und unsere alten Uber- 

Jegungen bleiben richtig, wenn 9 (s) zweimal und € viermal stetig differenzierbar ist. 
Es sei o(s) nun hélderstetig differenzierbar. Dann bilden wir in einem belie- 

bigen Punkte r,C & 

(6.13) U(%o, 0(s)) = U(to, 0 (s)) + U(to, o*(s)), 

mit 


(6.14) 


Die Funktion o*(s) geniigt wegen der vorausgesetzten hdlderstetigen Differen- 
zierbarkeit der Abschatzung 


(6.15) |o*(s)|<|s —so] lo'(so + 68 — 50)) — oso) S|s— sl", O< a1. 


Da o!(s) zweimal differenzierbar ist, bleibt U(r, o!(s)) — V(x, e!(s)) nach obigen 
Uberlegungen gleichmaBig beschrankt. Um auch U(r9, e*(s)) abzuschatzen, be- 
nutzen wir wieder im Punkte y, das Tangentensystem und zerlegen 


(6.16) U (to, 0 qe peo ne Td Tae 
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Die Integrale der rechten Seite von (6.16) sind beschrankt und kénnen folgender- 
maBen abgeschatzt werden: 


\a[ te? da 2 
(6.17) ) fers =o f yx See ee ) 


ey 1 
(6.18) 1 f et(s) Vas 


Es gilt daher gleichmaBig beziiglich ry 
(6.19) |U (eo, e*(9))| S Ba (se +0") S Ba, 


womit auch Satz 4 bewiesen ist. 
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Uber das Neumannsche Randwertproblem 
fir die Helmboltzsche Schwingungsgleichung 


ROLF LEIS 


Vorgelegt von C. MULLER 


In der Potentialtheorie sind drei Randwertprobleme des AuBenraumes be- 
kannt geworden: die Dirichletsche Aufgabe (der Wert der Lésung U ist auf dem 
Rand der geschlossenen, orientierbaren und glatten Flache Q als stetige Funktion 


vorgegeben), die Neumannsche Aufgabe (oe ist vorgegeben) und die dritte 
Randwertaufgabe (So +8u ist als stetige Funktion vorgegeben). In allen 


Fallen gelingt es, die Aufgaben auf Integralgleichungen zuriickzufiihren und die 
Lésungen als Potentialfunktionen von einfachen und Dipolbelegungen darzu- 
stellen [7]. 
Eine Ubertragung dieser Ergebnisse auf die Theorie der Helmholtzschen 
Schwingungsgleichung 
AU +kU=0 (k+= 0, Imk= 0) 


bietet jedoch charakteristische Schwierigkeiten und erfordert weitergehende 
Methoden, da die Theorie der Schwingungsgleichung nicht durch Ubertragung 
der Ergebnisse der Potentialtheorie zu gewinnen ist. So gilt beispielsweise nicht 
das Maximumprinzip; auch existiert kein Analogon zur Kelvintransformation, 
so daB die asymptotischen Gesetze im Unendlichen nicht aus dem Verhalten im 
Endlichen gewonnen werden kénnen. Andererseits interessieren gerade solche 
Loésungen der Schwingungsgleichung, die ein spezielles asymptotisches Verhalten 
zeigen, das durch die Sommerfeldsche Ausstrahlungsbedingung 


SEB U +o(—} fiir groBe 7 
beschrieben wird. 

Es zeigt sich, daB die Randwertaufgaben zur Helmholtzschen Schwingungs- 
gleichung unter der zusatzlichen Forderung dieser Ausstrahlungsbedingung ein- 
deutig lésbar sind. Fiir die erste Randwertaufgabe ist das von H. WEyLI [2] 
und C. MULLER [3] gezeigt worden. Die zweite und dritte Randwertaufgabe wird 
in dieser Arbeit durch Ubertragung des Beweisgedankens von C. MULLER be- 
handelt. Die beiden ersten Randwertaufgaben sind ferner von W. D. KuPRADSE [4] 
bearbeitet worden. Hier werden jedoch die Eigenwerte der auftretenden Integral- 
gleichungen nicht diskutiert, so daB zur Lésung wesentlich umfangreichere Me- 
thoden aus der Theorie der Integralgleichungen herangezogen werden miissen. 
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Wir betrachten in dieser Arbeit die Schwingungsgleichung fiir komplexe k 
mit Imk=0. Unsere Betrachtungen gelten sinngema4B auch im Falle Im k<0O, 
wenn wir dann fiir das Verhalten im Unendlichen die sogenannte Einstrahlungs- 


bedingung 
a = —ihU + 0(—| fiir groBe 7 
fordern. 

Im §1 der Arbeit werden die Eigenschaften der Potentialfunktionen ein- 
facher bzw. doppelter Flachenbelegung zusammengestellt. Der § 2 folgert die 
Eindeutigkeit unserer Aufgabe aus einem Ergebnis von F. Retiicu. Im § 3 
werden die Lésungen als Potentialfunktionen angesetzt und fiir ihre Belegungs- 
funktionen Integralgleichungen hergeleitet. Der § 4 diskutiert die Frage der 
Eigenlésungen. Im letzten Paragraphen schlieBlich wird der Existenzbeweis nach 


der Methode der Kapazitatskoeffizienten gegeben. 


§ 1. Vorbemerkungen 

Unseren Betrachtungen liege der dreidimensionale euklidische Raum zugrunde. 
Er enthalte das beschrankte Gebiet G; mit der geschlossenen Oberflache 2. Das 
AuBere dieses Gebietes bezeichnen wir mit G,. Die Flache Q sei orientierbar 
und glatt; wo es notwendig ist, die verschiedenen Seiten von 2 zu unterscheiden, 
geschieht das durch angehangene Indizes + (auBen) und — (innen). Der Nor- 
malenvektor von §2 weise nach auBen. Mit p bezeichnen wir einen beliebigen 
Punkt im Raum, wahrend s stets auf 2 liegen soll. 7(f, s) sei der Abstand zwischen 
f und s. Dann sind die Funktionen 


(1.1) U(p) = f{ w(s) as 
und ‘ 
(1.2) V(b) =f) 2 “as 


fiir alle CG, Lésungen der Helmholtzschen Schwingungsgleichung. Aus der 
Potentialtheorie tibernehmen wir fiir diese Funktionen folgende Ergebnisse [5]: 


Lemma 1. y(s) sez stetig. Dann ist U(p) iiberall stetig. Fiir die Ableitung in 
Richtung der Normalen gilt bei Anniéherung an Q die Sprungrelation 


- = — 27 u(s) Hyp als r) oa ae 

Ve =" 250'u(s) 4 fee) om ri ten as’. 
Lemma 2. 1. v(s) sez stetig. Dann gilt 

| ee Malan 28 Ot ae ers a as 

V_i=~+ 2av(s hia fats se me es 
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und es ist 


tkyr 
PrOz 7 ds! = 4 (V, +V.) 


any r(s, 8’) 
2 
holderstetig. 
2. v(s) set hdlderstetig. Dann ist 


tkyr 


0 , 0 é , 
S ed 
ds Jr ere 7(s,s°) a 
2 


vorhanden und hélderstetig. 

3. v(s) set hélderstetig differenzierbar. Dann existiert die Normalableitung von 
V(p) und ist beim Durchgang von p durch Q hélderstetig. 

Aus Lemma 1 und Lemma 2 erhalten wir 


Lemma 3. Dve stetige Funktion o(s) gentige einer Integralgleichung 
ikr 
o(S\=—A [ ols’ Le Ee 9 

2 


ony Y(s,s’) 


Dann existiert die Normalableitung der Funktion 
tkr 


ir nO Z , 
V(p) = [ o(s' Gaye 
Q 
und ist beim Durchgang von .p durch Q stetig. 


§ 2. Der Eindeutigkeitssatz 
In diesem Paragraphen beweisen wir 
Satz 1. Es gibt hochstens eine Funktion U(p) mit folgenden Eigenschaften: 
1. U(p) ist in G, zweimal stetig differenzierbar. 
2. U(p) geniigt in G, der Differentialgleichung 
AU +k?U =0 mit k== 0, lmk 0. 


3. U(p) nimmt an der Aufenseite von Q die stetigen Randwerte y(s) = ue + 
6(s)U(s) an. Die Funktion 6(s) ist dabei als reellwertig und stetig und im Falle 
Im k+0 kleiner oder gleich Null vorausgesetzt. 

4. U(p) geniigt gleichméfig fiir alle Richtungen der Ausstrahlungsbedingung 


aU . 1 is 
ae o(-] fiir groBe 7. 


Zum Beweis von Satz 1 benétigen wir folgendes Ergebnis von RELLICH [6]: 


Lemma 4. Die Funktion U(p) geniige fiir alle |p|=c der Gleichung AU + 
k2U =0. Es sei Lp eine Kugel um den Ursprung mit dem Radius R2c, und es 


gelte 
f\U}a@F =0(1) fir Ro. 
=R 


Dann verschwindet U identisch. 
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Es seien nun U und V zwei Lésungen unserer Randwertaufgabe. Dann 
erhalten wir durch W=U—V eine Funktion, die auf Q, den Wert 
ew =) 
annimmt; wir wollen zeigen, daB W identisch verschwindet. 


Dazu wahlen wir eine geniigend groBe Kugel 2X’ mit dem Radius R, die Q 
ganz im Inneren enthalt, und erhalten aus der Greenschen Formel wegen 


(2.2) f ( a W = )4F = 0, 
Q, = Ai 
(2.3) i (Ww AW—Waw)av = | (Ww = —W\aF, 
G 2, 


wenn G das von Q und » berandete Gebiet ist. Auf Grund der Ausstrahlungs- 
bedingung gilt nun mit k=k,+7k, 

Oa \({EeW . .tae 
(Fie W\(5,- + RW) 
_ |ew 
| én 


(2.4) 


ow 
on 


‘+ [RWP—ik(W Wa) +2hRe(W oe =) 


\ on RR 


Sodann ergibt sich aus der ersten Greenschen Formel 


on 
on Q 


(2.5) | wawav + f|Pwhav = [ w ov AF / - ow AF 
G G bs e 
oder 


(2.6) (—ke +3) [|wlav + {\7w 
G G 


2qV = Re [ We gee i WoOaR 
F on , on 
P2, 2 
Damit folgt aus (2.3) fiir groBe R 
: se see Pc Ce eee 
(2.7) }| ) dl +2h, | Wa dF =0(1), 


x G Q 


cow 
on 


RW|2+I7W 


+A) dF + 2h | ( 


Da nach Voraussetzung im Falle kg >0 


(2.8) ‘i Wears 
On a 
Q 


ist, folgt aus (2.7), daB W identisch verschwindet. Ist jedoch k,=0, dann erhalten 
wir aus (2.7) 


(2.9) {|W }P2aF =o(1). 


Nach Lemma 4 verschwindet W auch in diesem Falle, womit Satz 1 bewiesen ist. 


§3. Aufstellung der Integralgleichung 
Nachdem wir die Frage der Eindeutigkeit der Lésungen unserer Randwert- 
aufgaben behandelt haben, wollen wir uns nun mit dem Existenzproblem be- 
schaftigen und die Fragestellung zunichst auf Integralgleichungen zuriickfiihren. 


Neumannsches Problem fiir Schwingungen 105 


Dazu betrachten wir die Funktionen 
ikr 


(3.1) OU es OK N=. 9 


22 4¥(p,5) 
die fiir p+: s der Differentialgleichung 


(3.2) AU+RU=0 


geniigen. Die Funktionen Q und K sind Grundlésungen unseres Problems. Wir 
werden versuchen, aus ihnen durch Integration tiber 2 Lésungen der Rand- 
wertaufgaben aufzubauen, und bilden dazu mit stetigen Funktionen yu (s) und »(s) 


(3.3) U, (2) == i(s) Q(p, s) ds 
und ; 
(3.4) U; (2) a?) y(s) K(p,s) ds. 


Die Funktionen U,(f) und U,(p) geniigen in G, der Differentialgleichung (3.2) 
und der Ausstrahlungsbedingung. 


Zur Lésung unserer Randwertaufgaben machen wir den Ansatz U =U, und 
erhalten nach Lemma 1 


3.5) 2 +6U=y(s) =~pls) + f mls’) {KUs’, 8) +0(5) Qs’, s)}as' 
2 
wofiir wir symbolisch 


(3.6) y=u—-pwk 


schreiben. Die zugehdérige transponierte Integralgleichung bezeichnen wir ent- 
sprechend mit 


(3.7) y=u— Kg 


Damit haben wir unsere Randwertaufgaben auf eine Integralgleichung zuriick- 
gefiihrt, die wir nun diskutieren wollen. 


§ 4. Uber die Verteilung der Eigenwerte der Integralgleichung 


Die im vorigen Paragraphen aufgestellte Integralgleichung geht aus 


fiir A——1 hervor. Zur Losung von (3.6) untersuchen wir daher zunachst die 
homogene Gleichung fiir A = —1. 


Es sei ® bzw. H die Gesamtheit der Losungen der zueinander transponierten 
Integralgleichungen 


(4.2) y—Kyn=0; o—pK=0. 
Dann bilden wir mit einem 7¢€ H die Funktion 
(4.3) V(p) wae n(s') {K(B, s’) + 4(s’) Q(6, 8‘) as’, 


8* 
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die der Schwingungsgleichung sowie der Ausstrahlungsbedingung gentigt und 
zeigen zunachst, daB die Randwerte von @V_/én, die nach Lemma 3 ja existieren, 
zu ® gehoren. Es gilt namlich 


(4.4) V=—7y+Ky=0 
woraus nach Lemma 2 


folgt. Wenden wir nun die Greensche Formel mit den Funktionen V(g) und 
QO(p,9) (PEG,) auf das Gebiet G; an, dann folgt aus (4.4) 


— oV / / 
(4.6) o= [2 0@.s')as’ 
iia 
Durch Differentiation entsteht aus (4.6) bei Annaherung von # an Q, 
ce eVe Ee ; 
(4.7) 0=- 5 te K(s',s) ds’. 
Q 


Multiplizieren wir (4.6) fiir ~=s mit 6(s) und addieren das Ergebnis zu (4.7), 
dann folgt schlieBlich 


(4.8) ee + { {K(s’, s) + 6(s) O(s’, s)} ds’ 


Q 


womit die Behauptung bewiesen ist. Daraus ergibt sich 


. V_ : : ; 
Lemma 5. Sex My = a Dann ist die Abbildung gy =Mn umkehrbar ein- 
deutig. i 
Wegen 
OV naar Sas 
(4.9) rae 20n + os 


folgt namlich aus m =0 und (4.5) 
(4.10) a+ + OV, =0. 


Auf Grund des Eindeutigkeitssatzes erhalten wir daher 


(4.11) V=0 iniGes 
woraus sich wegen (4.5) 
(4.12) fea 


ergibt. Damit ist Lemma 5 bewiesen. Wir zeigen nun 


Lemma 6. Die Funktion U(p) geniige in G, der Differentialgleichung 
AU +kU =0 und der Ausstrahlungsbedingung. Die Grenzwerte von U und 0U [On 
ber Anndaherung an Q mégen existieren. Dann ist 

Im [us a= —CREk  miiG=6, 


nN 
24. 


wober C =0 nur dann gilt, wenn U(p) in G, identisch verschwindet. 
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Um Lemma 6 zu beweisen, wenden wir die Greensche Formel mit U und U 
auf das Gebiet G an, das von 2 und einer geniigend groBen Kugel » berandet 
wird. Es folgt 


443) f(C ee EA) a |p ie 7) ap 
R 


on 
ob, 


oder mit k=k,+7k, (k,=0) 


(4.14) Im [ U e dF = —2hyk, [ |U}aV —k, [ |UlPaF + 0(t) 
Du G Zp 


womit der erste Teil der Behauptung von Lemma 6 bewiesen ist. Aus 


(4.15) sf\U|?dF=o0(1) fir Reo 


folgt aber nach Lemma 4, daB U in G, identisch verschwindet. 

Setzen wir nun in Lemma 6 fiir U(p) die Funktion V(p) ein, dann erhalten 
wir wegen 
ney ee 
(4.16) V, =2n; a 


= Mn —26n. 


Lemma 7. Fiir alle yEH gilt 
Im fy-MyjdF =—CRek mit C20 
Q , 


wobet C =0 nur fiir 4 =0 moglich ist. 
Nach diesen Vorbereitungen beweisen wir 
. Satz 2. Es sei Imk+0. Dann ist A=—1 kein Eigenwert der Integralglet- 
chungen 
O=yptAypk 
O=e+AkKe. 
Zum Beweis betrachten wir die Funktion 
(4.17) U(p) = f e{K(p,s’) + 6(s’) Q(2,8')} as’ 
die auf 2 die Werte 


U,=etKe, U=—e+Ke 
. Us 0 
(4.18) Se lag 266 


annimmt. Aus (4.18) folgt 
(4.19) 4+a)U,—(1—¢a W=2(e+A1Ke). 


Es sei 4 nun reeller Eigenwert. Dann verschwindet (4.19), und es gibt zwei 
reelle Zahlen 

(4.20) cjA-1; fP=1A—A 

mit 

(4.21) aU — BU. =0: 
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Durch Anwendung der Greenschen Formel auf das Gebiet G; erhalten wir nun 


(4.22) [wa0—Uanar =f (v2 0 Far 
G; oe 
oder mit k =k, +7h, 
(4.23) Be gee: [\oPav =h,kD 
G; 


mit D0, wobei D=0O nur dann gilt, wenn U in G; identisch verschwindet. 
Zusammen mit Lemma 6 folgt daraus 


au, 
(4.24) —Im f {a¥, TUR ba =a hy CHBLE.D 
2 
woraus nach (4.18) und (4.21) 
(4.25) 2a fedU,dF =ak,C + Bk, k,D 
Q 
folgt. 


Es sei zunachst Re k =k,=- 0 und A= —1 Eigenwert. Dann ist ¢=7, und es 
folgt aus (4.25) D=0. Die Funktion U verschwindet also in G; identisch, woraus 
insbesondere 


(4.26) = 6) 


folgt. Nach Lemma 5 gilt dann aber ¢=0; in diesem Falle ist daher 2=—1 
kein Eigenwert. 

Sei nun Re k=k,=0. Dann folgt durch Anwendung der ersten Greenschen 
Formel 


(4.27) fous ue Mel Laas 
Q- 

oder 

(4.28) BR RS D,=0. 


Andererseits gilt 

(4.29) ) Lon lee A hy [ |UPaF + (1) 
ZR 

oder 


(4.30) Sats a mit C,=0. 


Wir erhalten daher an an Stelle von (4.24) 


(4.31) slit aU, ae Spl oe  }aF = aC,+BD, 
oder 


(4.32) 20 [6 0UdP =a Ge 8 De 
Q 
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Angenommen 4=—1 ware Eigenwert. Dann verschwindet D, und nach (4.27) 
U in G;, was wieder zum Widerspruch fiihrt. Damit ist Satz 2 bewiesen. 
In Satz 2 haben wir eine Aussage tiber die Verteilung der Eigenwerte unserer 


Integralgleichung im Falle Imk=E0 erhalten. Wir beschaftigen uns nun mit 
reellen k und beweisen 


Satz 3. Set k=-0 reell. Dann sind die reellen Eigenwerte der Integralgleichungen 
O=yp+AykK 
O=e+AKe 
einfach und entweder gleich +1 oder —1. 
Zum Beweis betrachten wir die Funktion 


(4.33) U(p) ee Q (2, s') ds’ 
die auf 2 die Werte 
ou, . aU_ 
(4.54) tt a yp ky; em +p ky 
4.34 
Ol RO pk, Svs) K(s’, s) ds’ 
annimmt. Aus (4.34) folgt 
ou, oU_ 
(4.35) (1 — A) oe (4+ A) ae =2(y +AykK) —2A6U. 
Es sei A nun reeller Eigenwert. Dann gibt es zwei reelle Zahlen 
(4.36) Cat baie p= t—A 
mit 
aU_ Bes 
(4.37) a Ui as pee aay AU, 
Aus (4.37) folgt analog (4.21) ff. 
aU_ av, re. 
(4.38) Im [ Ufa ———B aa} oF Bh,C 
oder 2 
(4.39) O=BRC=1—A/) RC. 


Angenommen U ist in G, von Null verschieden. Dann folgt A=1. Wenn U 
jedoch in G, identisch verschwindet, ergibt sich aus (4.37) 


Oe 
(4.40) (4 + A) ors O, 
das heiBt, es ist A = —1 oder U verschwindet iiberall identisch; dann verschwindet 


aber auch yp. A= -+1 sind also die einzigen reellen Eigenwerte. Wir wollen zeigen, 
daB sie auch einfach sind. 

Dazu nehmen wir an, A ware mehrfacher Eigenwert. Dann gibt es Funktionen 
y* und ¢%, die nicht identisch verschwinden mit [7] 


ptrAyK=0; pwtApypK =k 
e tAKe=0; e? 1) Ke” = K et. 
1 Definiert in Gl. (4.41) 


(4.41) 
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Wir wollen (4.41) fiir 2=-1 zum Widerspruch fiihren. Dazu wahlen wir die 
Funktionen 
(4.42) Oe a ara 
2 
und 
(4.43) Usa =f 8" {K(p, 8’) + 6(s') Q(6,s')} as". 
Q 
Nach (4.19) und (4.35) folgt 
(1—1) 2G 4 (4 4 2) 2 = — 20 
(4.44) ka. eee 
(1 — A) ae t (1 + A) ioe = ae 26(A U2 —U4;) 
und ao 
Ae AN Onl doe -=0 
(4.45) (1 + A) Us — ( ) Us, 


(4 + A) Ue,o+ — (4 —A\U aie ee Us3+ +U 21. 
Wahlen wir zunichst 2 =1. Dann folgt aus (4.44) 


ou, 
a =— dU, 
(4.46) Aen AU | 
2 ee = ae 26U,,+ 6U yj. 


Multiplizieren wir die erste dieser Gleichungen mit 2U,,, gehen wir in der zweiten 
zur konjugiert komplexen iiber und multiplizieren sie dann mit — U;,+, so folgt 
durch Addition 

aU, » CUT 


Hh He a z Me 
(4.47) 2 {pu Ty, — te Uy = — Fo Uy, — by Di. 


Integrieren wir (4.47) ttber 2, dann verschwindet die linke Seite des Integrals 
nach der Greenschen Formel, und wir erhalten 


(4.48) 0 =Im [ Ue oes 
on 
Qy 
Nach Lemma 6 ist das nur méglich, wenn U,, in G, identisch verschwindet. 
Dann folgt aber aus (4.46) 


(4.49) hes dik au, ,- 
on on 
Nach Lemma 1 ist daher entgegen der Annahme »,=0. 
Sei nun A= —1. Dann erhalten wir aus (4.45) 
(4.50) Uai- = 03) = 2U yy = Une 


woraus jetzt 
(4.54) 2 {Use 
folgt. Nach (4.18) gilt aber 

(4.52) Up + = 2e1, 
so daB aus (4.51) nach Integration tiber Q 


Uy 9 


ae Wie ae Uae —= OU 7 


(4.53) pao Uy, 20 dF 


2, 
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folgt. Gl. (4.53) ist wieder nur dann méglich, wenn U,, in G, identisch ver- 
schwindet. Dann gilt aber nach (4.50) 


(4.54) Up U21- = 2e'— 0, 


so daB e, entgegen der Annahme verschwindet. Damit ist auch Satz 3 bewiesen. 


§5. Auflésung der Integralgleichung und Beweis des Existenzsatzes 

Wir haben im letzten Paragraphen Aussagen iiber die Verteilung der Eigen- 
werte unserer Integralgleichung gewonnen und sind nach diesen Vorbereitungen 
in der Lage, den Existenzsatz zu unserem Randwertproblem zu beweisen. Wir 
zeigen 

Satz 4. Es gibt eine Funktion U(p), die auf Q die Randwerte y a +6U 
annimmt und im ubrigen den Voraussetzungen von Satz 1 geniigt. 

Im § 3 haben wir diese Randwertaufgabe auf die Lésung der Integralgleichung 
(5.1) y= — wk 
zuriickgefiihrt. Wie wir in Satz 2 sahen, besitzt sie fiir _Im(k)=+-0 keine Eigen- 
lé6sungen und ist in diesem Falle also stets lésbar [7]; wir brauchen uns daher 
nur noch mit dem Fall reeller k (k =k) zu beschaftigen. 

Wir setzen also fortan k reell (k =k,) voraus und miissen in diesem Fall die 
Eigenlésungen von (5.1) diskutieren. Dazu zeigen wir 


Lemma 8. Es sez k reell. Dann gehort mit y auch g zu ®. 


Der Beweis folgt analog zu (4.4—8), wenn wir statt V die Funktion V ver- 


wenden. 
Es sei nun 7’ eine Basis von H. Dann erhalten wir nach Lemma 5 durch 


y'=M7' eine Basis von ® und beweisen 


Ce) eat 7) 


ist symmetrisch und nicht entartet. 
Denn ware det C’*=0, so gabe es ein nichttriviales 7*€H, derart, daB fiir 
alle pE@® 


Lemma 9. Die Matrix 


(5.2) (p,7*) =0 
galte. Da nach Lemma 8 auch M7* zu ® gehért, folgte 
(5.3) (n*,M7*) =0. 


Das ist nach Lemma 7 nicht méglich. Um die Symmetrie der Matrix GF 7 
zeigen, setzen wir 


(5.4) V' (6) ap ie s) + 5(s) Q(p,s)}ds. 
An der AuBenseite von 2 gilt dann 

(5.5) Vi= 2m; 

(5.6) Y= Mri — 267, 
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und durch Anwendung der Greenschen Formel folgt 


{vs ye— yr 2-V'\ aF =/ (Vs vt ve o_V') aF 
x on on \ on 
(5.7) 22 a, 


= 2(n', Mn*) — 2(n*, My’) =2(C**— CC). 


Wegen der Ausstrahlungsbedingung 


Lene 1 
(5.8) gp Vixikv Hole 
und 

ae 
(5.9) Vi =0() 
gilt jedoch fiir groBe R 
iO WR nr 8 Wy = 

(5.10) {v 2_V*— Vb2_V\dF =olt). 

=R 
Damit ergibt sich aus (5.7) 
(52441) CHE CU — O04) tar Jehu sce 


und Lemma 9 ist bewiesen. 


Diese Vorbereitungen erméglichen es uns, Satz 4 auch im Falle reeller k zu 
beweisen. Dazu bestimmen wir zunachst zu y die durch die Bedingungen 


(5312) (y7,y—n*)=0 fir alleneH 


eindeutig festgelegte Funktion 4*€H. Da die Matrix C** nicht entartet ist, 
k6énnen wir zu 7* eindeutig ein g*c ® mit 


(5.13) (7,n* — ot) =0 ftir allene€H 

bestimmen. Fiir dieses m* gilt dann auch 

(5.14) (7,7 —@*) = 0" — fir valleaeh: 

Die Integralgleichung 

(5.15) u—-ukK=y—g*; O0=(u% 9) fir allepe®@ 

ist daher eindeutig nach w auflésbar. Bestimmen wir noch 7’ aus der Gleichung 
(5.16) g* =M7', 

dann erhalten wir die gesuchte Lésung unserer Randwertaufgabe in der Form 
(5.47) U(p) = J #(s) Q(p,s) ds ae n'(s) {K(p, s) + 6(s) Q(6,s)}ds. 


Es gilt namlich an der -AuBenseite von 2 wegen (4.16) 


(5.18) su 4+ 20 =—pw+ukK + M7 =y, 


womit die zweite und die dritte Randwertaufgabe des AuBenraumes fiir die 
Helmholtzsche Schwingungsgleichung gelést sind. 
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Die Anwendung des Operatorenkalkiils 
von Jan Mikusitski auf lineare Integralgleichungen 
vom Faltungstypus 
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Vorgelegt von A. ERDELYI 


1. Orientierung und Problemstellung 


Die Operatorenrechnung, die JAN MrkusINskI [10], [17] eingefiihrt hat, erlaubt 
Anfangs- wie Randwertprobleme bei gewéhnlichen und partiellen Differential- 
gleichungen mit konstanten Koeffizienten ganz einfach kalktilmaBig oder sym- 
bolisch zu lésen. Nur solche Differentialgleichungen werden betrachtet, deren 
Lésungen in einer geschlossenen Form zu erhalten sind. Bisher sind verschiedene 
Methoden bekannt, um diese Gleichungen zu lésen, wie z.B. die Methode von 
HEAVISIDE, die der Integraltransformationen (LAPLACE, FOURIER, HANKEL, 
MELLIN usw. [3], [23], [20], [21], [25]), die Methode, die sich auf die Theorie der 
analytischen Funktionalen von L. FANTAPPIE stiitzt, und in neuester Zeit der 
Distributionenkalktil von L. ScHwarvz ({J/a], [17], [19]). 


Bekanntlich laBt sich die Operatorenrechnung im elementaren Fall der gewohn- 
lichen Differentialgleichungen auf algebraischem Wege (ohne Laplace-Transfor- 
mationen) begriinden ([3] Bd. II). Mit anderen Worten: in diesem Fall laBt sich 
der Heaviside-Kalkiil in einem gewissen Sinne mathematisch rechtfertigen. Dage- 
gen war im transzendenten Fall der partiellen Differentialgleichungen bis vor 
kurzem die Laplace-Transformation das einzige Mittel, um dieses Problem exakt 
kalkiilmadPig zu losen. Nun ist einer der vielen Vorteile des Mikusinski-Kalkiils 
die exakte selbstandige Fundierung des letzten Problems auf algebraischem oder 
symbolischem Wege. 

Auch ist die Theorie von MixusiXsk1 allgemeiner und einfacher als die auf 
der Laplace-Transformation beruhenden Darstellung, da z.B. das Verhalten der 
gesuchten Funktion im Unendlichen bei der Laplace-Transformation eine Be- 
schrankung auferlegt, die zu gewissen Schwierigkeiten bei den Eindeutigkeits- 
satzen fiihrt. Sogar wird die Kenntnis der Funktionentheorie nicht vorausgesetzt. 

Es ist zu erwahnen, da der Verfasser [2] dieses Problems mittels der Laplace- 


Transformationen vom Standpunkt der Funktionalanalysis und insbesondere 
der Halbgruppentheorie betrachtet hat. 


Die vorliegenden Resultate wurden am 10. April 1958 in einem Vortrag anla@lich 
der Gamm-Tagung an der Universitat des Saarlandes mitgeteilt. 
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Im Vergleich mit der Theorie von Fantraprié! oder ScHwartz? braucht sie 
nicht den groBen topologischen und funktionalanalytischen Aufwand. Die Er- 
weiterung des Begriffs der natiirlichen Zahlen zu den rationalen und den reellen 
Zahlen ist ein Modell fiir die Methode von MixusrNsk1. Kurz gesagt: die stetigen 
Funktionen einer reellen Variablen ¢, 0, bilden bei der gewohnlichen Addition 
und der Faltung als Multiplikation eine kommutative Algebra, die zu einem 
Quotientenkérper K erweitert werden kann (da diese Algebra Null-Teiler frei ist), 
dessen Elemente dann Operatoren genannt werden. Dieser Operatorenkérper K 
enthalt die komplexen Zahlen, die stetigen Funktionen, den Differential- und 
Integraloperator. Der Dirac-Delta-Operator ist das Einheitselement dieses 
Korpers K. 


Um iiber die Algebra der Operatoren hinaus zu einem Infinitesimalkalkiil 
vorzustoBen, wird zunachst in der kommutativen Algebra ein Konvergenzbegriff 
eingefiihrt und dieser dann in geeigneter Weise auf Folgen von Operatoren tiber- 
tragen. 


In dieser Arbeit wird die Methode von MixusINski auf lineare Volterrasche 
Integralgleichungen vom Faltungstypus von erster und zweiter Art angewandt. 
Auch werden Systeme von linearen Integralgleichungen sowie Integrodifferen- 
tialgleichungen und singulare Integralgleichungen zusammen mit verschiedenen 
Beispielen betrachtet. Uber nicht-lineare Integralgleichungen l4Bt sich auch 
etwas aussagen. Integralgleichungen, die ,,verallgemeinerte’’ Funktionen wie 
z.B. die Diracsche Deltafunktion enthalten, lassen sich mittels dieser Theorie 
auflédsen. Diese kommen in der Theorie der elektrischen Netze (z.B. der Kurz- 
schluBstrom) und des Balkens vor. Zum Beispiel lassen sich die Einzelkrafte 
und Einzelmomente in der Theorie des Balkens nicht durch gewohnliche Funk- 
tionen wiedergeben, kénnen jedoch durch Operatoren dargestellt werden. In § 5 
wird ein anderes Beispiel einer Integralgleichung angegeben, welches sich auch 
nicht mittels der Laplace-Transformation lésen laBt, weil das entsprechende 
Laplace-Integral divergiert. 

Diese Arbeit ist absichtlich einfach abgefaBt, so wie es in der Operatoren- 
rechnung wblich und wiinschenswert ist. 


2. Vorkenntnisse aus der Operatorentheorie von Mikusinski 


Es sei C die Klasse derjenigen Funktionen, die im Intervall 0St<-+ oc 
definiert und stetig sind. Wir unterscheiden zwischen der Funktion f ={f (t)} 
und dem Wert f(t) der Funktion im Punkte ¢ (eine Zahl). Wir definieren die 


1 Eine zusammenfassende Einfiihrung in diese Theorie mit Anwendungen auf 
Differentialgleichungen findet man in F. PELLEGRINO [12] oder [13]. Jedoch ist zu 
bemerken, daB die topologische Grundlegung der Theorie von L. FANTAPPIE, die in 
verschiedener Hinsicht recht unbefriedigend war, erst in jiingster Zeit in der Habili- 
tationsschrift von H. G. TILLMANN [24] dargelegt wurde. Fiir den Zusammenhang 
zwischen dieser Theorie und der auf der Laplace-Transformation beruhenden Methode, 
siehe auch [7]. 

2 J. KorEvAaAR [7] gibt vom Standpunkt der angewandten Mathematik aus eine 
Einfiihrung in diese Theorie, die sehr rasch zu den Anwendungen fiihrt. Siehe auch 
J. MrxusiNsx1 und R. Sikorski [//a], sowie E Siva [19]. Eine algebraische Be- 
griindung der Distributionentheorie hat H. K6nic [6] angegeben. 
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Addition und das Produkt (Faltungsprodukt) zweier Funktionen f und g durch 
(21) ft+e=fO +s} dew. fe=(O}{gO} = {fre — u) g(u) du; 


das Skalarprodukt durch « {f(#)} ={a/()} wo eine Zahl (reell oder komplex) 
ist. Dies fiihrt zu einer kommutativen Algebra (ohne Einheitselement), die nach 
einem Satz von TrrcumarsH ([23] S.328 oder [11] S.9—17) Null-Teiler frei ist. 
In anderen Worten ist eine Erweiterung zu einer kommutativen Divisionsalgebra, 
der sogenannte Operatorenkérper K, méglich. Die Elemente dieses Korpers, die 
Operatoren, sind Briiche //g fiir 7, g€C, g= {0} (natiirlich versteht man hierbei 
nicht die gewohnliche Division, sondern die zur Faltung inverse Operation) und 
die Rechenregeln entsprechen denen der rationalen Zahlen. Dieser Kérper K 
enthalt die Funktionen / € C (weil sie in der Gestalt /g/g, wo g+- {0} eine beliebige 
Funktion aus C ist, darstellbar sind) wahrend nicht jeder Operator eine Funk- 
tion darstellt. So ist ein Operator eine Verallgemeinerung des Begriffs einer 
Funktion. Auch enthalt K die komplexen Zahlen (die in der Gestalt «=«g/g, 
g++ {0} darstellbar sind). So wird der Begriff der komplexen Zahl dem Begriff 
des Operators untergeordnet. Die Funktion /={1} mit der Eigenschaft / {f (¢)} = 


, t 
{ fils) a u} wird Integrationsoperator genannt. Im allgemeinen gilt fiirm=1, 2,... 
0 


j* — fe—at 

(m—1)!J’ 
und die Faltung von {/(/)} mit /" bedeutet n-malige Integration von {f(¢)} tiber 
das Integral von 0 bis ¢. Wieder gehért /” zu K (sogar zu C). Bedeutet L die 
Klasse derjenigen Funktionen {7(d)} die lokal summierbar sind (im Sinne von 
LEBESGUE), d.h. die auf jedem Intervall O<¢<¢, summierbar sind, so kann 
auch die Klasse L in den Kérper K eingebettet werden (in dem man / mit 


t 
+ Jt (w) du} als Quotient zweier Funktionen aus C, identifiziert). Hierbei ist 
0 


4 das Einheitselement 7 von K, welches der Diracschen ,,Deltafunktion 6 ent- 
spricht*. Auch enthalt K den Differentialoperator s =1/I. 


Ist {f(t)} eine total-stetige® Funktion fiir t= 0, so gilt (| 20).S.47, [LLies. 103) 
(2.2) sf=f'+ (0), 


wobei f(0) der Wert der Funktion f im Punkte ¢=0 ist (ein Zahlenoperator) 
und /’ ={/'(0} die Ableitung von f ist. Ist /(0)=0, so bedeutet sf=/’ einfache 
Differentiation. Aber das Produkt sf hat stets einen Sinn (sogar wenn / nicht 
differenzierbar ist). Allgemeiner, ist die (m—1)-te Ableitung /"-) total-stetig 


® Der Einfachheit halber, kénnten wir annehmen, daB {/(¢)} nur endlich viele 
eee Sean eo in jedem endlichen Intervall besitzt und dort absolut integrier- 

ar ist. 

4 Es ist zu bemerken, daB vor kurzem C. ScHMIEDEN und D. Lauewirz [16] auch 
eine Theorie dargelegt haben, die u.a. zu einer Erweiterung des Funktionsbegriffes 
fiihrt und so auch die Delta-Funktion enthalt. 

> Einfachheitshalber kénnten wir die starkere Voraussetzung f/€C, wobei die Ab- 
leitung /’ zur Klasse gehdrt, die in FuBnote 3 erklart ist, benutzen. 
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und /® = {7 (} so gilt 
(23) sf =fOF 9 (0) + sf*-2(0) +--+ + 524 (0) 
wo /(0) den Wert von /* im Punkte t=0 bedeutet. 

Als Anwendung dieser Uberlegungen, erhalt man z.B. 


(2.4) Ne me oe {sin at} = | a za} 


Ist f={f(A)}€L, so hat Cz. Rvtt-Narpzewsk1 ([10] S. 62, [11] S. 322) gezeigt, 
daB 


CO 


(2.5) fes*Fa)da =(f(O} 


0 


gilt, wobei dieses Integral als ein abstraktes Integral (ein Operator) im Sinne 
von MixusiNsxt ([10] S. 61, [11] S. 306) zu verstehen ist (hier bedeutet e~°?, 
A>0, den Verschiebungsoperator). Ist {7 (2) (EL und betrachtet man s als eine 
komplexe Veranderliche (nicht als den Differentialoperator) so kann man das 
obige Integral als ein gewohnliches Integral auffassen und als eine Transformation 
die jeder Funktion /€ LZ der reellen Veranderlichen t= 0 (fiir welches das Integral 
konvergiert) eine analytische Funktion der Veranderlichen s zuordnet. So haben 
wir eine formale Beziehung zwischen der Laplace-Transformation und den 
Operatoren. Aber beide Methoden sind nicht aquivalent, da die Laplace-Trans- 
formation nur anwendbar ist, falls das Integral (2.5) konvergent ist. Dagegen 
hat das Integral (2.5) im Sinne von MIKUSINSKI immer einen Sinn falls nur /€ LE. 
So ist es auch méglich die bekannten Tabellen zur Laplace-Transformation zu 
verwenden. 

Eine Folge von Funktionen /,, /;,...¢ C, ist C-konvergent gegen f€C, oder 
C-lim f/,=/ genau dann, wenn /,(t)—/(i), noo, gleichmaBig in jedem abge- 
schlossenen Intervall OStS%,. Eine Folge von Operatoren /,,/.,...€ K ist 
K-konvergent gegen {C K, K-lim /,=/, dann und nur dann wenn es ein Element 
qe C, ¢+{0} gibt, so daB gf, C-konvergent gegen q/ fiir »— co ist. Zum Bei- 
svete K-lim ae ea == 4, 


n> CO 


welches eine Approximation des Einheitsoperators von K durch eine Folge 


gewohnlicher Funktionen aus € darstellt®. Man sagt eine Reihe }'/, sei gegen h 
n=1 


konvergent, falls bm a Ar fot-:-+/7,) =h gilt (wobei beide Konvergenzarten 
moglich sind). 
3. Drei Sadtze iiber Potenzreihen von Operatoren 
Bekanntlich gilt ({10] S. 64, [11] S. 155, [5] S. 39) 
Satz 3.1. Ist FEC, so folgt 


CO 
mn if it 2 
= wo = Usw. 
(3-1) Diao ght athe). 
6 Eine entsprechende Eigenschaft gilt in der Distributionentheorie, z.B. [/8}, 


Livre II S. 29. Der Zusammenhang zwischen dieser Theorie und der von MIkusINsKI 
ist z.B. in TEMPLE [22] betrachtet. 
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Vollstandigkeitshalber skizzieren wir den Beweis. Wahle ein beliebiges festes 
fy>0. Da fC ist, existiert ein M mit |/@)|SM fiir OS#St. Wir haben 
unter der Bezeichnung jee AS A eee 
t 
AfOH=10, PO =Sie—w ho) du 
0 
und durch Induktion gilt 


= pS 
POLS Goa Tate) Ostet. 


co 


Hieraus folgt die gleichmaBige Konvergenz der Reihe x /*(@) in OStS4, die 
n=1 foe) 
dort eine stetige Funktion darstellt. Da aber 4) beliebig ist, ist die Reihe >) /” 


n=1 
im Sinne der C-Konvergenz zu einem Grenzwert CC konvergent. Um h zu 


finden setzen wir Shei S0-gilt 
(k= U-A P= LP- LPs, 


ake n=1 n 


woraus h=}/(1—/) folgt, w.z.b.w. 
Im Falle, daB f quadratisch integrierbar ist (oder im Sinne von LEBESGUE 
daB fC L? ist) beweisen wir nun 


Satz 3.2. Gilt f€ L’, so folgt 


(3.2) [ie one air A ine WEL be. 


n=1 


Beweis. Wir nehmen irgendein ¢,>0 und halten es fest. Da {€ L? ist, setzen 
wir 


AC) du = M,, fiftalean ae) 


Die folgenden Ungleichungen gelten alle auf dem Intervall O<t<¢. Zuniachst 
gilt |/f()|SM,. Nach der Cauchy-Schwarzschen Ungleichung gilt 


7] = sie — u) f (uw) du < [Sin — wlan {f@)[du) My, 


woraus |/ f?(t)| <M, ¢/1! folgt. 
Auch gilt 


LO) sn — u) Pu) du| <M, f | (¢ —u)|\du=M,M,, 


und so ist |/ /#(¢)| <M, M,#/1!. 
Ferner ist 


(| =|PPO| SMBS du = Mey, 


t 
und so gilt |/ /*(@)|< My f udu =M? 2/21. 
0 
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Auch gilt F 
1@| =|PPO|S MEM, f du = M2 M, 1/11, 
0 
woraus |/ f° (#)| < My M, #/2! folgt. 
Im allgemeinen (durch Induktion) gilt 


[POLS Mpe-y@—1)!, [PPO] SMe en —1)! 
mit 
[2P"oQ|SMprint, lle) <MeM,eyn!. 
Nun ist 5 if" (t) eine (gewdhnliche) unendliche Reihe von stetigen Funktionen 


die im Tatervall OStS7%, durch die Reihe 
Sage = (1+ M,) oe 


majorisiert wird, und also gleichmaBig konvergent in O<#S@, ist. Also stellt 


die Reihe >'//"(é) eine stetige ee dar (alle Glieder sind stetig). Folglich 


n=1 © 
existiert C-l >’ f” und daher K- s f/f". Um den K-Grenzwert zu finden, ees gs 
n=1 n=1 
man dieselbe Methode wie im vorhergehenden Beweis. Auch ist ¥ Pec. 


(oe) n=2 
Infolgedessen ist sowohl h=)'/”C€ HE als auch €L? da 


n=1 


s|/e+|ardr|+ 


Damnit ist dieser Satz bewiesen. 


Nun betrachten wir auf dem Intervall [0, oo) definierte Funktionca, die nicht 
dort beschrankt sind, aber welche eine Singularitat besitzen, so daB | / (¢)|< M//t*, 
O<a<1 gilt. 


Satz 3.3. Gilt |f()|<M/t, 0<a<1, so folgt 


paras | 
ys a = 
n=1 
Beweis. Fir 1” (t) gilt die folgende Abschatzung: 


MIA a) Bao 
(3.3) WQS eaceay | PE 


Offenbar ist diese Ungleichung giiltig fiir 7 =1. Der Beweis wird mittels Induk- 
tion gefiihrt. Zunachst gilt 

n( [MI — aj" @a-9—1 
(3.4) I” Ol Ss 


=f yeaa 
I((14 — a) 4 * 


(wo I'(«) die Gamma-Funktion bedeutet). Um (3.4) zu beweisen gilt 


Sao Mr esi N pnN—na—1 
PeilghGlayae resus oe 
0 
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und setzen wir u —tv in das Integral ein, so gilt 
1 
Panta fu — y) ty" "AT dy 


I'(n—n«) 
0 
MS M*t17"(1—«a) eee ida: ta eee (n— no) 
T(n—na) TD (n—na+1—2) 
ee =a) rt Unie “—1 


T'((n+1) )(1—a)) } 


woraus der Beweis von (3.4) und (3.3) sofort folgt. (Fiir ahnliche Abschatzungen 
siehe MIKHLIN [9] S. 18.) Wir nehmen irgendein festes fj >0. Um zu beweisen, 


daB die Reihe J /"(¢) gleichmaBig und absolut auf dem Intervall 0OS/S% 
n=1 


konvergiert, geniigt es zu zeigen, daB die Reihe, die a, als allgemeines Glied hat, 
konvergent ist. Hier ist a, das Glied auf der rechten Seite von (3.3), wobei ¢ 
durch f, ersetzt ist. Hierzu braucht man nach Caucuy’s Kriterium’ nur zu zeigen, 
daB lim Van<1 ist. 

n— CO 

In der Tat, mittels der so eats ro gilt 


Pie —— =; Bt exp(— B+ a) O<b<1. 


Hieraus folgt 


eel 4 ETS) 
n ~ T3 (n(1—a) +1) 


M t-*I'(1—«) [2 (1—a) +22)!" exp( 1 a4 : atta) 


[m (4 —«x) 44 ]i-a+in 
und da dieser Ausdruck mit 1/m gegen Null strebt, ist der Satz bewiesen. 
Diese drei Satze bilden die theoretische Grundlage zu den folgenden Betrach- 
tungen. 
4. Integralgleichungen erster Art 
Die Volterrasche Integralgleichung erster Art hat die Form 


(4.1) The May aioe 


wo fiir #20, k und g gegebene Funktionen sind und / gesucht ist. Hier hangt 
der Kern k(¢—w) nur von der Differenz t—wu ab. Es sei k, g, fC L. Die Opera- 
torenmethode fiihrt die Gl. (4.1) in eine gewhnliche lineare algebraische Gleichung 


(4.2) kf=g 


iiber, die im Operatorenkérper K die Lésung f/=g/k hat. AnschlieBend findet 
man die Lésung im Originalraum (falls diese existiert). Wenn die Integral- 
gleichung (4.1) eine Lésung hat, so ist sie eindeutig bestimmt. Hatten wir nam- 
lich zwei Lésungen /, und /,, so ware kf, =g und kf,=g, woraus kh (Gas) ={ot 
folgt. Der Satz von TitcHMarsH ergibt },=/, oder k= {o} was ee ware. 


aa, B Le 
dass ae eispie Knopp, Theory and application of infinite series, S. 116. Lon- 
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Ein erstes Beispiel fiir die Anwendung dieser Methode wird geliefert durch 
die Gleichung 
t 
(4.3) J sin(t — u) f(u) du =e(), 
¥ 0 
wo g gegeben ist, f unbekannt, ‘20 ist. Nach (2.4) § 2 nimmt die Gl. (4.3) die 
Form 
1 


vi! 8 


im Operatorenk6érper mit der Lésung f=g-+s?¢ an. 


Wir haben also eine ,,verallgemeinerte Lésung“‘ der Integralgleichung erhalten, 
d.h. eine Lésung, die ein Operator im Sinne von MixusINskI ist 8. Diese Lésung 
1a8t sich auf eine gewohnliche Funktion reduzieren, wenn s?g eine gewohnliche 
Funktion ist. Ist die Ableitung g’ total-stetig fiir ‘= 0, so ergibt die Gl. (2.3), daB 


s?g =g" + g'(0) +s¢(0) 


ist. Ist nun auch g(0) =g’(0) =0 so reduziert sich s*g auf die Funktion g’’ und 
wir haben als eindeutige Lésung von (4.3), /(¢) =g(¢) +-g’’(¢). In der Tat ist dies 
eine Lésung von (4.3) wie die Substitution von f(t) in (4.3) zeigt. 

Ein zweites Beispiel liefert die entsprechende Gleichung mit dem Kern 


(4.4) k(t) = 


Die Operatorenform dieser Gleichung ist /”f=g mit Lésung f=s"g. Diese 
Losung existiert immer als ein Operator falls g nur €E ist. Ist aber auch die 
(1 —1)-te Ableitung g~» total-stetig und g(0) =g’(0) =---=g“—)(0) =0, so 
hat die Integralgleichung (4.1) mit dem Kern (4.4) die eindeutige Loésung 
/(t) =g™(t) was eine gewohnliche Funktion ist (s. SCHMEIDLER [15], S. 210 fiir 
ein ahnliches Beispiel). 


Die Integralgleichung 
t 
SIo(t—u)t(udu=e), t>0 
0 


wo Jy(t) die Besselsche Funktion vom Index 0 ist, hat die eindeutige Lésung® 
t 
y= , 
HO) = {| 2M eu) du +e’ 
0 
falls g’ total-stetig ist mit g(0) =0. In der Tat, der Ubergang zur Operatoren- 


form ergibt 
(oO}FO} =8- 


8 Die klassischen Methoden liefern keine ,,verallgemeinerte Losungen*’. Obwohl 
diese keine gew6hnlichen Funktionen sind, ist es jedoch 6fters méglich, sie zu inter- 
pretieren, wie es z.B. bei der Diracschen Delta-Funktion méglich ist. Der Distribu- 
tionenkalkiil von Scuwarrz [17] liefert auch Losungen, die Distributionen sind, 
welche sich in gewissen Fallen auf Funktionen reduzieren. 


9 Siehe auch: ScHwartz [18] Livre VI, Transformations de Laplace, S. 21. 
Q* 
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1 (2 
Da sv aa yep ae ae 


({11] S. 163) ist, gilt 
ee {hO\e +s¢, 


und der Beweis folgt wie in den obigen Beispielen. 

Ist A(t) und g(t) n-mal differenzierbar (n=1) fir 20 und k(0) =2'(0) = 
. = K"-2)(0) =0 und #"-Y(0) £0, so ergibt sich (s. auch [3] Bd. III, S. 154) 
durch n-malige Differentiation von (4.1) falls eine fiir ¢>0 stetige Lésung f(t) 
existiert, 


J WOO (E — 1) fee) de + RO-(0) F() = CW 


fiir <>0. Man erhalt nach Division durch k%~)(0) hieraus eine Integral- 
gleichung von zweiter Art, die nach der Methode von §5 behandelt wird. Die 
obige Methode versagt z.B., wenn k(¢) in t=O keine Ableitungen besitzt, wie in 
dem Fall k(t) =¢~*, «>0, was uns zu den singularen Integralgleichungen fihrt, 
die in §8 betrachtet werden. 


5. Integralgleichungen zweiter Art 
Diese hat die Gestalt 


(5.1) fre —w) fu) du +80) =i 6 


wo wieder 7 die gesuchte Funktion ist und k, g¢€ Z bekannt sind. Diese Gleichung 
nimmt in der Operatorenform die Gestalt 


kf+g= 


an, mit Losung 

(5.2) Pei ks 

Die Lésung der Gl. (5.1) ist wiederum eindeutig bestimmt falls sie existiert. 
Gabe es zwei Losungen, so wiirde ihre Differenz die Gleichung k/,=/, befriedigen, 
so daB fs = 10) oder k =1 ware, was beidemal trivial ist. Da in der Gestalt (5.2) 


die Lésung meistens nicht unmittelbar! zum Originalraum zuriickzufiihren ist, 
schreiben wir die Lésung in der Form 


[ies Sir eee c 


und entwickeln k/(1—k) in eine Potenzreihe, indem wir einen der Satze aus § 3 
benutzen, je nachdem kE€C oder k€ LE? ist. 


Satz 5.1. Es set f,g€L und REC. So hat die Integralgleichung (5.1) die 
Lésung 


(5.3) {) =¢@ + J hla) gi ae, 


*° Dagegen erhalt man unmittelbar /(t)=¢ als Lésung der Gl. (5.4) mit & (js 
(0) Sint CoO, 
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wobet 


eine Rethe ist, die in jedem endlichen Intervall OStSt, absolut und gleichmafig 
konvergent ist und eine stetige Funktion darstellt. 


Der Beweis folgt sofort aus Satz 3.1 und es laBt sich zeigen, daB diese Lésung 
die Gl. (5.1) erfiillt. 


Die Reihendarstellung von / ist in der allgemeinen Theorie der linearen 
Integralgleichungen als die ,,.Neumannsche Reihe‘“ und die hk” sind als die 
,iterierten Kerne‘‘ bekannt. 


Auch gilt 


Satz 5.2. Wenn g und k€ L? sind, so hat die Gl. (5.1) die Lésung (5.3). Hierbet 
ist hE L? mit fe L. 


Als Beispiel hierzu betrachten wir 
t 
fel f(u)du+2t—1)=f(), t>0. 
0 
In der Operatoren-Schreibweise gilt dann 


fe} {FO} + 2{(¢ —1) } = {fF O}. 
Hier ist 
{?—-1 


k= {ol}, i= je ot, fiir m= 2,3,0.. 


und 
baler rhs ty}. 
Nach Satz 5.1 folgt, 
f(t) =2(¢—1) & + fern 2(u —1) edu 
= 2(t — 1) e* + e?!(e"—?#— 1) 
= (2t—1) e*— 


und in der Tat erfiillt diese Losung die gegebene Gleichung. 


Es ist zu bemerken, daB in diesem Beispiel die Anwendung der Laplace- 
Transformation nicht unmittelbar™ modglich ist, da das Integral 


fe~S*(t—1) edt 
i) 
divergiert. 
Beispiele die SCHMEIDLER ([25] S.250) und TitcHmarsu ([23] S.313) mittels 
der Laplace- bzw. Fourier-Transformationen lést sind auch vermége der obigen 


Methode zu losen. 


1 Oder man miiBte bei der Behandlung dieser Integralgleichung das_,,Fort- 
setzungsprinzip“ (s. [3] Bd. II, 5S. 259) benutzen. 
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6. Systeme von Integralgleichungen 


Wir betrachten ein System Volterrascher Integralgleichungen zweiter Art 


AW —D fhylt—w) fw) au = a(t) 
(6.1) MN Reta Bowe 
if (t) = > fi ky v (t a. u) is (w) du = En (t) 


y=1 0 
wo die Funktionen {7, (é)}, {/o(é)}, ..- {f, (t)} unbekannt sind. Dieses System 1aBt 
sich mittels dieser Theorie in eleganter Weise auflésen, indem man die Matrizen- 
rechnung benutzt. Setzen wir 


f(t) x(t kyi(t) Ryo(t) - k, ,(t) 

eG eee a 

% (1) = fa(0 6 () = 82 (2) R() = hoi (0) Roald) an(f) 
i, (t gn) Fnrld) Bao(t) -»- ban (® 


wo §(t) die Kernmatrix des Systems (6.1) ist, so laBt sich das System in der 
Gestalt 


t 

(6.2) $(t) —[ RE — u) F(u) du = Gt) 

0 
schreiben. Die Matrizen sind Funktionen von ¢ und das Integral einer Matrix 
ist diejenige Matrix, deren Elemente die Integrale der entsprechenden Elemente 
der gegebenen Matrix sind. Auch ist ® (t — u) & (wu) im Sinne des Matrizenproduktes 
von ®(¢—u) mit (uw) zu verstehen. (Auf die Nichtkommutativitat dieses Pro- 
duktes ist zu achten.) 

Nun lést man die Gl. (6.2) in ahnlicher Weise wie eine einzige Volterra- 
Integralgleichung zweiter Art auf. Ist F$={§ (‘)} die Spaltenmatrix, deren 
Elemente nun /; = {f;(¢)},7=4, 2, ..., n sind, und gilt entsprechend fiir 6 = {6 (#)} 
und Rae (i)$, so nimmt das System in der Operatorenform die Gestalt 


(6.3) F—KF=G 


an, wobei 


ist. Besteht die quadratische Matrix 8 bis auf die Elemente 1 (der Einheits- 
operator in K) in der Hauptdiagonale, aus lauter Nullen {0}, so ist die Lésung 
von (6.3) formal durch 


5 =(3— K)+G 


gegeben, falls die Matrix (8—S) nicht-singular ist (d.h. die Determinante 
| 3 — ®| + {0} ist), und (S — ®)+ die Kehrmatrix bedeutet, die durch die gewohn- 
lichen Methoden (z.B. Kowatrwskt [8] S. 130) erhalten werden kann. An- 
schlieBend findet man die Lésung von (6.2) im Originalraum, falls sie existiert. 
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7. Integro-Differentialgleichungen und nicht-lineare Gleichungen 
Technische” Aufgaben fiihren zuweilen auf Gleichungen, in denen die ge- 
suchte Funktion auBer unter einem Integral auch noch in ihren Differential- 
quotienten vorkommt. Die allgemeine Gleichung [4] dieser Art (vom Volterra- 
schen Typ) 


TT Ona Sa HO) $+ Hooft) + LAE 4) flu) du =e 


ah = = 
TG 


1a8t sich mit den vorstehenden Mitteln lésen. Anstatt eine allgemeine Theorie 
zu entwickeln, illustrieren wir die Methode durch ein Beispiel, in dem wir fol- 
gende Gleichungen lésen: 


f® +1 + few f(u)du=a  (£>0). 
Ist {f(#)} total-stetig und /(0) =0, so gilt 
sf+f+——f=al 


oder f =a(s —1)/s3 woraus f(t) =« ( = 5} folgt, welches tatsachlich die Gleichung 


erfiillt. fs 

Die Operatorenmethode ist nicht nur auf lineare Integralgleichungen be- 
schrankt, sondern auch anwendbar auf Gleichungen, in denen die unbekannte 
Funktion f mit sich selbst und mit gegebenen Funktionen mehrfach gefaltet 
vorkommt. Wir betrachten folgende quadratische Integralgleichungen (s. auch 
Port [14]): : 
2f(t) —f f(t —u) f(w) du =sint. 

0 


Die Operatorenform dieser Gleichung 


1 
2 | foxy 
/ 2i T s2t1 OF 


eine algebraische Gleichung zweiten Grades, hat die zwei Lésungen!?* ([7] S. 159) 


fatto, =1 4 s{hi0}. 


Da s {Jo(t)} ={Jo ($+ Jo(0) und Jo(t) =—S (4), Jo(0) =1 ist, sind die Lésungen 
f=1 Sis (— {/,()} +1). 


Die eine Lésung /={/J,(é)} ergibt die gewéhnliche Funktion /(t) =/,() wahrend 
die andere, namlich f=2—{ J, ()}, nicht auf eine gewohnliche Funktion™ redu- 


zierbar ist. 
Weitere Beispiele, die auch mittels dieser Methode lésbar sind, hat L. Port[14] 
durch Laplace-Transformationen geldst. 


12 Zum Beispiel Leon LICHTENSTEIN: Vorlesungen tiber einige Klassen nichtlinearer 
Integralgleichungen und Integrodifferentialgleichungen. Berlin 1931. 

12a Nicht jede algebraische Gleichung #-ten Grades ist im Operatorenkorper K 
losbar (z.B. Beispiel in [17], S. 159). 

13 Hier entspricht der Operator 2, der Diracschen Delta-Funktion (mit der Zahl 2 
multipliziert). 
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8. Singulare Integralgleichungen 


Hat der Kern einer Integralgleichung die Form k (#) =7(#)/, «>0, wo z.B. 
j (t) stetig ist, so wird die entsprechende Gleichung singulay genannt. Nehmen wir 
zuerst den Fall 0<«<1, so hat das entsprechende Integral trotz des Unendlich- 
werdens des Integranden einen Sinn. Zuerst betrachten wir als Beispiel die 
Abelsche singulare Gleichung erster Art'* 
t 


i io du = g(t) (0 <0 =< 1) 


¢— uy 


0 
fiir g eine total-stetige Funktion fiir 20. Als Lésung ist 


fy) = 2224 | fe — wy ew) du + 6(0) I]. 


0 


In der Tat, hat die algebraische Gleichung 


SBE on (0< a<-1) 


die Lésung 
— 1 (od 
Nea 
1 aye 
i Rareer te [{g’()} +8 (0)]- 
Da [(4—«) I’(«)]1= (sin «)/z ist, folgt die Behauptung. 
Die Abelsche Gleichung zweiter Art 


t 
Ap fuk Ces. hx a8 Ne 
va | ia ae eect 


ist in der Operatorenform 


Bet haf 


darstellbar. Da die Funktion /? die Voraussetzungen von Satz 3.3. erfiillt 
nimmt die Lésung die Form 


Feet | eae tay 
{ 4 


=p y yet 


n=1 


‘eta Pn ae 
an. Da weiterhin ae fiir Re A>0 ist ([ZJ] S. 101), gilt 
i(t) ae aig ! 1 7 {2/n 
Vat ~ Vi r(*+*) ; 
n=1 2 


Fis NE See ha canes age an, daB f in jedem endlichen Intervall 
1S? St, beschrankt ist. Dann existiert das Fal i ii i 
Dea ae altungsintegral fiir alle ¢ und ist 
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Nun gilt bekanntlich (z.B. [26] S. 507) 


co 


Pleaa al 4 erf yi) 


n=1 2 


4 
wo erf t= (2/|/z) f e~“ du, das Fehlerintegral bedeutet. Die Lésung der Integral- 
0 


gleichung nimmt also die Gestalt 
1 
t) =——+é(1+ erf jt 
MO) = 7+ (1 + erf 2) 


an. Diese Lésung erfiillt die gegebene Gleichung. 


Bisher war O<a<1. Hat aber der Kern die Form 7 (é)/¢*, mit « =1, so miiBte 
man das Integral der Gleichung, damit es iiberhaupt einen Sinn hat (4/¢ ist nicht 
integrierbar), im Sinne des Cauchyschen Hauptwertes z.B. verstehen. Fiir den 
Fall «>1, kénnte man z.B. das Integral im Sinne des Hadamardschen ,,partie 
finie’’ interpretieren (z.B. ScHwarz [17] Bd. I, S.38). Der Verf. hofft in einer 
spdteren Arbeit singulare Integralgleichungen dieser Art in Betracht zu ziehen. 
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Asymptotic Theory of Second Order Differential 
Equations with Two Simple Turning Points 


NIcHOLAS D. KAZARINOFF 


Communicated by L. CESARI 


1. Introduction!. In 1931 Langer initiated and gave the first of numerous 
contributions to what has become a successful theory for asymptotic expansions 
of the solutions of a differential equation with a turning point. This theory has 
been extended and applied to a great many questions by him and by others. 
An extensive list of references may be found in L. Cesart’s book, Asymptotic 
behavior and stability of solutions of differential equations, Chapter IV, Springer- 
Verlag, Berlin, 1958. 

In this paper we present a general asymptotic theory of ordinary second order 
linear differential equations with two simple turning points and containing a 
numerically large parameter. In particular, we are concerned with the asymptotic 
expansions with respect to complex / of solutions of differential equations of 
the form 


(1.1) YG? P(s, 1) y= 0, 


We consider this equation for s in a closed, simply connected, perhaps unbounded 
region Y, of the complex plane. We assume that for |A|>N and for sCQ, the 
coefficient P(s, A) 1s of the form? 


(1.2) P(s, A) = DP, (s) A, 


where each p,; is analytic, and, most importantly, that py(s) has precisely two 
simple zeros, « and B, in the interior of Z,. 

In a region which includes a turning point, 7.e., a zero or singularity of f(s), 
the solutions of the differential equation (1.1) depend upon 4 in so intricate a 
way as to have quite distinct asymptotic forms in different parts of the region, 
being dominant (exponentially large) in some parts and subdominant (exponen- 
tially small) in others. The asymptotic series for solutions of (1.1) over a region 


1 This research was supported by the United States Air Force through the Air 
Force Office of Scientific Research of the Air Research and Development Command 
under contracts No. AF 18(600)-1481 and No. AF 49(638)-192. 

2 The letters MM and N are always to be used as generic symbols for positive 
constants. 
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which includes just one simple turning point are known and are based upon 
Airy functions [2,4]. Thus, it would appear that the behavior of solutions of 
(1.4) over the entire region Y,, which contains two simple turning points, 1s 
obtainable by the familiar procedure of evaluating the coefficients in the depend- 
ence relations which connect solutions whose behavior is known about one of 
the turning points with those solutions whose behavior is known about the other 
turning point. Unfortunately, this evaluation is not possible in general; two 
linearly independent solutions may have the same dominant asymptotic form, 
which makes the inference of the identity of two solutions from the identity of 
their asymptotic forms invalid. Thus, a new theory is necessary if we desire to: 
have uniform asymptotic expansions of solutions of (1.1) over all Z,. LANGER [4] 
has derived the leading terms of such asymptotic expansions in the special case 
where s is a real variable on a bounded interval, f(s) is real-valued, and #,(s) =0. 
We derive the asymptotic expansions to +1 terms, where m is an arbitrary 
non-negative integer, of the solutions of (1.1) under the general hypotheses set 
forth in the first paragraph above together with some others of a more technical 
nature to be set forth later. 

Interest in the problem discussed here stems mainly from possible applications 
for the theory derived. In certain regions, the differential equations for the 
angular and radial spheroidal functions are of the type (1.1). This is also true 
of the Whittaker equation for certain configurations of its parameters. The 
spheroidal functions are important in problems concerning scattering by a prolate 
spheroid. The Whittaker functions, disguised as Coulomb wave functions, occur 
in quantum mechanics. Equations of type (4.1) are also of interest in other 
problems of wave motion and diffraction. 

We have divided the discussion below into three parts. In Part I we transform 
the differential equation (1.1) into one more suitable for analysis. We call this 
canonical form the given equation. We then give an algorithm for the construction 
of a related equation whose coefficients resemble the coefficients of the given 
equation to an arbitrarily prescribed degree. In Part II we study the solutions 
of the related equation. These involve Weber functions of large complex order 
and argument. They have been studied by ERDELYI, KENNEDY & McGREGor [8]. 
We make considerable use of their results and have shown that their asymptotic 
representations hold uniformly in arg y over a finite range of arg, where y is 
the order of the Weber functions involved. We also give an algorithm for re- 
cursively determining the terms in the asymptotic expansions of these Weber 
functions. In Part III we prove that solutions of the related equation are asympto- 
tic expansions to 7 +1 terms of solutions of the given equation, being any 
non-negative integer. The method of proof is, as usual, to transform the given 
equation into an integral equation of Volterra type, whose kernel involves solu- 
tions of the related equation, and to solve this integral equation by the Picard 
method of successive iteration. We also give approximations for derivatives of 
solutions of the given equation. Our main results are stated as Theorems 1 and 2: 
of §§ 12 and 13. An especially interesting feature is the occurence of a denumerable 
number of characteristic solutions of the given equation. These are bounded and 
oscillatory on certain curves joining the turning points and exponentially small 
on their extensions beyond the turning points. 
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Part I: Construction of the Related Equation 


2. The given equation. The analysis to follow involves certain functions 
and a mapping connected with #,(s). To simplify their form we adopt a normalized 
form of the differential equation (1.1) as the basis of the sequel. This normaliza- 
tion was used by LANGER in [4, §§ 2, 3]. Its adoption here will be seen to entail 
certain assumptions on 9(s). 

We consider the mapping from the given region Y onto a region Y of the 
z-plane defined implicitly by the equation 


(2.1) J (®@—1)8dt— f pe) dt =0. 
= 8 od 

An immediate question is whether or not (2.1) defines a mapping at all, and if 
so, is it a schlicht mapping? We can show that (2.1) defines a mapping and that 
at each point of ZY it defines a locally one-to-one mapping. The function which 
is the left member of (2.1) is analytic in z except at +1 and ins except at « and f. 
Its partial derivatives do not vanish except on the lines z= +1 and s=z or f. 
Thus by the implicit function theorem for analytic functions, there exist in 
neighborhoods of all points of Q, except possibly « and #, analytic solutions 
2(s) of (2.1) with inverses s(z). Further, ds/dz=(z®—1)!p)#(s); and hence, 
s’(z)==0 except perhaps at 4: 

We now examine what happens at the exceptional points « and f. A com- 
putation shows that we may write 


(2.2) f(@—1)$dt=(2—1)8@,(z) and f ph() dt =(s — a) B(s), 
1 a 

where ®, and P, are analytic in neighborhoods of +1 and «, respectively, and 

neither ®,(1) nor P,(«) is zero. Thus, we may write 


(2.3) F(z, s) = (2 —1) OF(2) — (s — a) Bis) =0, 


where F is analytic in a neighborhood of (1,«) and neither 0f/0z nor dF/ds 
vanishes at (1,«). Therefore, by the implicit function theorem, there exists a 
solution z(s) of (2.3), and hence of (2.1), with inverse s(z), which is analytic in 
a neighborhood of « and such that z(a) =1 and s’(1) +0. 


In order to apply this same technique to (2.1) at s =f, we must assume a 
=i 


B 
normalization of the coefficient ~). The integrals f (—1)#dt and f p3(t) dt 
1 


a 
are independent of path provided the paths do not encircle the singularities of 
the integrands. However, it may be that the latter integral vanishes. We assume 
the contrary and choose the parameter so that? 
B =i 
(i) Spe) dt =f (@—1)8dt =— pat. 
e4 a 

We can now apply considerations analogous to those above to reach the con- 
clusions that (2.1) has an analytic solution z(s) at (—1,f) with inverse s(z) 

3 We do not consider the limiting behavior of the solutions of (1.1) as the turning 


points a and f# approach each other. F. W. J. OLvER has informed the author by a 
private communication that he is examining this limiting behavior. 
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whose derivative s’ is nonzero at —1. If arg #9 is constant on the line segment 
joining « and f, the hypothesis (i) 1s clearly fulfilled. This is the case treated 
by LANGER [4]; and since, in this instance, the mapping defined by (2.4) is one 
to one on the line segment joining « to f, it is schlicht in an open neighborhood 
of this closed segment. 

We have shown that the implicit relation (2.1) defines a mapping from J 
onto a region J and that it defines a locally one to one mapping of some neighbor- 
hood of any point of Z onto a neighborhood of ZY. We require more than 


this and assume that 
(ii) The relation (2.1) defines a schlicht mapping 
2(s): RoW, 
with [—1,1] CZ. 
Under this assumption, the change of variables 


(2.42) s=s(2), 


(2.40) y(s) = {5b ula), 


es Hu _ 720(2,2) u=0, 
with 
OG =P) As*+ilger— az} (ae) 


where Q is analytic for |A|>N and for z¢Z. We write 


Q(z, 4) = 2 q; (2) Av?. 
In particular, 
go(2) = (2 —1). 


We henceforward refer to the differential equation (2.5) as “‘the given equation’. 


3. The first approximating equation. The analysis of the given equation 
is based upon the construction of an equation which resembles it up to terms of 
the form A~"~* O(4) in the coefficient of wu, where n is any non-negative integer 
and O (1) denotes a function of z and 4 which is bounded for z€ Z and for |A|>N. 
The algorithm for the construction of this related equation is similar to the algo- 
rithm of McKELVEy given in [6]. The matter at issue in MCKELVEY’s paper 
is the approximation of solutions of an equation with a single turning point of 
order two. Such a turning point may be thought of as the confluence of two 


simple turning points so that similarities in the analyses of the two cases are 
not unexpected. 


The point of beginning in both instances is WEBER’S equation, which may 
be written in the form 
av 


v t2 
ae 1 


2anrG 


|V =o. 
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We make the transformation 
t= (22)87, “Vir =v), 
and WEBER’S equation becomes 


“1 Ade —9) 0 =0. 


Let a non-negative integer 2 be chosen once for all. We now choose » so that 
(3.4) HS) Ore) i as 
0 


where the (arbitrary) constants c and c,; are yet to be determined. This brings 
WEBER’S equation into the form 


(3.2) ooo A? (22 —1) + Ac cpg) v=0 or v'—Rliz,A)v=0, 
0 


which is our first approximating equation. 


4. The second approximating equation. We next construct an equation 
resembling the given equation in both the 4? and 4 terms of the coefficient of w. 
Formally, this step is almost identical with the corresponding one in [6]. The 
idea is to make a change of dependent variable 


(4.1) LE gt + iva 


in the first approximating equation (3.2) and to determine the functions j1, 
and yw, in such a way as to make the new differential equation in Z more closely 
resemble the given equation. By differentiating Z twice, employing (3.2) at each 
step to replace v’’ by Rv, and by constructing the eliminant for v and v’ among 
the three relations connecting Z, Z’, and Z’’ with v and v’, we find that the dif- 
ferential equation satisfied by functions Z of the form (4.1) is 


Woe und At ye Mies 
ony ZF itt Dey 22 tic ai 
wherein 
Mo Ho + py RA | 
4.3a DZ A) = ‘ 
(4.3) 0 (2,4) Sve Renan 
bo fo + 2m, RA2+ py, R'/A | 
4 b H(z,A —— , ’ , 
ae?) ae) y/A 2Mo + ya/A 
and 
(4.30) Tea) — [M0 THOR + 2m RAM + RIA ahem 
. Quo + or At + py RIA fo + wild 


In so far as A is concerned, the functions Dy and H are bounded for |A|>N 
with Dg (z, 0) ="0— 141 9o- 

We shall determine jup and jz, by the condition that the terms in the coefficient 
J [Do in (4.2) which become infinite with A have sum A?q9(z) +49(z). We recall 
that the choice of the constant c in (3.2) is in our hands. Now, a simple division 
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shows that 


(4.4) 2 =R+Dp" 


0 


Ho + 24 RAA+ wy RA potty a 

2p + pr /A Mo + palA 
Recalling the definition of R from (3.2), we observe that the coefficient of A? in 
J|Do is q- The coefficient of A may be computed and is 


211 Jo + M4 qo 14 %0 


Cc Ds Z, 00 ; 
+ Do (, e) au * 


’ 


or 
Do" (z, 00) [(c Mo + 2Lo [4 Jo + Mola qo) meee Ma Jo + 204 qo) | , 


Therefore, a sufficient condition that this coefficient be equal to q, is that srg 
and py, satisfy the following system of differential equations: 
211 Jo + 1190 = Mo lM — ©), 


210 = ta (%1— ¢); 
with boundary condition that 


10 — 144 Go = Do(z, 20) = 1. 


A possible solution of this system is 


(4.5a) Me = COS, a= = e ; 

where 

(4.5 b) Q(z) == (22 4)\), = are p= OP for oe 
and 

(4.5) B(2) = Data oo ge 


27 (s) 


We deduce from these formulas that wy and m4, are analytic in Y except at 
%=—1 provided yw, is defined so as to be continuous at z=1. Without proper 
choice of the constant c, #y and yw, will fail to be analytic at —1. In order to 
render them analytic there, we choose c so that 


(4.6) C= eh 


Then #(—1)=0, and # may be written in the form 


27(s) 


0 (z) — n(s) ae, ds 


Provided 1, is now defined so as to be continuous at — 1, this fact and computa- 
tion show that 4) and mu, are analytic throughout Z. 
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The choice of fv) and jw, determine Dy, and from (4.5) it follows that 
(4.7) Ho — Ma Go = 1- 
The relations (4.7) and (4.3) and the definition (3.2) of R imply that 


(4.8) Dy (2, A) = 1 = AA (tg Ya — Ho Hu =e mt) — DG A *. 
0 
Thus if Y is bounded, it is clear from (4.8) that Dy is bounded away from zero 
for |A|>N and for z€Y. If Z is unbounded, we assume 
(iii) Do ts bounded away from zero for z€ QB, and for |A|>N. 


We observe directly from (4.3) that H=D . Therefore, we may make the 
change of variable 


(4.9) ¢=D,iZ 
to remove the first derivative term in equation (4.2). It then takes the form 
(4.10) © — [22 =1) SAG ee Ac = 0. 
where 
2 (52 aeons 
(4.10b) BGA) hg Pst ee on 


The function T is analytic in Yas R, wo, and py are all analytic there. Further, 
T is analytic in A for || >WN; hence, we may write 


(4.11) F(a, a) = Sa (—. 


We call the differential equation (4.10a) “the second approximating equation”’ 
for the given equation (2.5). It is important to note the way the constants c,, 
which were introduced in (3.1), enter into the functions ¢;. The c; appear in T 
wherever R and R’ do. From (4.4) and (4.3a), we observe that aside from the 
leading term R in the right member of (4.10b), R and Kk’ have coefficients of 
order A+ or smaller in this right member. It follows that the constant c; occurs 
linearly in ¢;, always with a coefficient of +1, and is absent from every ¢, with 
Te 

5. The related equation. We are now in a position to construct a dif- 
ferential equation whose solutions are known and whose coefficients are identical 
with those of the given equation up to terms of order 4~”~'. Inasmuch as this 
construction is already well known [4,6], we present only the results. 


With ¢ standing for any solution of equation (4.10a), we set 


1 Biel 
(5.1a) w =Dz#(AC + iT), 
where 
AZ A=), 4,2) 72, 
(5.1 b) : 
B(z,a) =D 8;(2) 47, 
0 
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and a 
Ae 

(5.1c) D, (z, 4) = 5, lh tea” 
ru ae 


and choose the functions a; and b,; so that w satisfies a differential equation of 
the form 
Q(z, 4) 


(5.2) oe | 0, =e |@ =O. 


where Q(z, 4) is bounded for z€Q and for |4|>N. The proper choices of the 
a;s and },’s are: 


; 2 —' 
A a ae te : | a 


, d 
a,(z) =0, — 29, by — bo) sory 
(5.3) and for7=2,..., 
Cig = 
a; (z) = sf . +> (Jj—e = t;_p—2) by| ds, 
: 0 
1 : j , , 
b,@) = —— i bs (Qj—-r+e— be) Oe — (2% 8-1 + 91 51) — 
eI 14 
hee ds 


ard (Dls paca Opt yp Oa ian a 20° 

0 
In order to render the a,’s and 6,’s analytic in Y, the constants c; occurring in 
their defining expressions must be properly chosen. This may be done recursively 
since c; occurs linearly in ¢; and is absent from each ¢; with ;<7. The correct 
determinations are 


1 
aoe = aie Seen a eat ds 
J 20 Cy yJ CE 1 150 — 41 Bo] 20” 
(5.4) and for 7=2 


1 


Parad 
=y | bx Qpthie = BS 5) a, — (29 0;-1 + 91 5,1) t (Gj+2 t;*) 


-1 


7-2 
—a 2 (2t;_4-2 Dp A Gites by) = a; ss 


where 


We now observe that the division by D? in (5.1) is legitimate if QZ, is bounded, 
since D,(z,0)=1. If ZY is unbounded, we assume 
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(iv) D(z, a) 1s bounded away from zero for z€ B and for |A|>N. 
The function £2(z, 2) which appears in equation (5.2) is a function whose specific 
form may be computed, although it is of no interest to us. It involves the several 
functions 4, B, Q, R, gy, and q,: It is, however, important to note that Q(z, A) 
is analytic in A for |A|>WN and is analytic in z for z€Y. We henceforward 
refer to the differential equation (5.2) as the related equation. In Part II we 
shall single out certain of its solutions and describe their behavior for z€Z 
and |A|>N. 
Part II: Solutions of the Related Equation 


6. Introduction. The construction just given of the related equation enables 
one to know the behavior of its solutions only as well as one knows the behavior 
of the Weber functions which are solutions of the first approximating equation 
(3.2). These are Weber functions of large complex order and: unrestricted argu- 
ment, ¢.g., Ds(,-3)(/24z). Their behavior is not yet completely known. ERDELYI, 
KENNEDY & McGREGOR have derived the results most useful to us [3]. They 
have established approximations to the various functions 


(6.1) Vin (Xs ”) = D(ov—1)/2 (V2» eae 7) , 0) == (— 1%, 


with an error involving y+. We strengthen their results by showing that their 
asymptotic representations hold uniformly in argy over sectors of the v-plane. 
We also extend their representations to asymptotic expansions, the terms of 
which are determined recursively by a quadrature at each step. 


We obtain the behavior of solutions of equation (3.2) by letting xa|/4e 


in the expansions for the functions (6.1). We then determine the asymptotic 
expansions over Y of the solutions of the related equation. These are explicit 
up to terms involving y+. While our results are limited by the lack of precise 
information on the structure of Weber functions of large complex order, ap- 
proximations to terms involving y+ are usually adequate in applications. When 
they are not, numerical methods for the calculation of the terms of the asymptotic 
expansions given in § 8 may be employed to give more precise approximations. 


7. Domains in the x-plane. The determination of the regions of validity for 
the approximations to solutions of the given equation, which are derived in 
Part III, depends upon the character of a function @(x). It is defined by analytic 
continuation from its positive values on (1, oo), which are given by 


(7.1) D(x) =e 1)2dt =F ol dt. 


It is shown in [3, pp. 469, 470] that ® is a schlicht mapping of a Riemann sur- 
face X over the x-plane with branch points at +1 onto a Riemann surface 8 
over the ®-plane whose branch points are the images of the branch points in &. 


Clearly, 
(7.2) P(x) = 3 {x (x? — 1)! — In[x + (x? —1)4]}, 
and for large |x|, 
(73) P(x) = + }2[1 40(29)]. 


10* 
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The Weber functions (6.1) are entire functions of their argument and their 
order. Hence, in discussing them we may assume that 


(7.4) The x-plane is cut from —1 to +1 andlargy|Sa+e, e>0. 


These assumptions enable us to avoid considering more than a finite number of 
sheets of the logarithmic Riemann surface 3. 


Let X, denote the x-plane cut in accordance with (7.4) and the condition 
|arg(x—1)|<a. The corresponding portion 3, of 3 consists of two sheets joined 
along the ray arg®=0 (see Fig. 1). We now 

ag¢=--2 define certain regions of X on which the be- 
havior of the Weber functions (6.1) is con- 

veniently describable. Let the index m, which 

appears in (6.1), have range 0, +1, +2, +3, 


a 2L<-7,0°9 L=— 


o= arg @ =0 and let 
(7.5) E(x, ») =v O(a). 
Or arg B=0 
. Let y,, be the point at infinity on 3 in the 
direction arg ®=—ma—argy, and let ¢ be 


any positive number small enough to make 
the following definition meaningful. We define 
&,, to be the largest closed, simply connected 
Fig. 1. The surface 3, region ® of X such that 


= — arg D=2N 
Lé-1,019XL =H 


(7.6a) —#?a+eSargét+mnaSin—e, xe@, 
(7.6b) |x +1/Se, EG, 


(7.6c) |P(x)|=2e, x€@, except that x=1 and a neighborhood of x=1 may 
be in &. 


(7.6d) Each point x in & may be joined to the image of y,, by an analytic, 
simple curve lying in ®& and on which R& is monotone (monotone in- 
creasing from the image of y,, to x if m is odd, monotone decreasing if m 
is even). 


The regions ©,, are illustrated in [3, p. 482]. In general, except for sectors with 
vertices at x =-++1 and widths 2¢ in arg(« +1) and the interior of the circle 


at *——41 deleted by hypothesis (7.6b), each region ,, covers the x-plane. 
Exception occurs, for example, when m=O and arg» =0. 


We denote by G;, the reflection of &,, in the origin; that is, 


(7.7) Gr={ale"'xeG,, if argxZo, ec xCG,, if arg x So}. 
Now 
(7.8) Der" x) = P(x) Sie ini, 


so that ®F.20G,, is precisely &,, plus all but a small sector with vertex at 
*=—1 and width 2¢ in arg(x%—1) of the neighborhood of x= —4 not in G,,- 
The definitions of &,, and ®% are essentially the same as those given in [3]. 
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8. Asymptotic expansions of Weber functions. In this section we describe 
the asymptotic behavior of the Weber functions (6.1) over the regions &,, and 
(7. These functions are solutions of the transformed Weber equation 


7 ad? 
(8.4) 7 v2 (x2—1) y =0. 


Since they are entire functions of x, y,,=y, if m=1 (mod 4). The Wronskians 
of certain pairs of these solutions are [/, p. 42]4 


Fe alg 
(8.2a) W (Vins Vice) eseeee 5 
2 
(8.2b) Wms scr» 6) = OTRO YZ 5 ebm, 
where 
dg d 
(8.2¢) We fe ge 
and 
(8.2d) Go =(—1)”. 


We henceforth adopt the convention that in formulas in which the double signs 
-++ or -- appear, all upper signs or all lower signs are always to be used. It follows 
from (8.2) that any two of the four distinct solutions y,, are linearly independent 
except for y =O and certain other integral values of v. If y is not zero, and we 
always consider |»| to be large, the solutions y,, and y,,., are linearly independent 
without exception. 

The asymptotic representations for the functions y,,(x%,v) involve the Airy 
function Ai(¢), which is defined by the relation 


(8.3) Ai(t) = t Re (= i) 


in which K, is a modified Bessel function of the third kind [3, p. 463]. For each 
admissible m, the function 


(8.4) X q(x) = 9 (a) 


2 nt QE 2 4 
6) x i Ai]2( OSs yom, i) = as 
42 (*) 42 


has a simple asymptotic expansion for x€,,. This is 


wa 8r ml 


(8:5) X,,(%, = 2 ta ty bey 8 (x) e-9F- apie: i ae aco ae Oa 


AO 0 — On Se 2, a 


The asymptotic series for dX,,/dx for «€®,, may be obtained from (8.5) by 
differentiation. ErbDELy1, KENNEDY & MCGREGOR have shown that the func- 
tions (8.4) are asymptotic representations for the Weber functions (6.1) [3, p. 479]. 
These representations may be extended to asymptotic expansions by using the 
algorithm given by LANGER in [4] or similar algorithms. We use LANGER’s. 
In the notation of [5], it enables us to construct certain functions «; and f; 


4 The formulas in [/] appear to contain misprints. 
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which are involved in the coefficients of the asymptotic expansions for the Weber 
functions (6.1). With reference to equation (8.1) the defining expressions of a; 
and f; specialize. They become: 


ola) =1, fos) = f AO at, 


ie 29 (2) 
ae a(x) =0, B(x) =0, 
and for 7= 2, 
aj(x) =S* f [Bi-al) +20 B,2(0) a2, 
(8.6b) ore ’ 
Bi) =a f {ai +20 los) + 26;-201 —¥ OB, 20} a0» 
wherein 
(8.6) k(x) = and W(x) = G4 (x) p-8(2). 


W (x) 


Since %,=f,=0, O%;+1=f2;41=0 for all 720. The following order relations are 
easily established. They are important in the extension of the asymptotic 
representations in [3] to asymptotic expansions. For || >N and 720, 


(8.7) O%;(4) =O(1), and  f2;(x) =O(x™). 


The asymptotic expansions given below extend the results in [3]. We omit 
their derivation to avoid repetition of work in [3]. The «,’s and £,’s are suf- 
ficiently small at x = oo to guarantee uniform convergence of the integrals which 
enter into the derivation. This is the point in the extension of the analysis in 
[3] that requires careful attention. The bounds for the error terms in our 
expansions are uniform with respect to arg v, if arg y lies on a bounded interval. 
This also applies to the results in [3]. The reason for this is that the limit a,, 
of integration in [3, equation (4.8)] may be fixed provided argy has a range of 
length less than a. For example, if |argy|<$a—e, ay may be chosen as 
x =e". 00; and if |argy—in|<ia—e, ay may be chosen as x =e7!%*- oo, 
For each such range of argv, a bound on the integral in [3, equation (4.8)] may 
be found which is independent of argv. It follows that for arg »y bounded, bounds 
for the error terms independent of arg» can be found. This reasoning does not 
imply, however, that in the x-plane the regions of validity for the asymptotic 
expansions below are independent of argy. They are not, although in the 
neighborhood of infinity the boundaries of the regions of validity in the |vx- 
plane are asymptotically independent of arg ». 


Na, PeeaIC), Da (x. hho xe 


(8.8a) c 
x {A¥ (a, v) X,,(x,») [1 +O(-7-)] +2 ee ”) ule) 


ax 
for x €G,,; and 


Ym (%,¥) = Cy, [Dy (x, »)]~* x 
x {AM DX 2) = eee ar FOO, 


y dx 


(8.8b) 
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for all x €{x| |x —1|<e} and |E(x, v)|<N. In these relations 


(8.9a) Ep. — Qt ypilr—d ¢ 4 [ov (1—mai)+ aL 
Di (x, ») 
(8.9b) 
_ |A*(x, »)- SSN) 
A*"(x,v)—(x?—1) B*(x,v) +k (x) B¥(x,v) v2 A*(x,v) +"? B¥'(x, 2) : 
and 
mas : qa : 
(8.9c) AeY) == OA) | WET (X,?) => B (ao «. 
0 0 


The functions «;,8;, and k are defined by the relations (8.6), the functions 
X,,(%, v) are described by the formulas (8.4) and (8.5), and o=(—1)”. 


One uses the identity 

(8.10) Vmn(%,¥) = Va (er** x, 9) 

to derive the form of y,,(x,v) when x€@%,.. Taken together the regions G,, 
and ®F, (m=0, +1, +2, +3) cover the x-plane, each point being in two regions 
having indices which differ by 1 or 3. Therefore, the expansions (8.8) yield the 
asymptotic behavior of two linearly independent solutions of the differential 
equation (8.1) for each x. We observe that zeros of the Airy function are not 
excluded from the region where the relation (8.8b) is fulfilled. This is the reason 
for the altered character of the error term in (8.8b) as compared with the error 
term in [3, equation (9.6) |, for example. The behavior of dy,,/dx may be found 
via the relation 


(8.11) 


from the expansions (8.8). 


dD, (« 1 
ee 

It is interesting to note that if » is real, all solutions of equation (8.1) are 
oscillatory on the interval —1<*<1. This occurs since, for the configuration 
of ~ and y cited, #(é(x, v)) =0. 

9. Solutions of the related equation. The first approximating equation (3.2) 
has solutions 


(9.1) Um (2, A) =D aman err. 2) (m =0, +1, +2, pet is 


where D,(z) is the standard parabolic cylinder function of order y [7], and where 


(9.2) PeaA=e CA 
0 
Of course, v,,=v, if =m (mod 4). The formula 
[2 
(9.3) Uy (5 A) = Im( | 2) 


relates the solutions of equations (3.2) and (8.1). In view of this connection, 
the asymptotic behavior of the solutions v,,(z, 4) may now be easily described 
with reference to the x-plane. Since 


(9.4) Ly \4 2, 


the transition from regions on the x-plane to their images on the z-plane is 
elementary. 
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By directly applying the theory in [4] to equation (3.2), we could have 
described the behavior of certain solutions of (3.2), other than the solutions (9.1), 
on regions in the z-plane identical with the regions &,, and @*. However, the 
slight simplification in the description of the regions of validity of the asymptotic 
expansions, which is gained by this approach, is offset by an increase in dif- 
ficulty in deriving these expansions, since equation (3.2) is superficially more 
complicated than equation (8.1). For this reason and to take fullest advantage 
of the work in [3], we proceeded as we have done. 

We compute the Wronskians of the solutions v,,(z, 4) using the formulas (8.2), 
(9.3), and (9.4). For future reference, the computed values are: 


2VaAe8™** 
Wm Uma 2) = —ov-+1)\ 
Pe 


(9.5) 


LA aki? spr etini(ltor) \22 ghmni 


To each solution v,,(z, A) of equation (3.2) there corresponds, through the 
relations (4.1), (4.9), and (5.1), a solution w,,(z, A) of the related equation. The 
solutions w,,(z, A) may be written in terms of the functions v,,(z, A) as follows: 


(9.6) Wy (2, A) = Eq(Z, A) Um (2, A) + Ey (2, A) Un (2, AA, 
. where 
3 B R “ D 
ue Egle, A) = (DyD3)4|A iy +5 (PR 4 oe — oe 
: ar B f Ds 
Bx(@,4) = (DoD:)-* [An + (to + Be 
Since 
(9.8) Ww, Wms 2) aay Wr, On 2) ? 


the linear independence of solutions w; and w,, depends upon that of their cor- 
respondents v, and v,,. 

For the proofs to follow it is necessary to have at hand the asymptotic 
behavior, in first approximation, as |A|-> oo, of the solutions w,,(z, 4) of the 
related equation together with the behavior of their first derivatives. One finds 
this behavior by considering the relations (9.6) in conjunction with the asymptotic 
expansions (8.8) and (8.5) and the relation (8.4). When | € (x, »)| <M and |x —1|<e 


(9.9a) Wy (2,4) =C,O(1), and wy,(z, A) =C,,0 (A). 
When |&(— x, »)| SM and |x +1|<e, 
(9.9b) (2,4) =CmisO(1), and w),(z, 4) =Cy 1.0 (Aa). 


If z is bounded, *€6,, Gx ,2, 


E(x, »)|>M, and |&(—x,»)|>M, then 


(9.10a) Wy (2, 4) = Ay, (2, A) O (AA), 
(9.40b) Wm (2, A) = Wy, (2, A) O(A), 
where 


(9.10¢) Ay, (2, A) = phere hor (l—mani) p—o8 
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If z is not bounded, then the functions EZ, and E, defined by (9.7) may be un- 
bounded. In this case, the relations (9.10a) and (9.10b) are invalid. We therefore 
introduce a function }(z) which is continuous, which is nonvanishing for |z| >M 
and equal to 1 for jz|SM, and which is the “largest” function such that for each m 


(9.11) |b(z) A;,*(z, A) w,,(z,A)| <<N|v|-* and |b(z) A;,1(z, a) LAGER <N|y|-4, 


when *€6,,  Gy42 and when |é(x, »)|>M and |é(—x,v)|>M. Thus when the 
relations (9.11) apply, 


(9.12a) Wy, (2, A) 5) 
and 


Part III: Asymptotic Expansions for Solutions of the Given Equation 

10. Final hypotheses on ZY. Our objective is to determine the asymptotic 
behavior throughout Y of a pair of linearly independent solutions of the given 
equation (2.5). This will be done in the sections to follow. We first link the 
behavior of solutions of the given equation to that of solutions of the related 
equation by using the familiar method of variation of parameters. Its application 
to the given equation yields the integral equation 


z 


eal Wa (2) Wy (t) — Walt) wy (2) Q(t, 4) 
(10.1) u=w- if ae Asai u(t) dt, 


which is equivalent to the given equation in the sense that an analytic solution 
of either one is a solution of the other. In equation (10.1), w, and w, may be 
any pair of linearly independent solutions of the related equation, w may be any 
solution of the related equation, and z, may be any point in Y. The kernel 
in equation (10.1) is, of course, independent of the choice of w, and w, so that 
a unique solution w of the given equation is determined by specifying w and 2, 
in equation (10.1). 

The variables z and x and the parameters A and vy are always considered to 
fulfill the relations (9.4) and (9.2). Thus, given z and A, y and x are determined. 
In order to keep the investigation to a reasonable length and to avoid further 
notational complications, we make the assumption that 


(v) JargA| Sa. 
The second of the conditions (7.4) is thereby fulfilled when || is sufficiently 
large. There is no loss of generality in assuming that the first of these conditions 
is also satisfied. In what follows we fix 4, once for all, with | A| sufficiently large 
to fulfill all requirements placed upon it. 

The following lemma will be used in the discussion of solutions of (10.1). 


Lemma. // 
a) &@ is a closed, simply connected region (whose boundary may depend upon 
a complex parameter i), 
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b) f(z, 4) is continuous and bounded for ze and for |A|\>N, 

c) K(z, t, A) ts continuous for z€&f, for te, and for | A| SN, 

d) there exist a constant M, a point s EL, and a set of rectifiable analytic curves 
oining the points of to s such that for all z CAH 


f |K(z,6,a)| [dt] <<, 
then the solution of the integral equation 


g(z, A) =H, d) FAAS Ket, A) g(t,a)dt  (k>0) 


has the form 
g(z,4) =f(z,4) +A “*O(1), 

where O(1) denotes a function of z and ) which is uniformly bounded for zEL and 
for |A|>N. 
The hypotheses of the Lemma guarantee the uniform convergence of the series 
which is obtained by successive iteration of the integral equation, and from this 
the proof of the Lemma follows. 

For the Lemma to apply to the integral equation (10.1), it is necessary that 


the kernel in (10.1) be bounded. We now make assumptions on Y sufficient to 
establish the boundedness of the kernel. For each 4, let 


(10.2) Ry, = {2| %(2) EG, Grae} (m=0, +1, +2, +3). 


In this definition and in the sequel, when a region with subscript larger than 3 
or less than —3 appears, the region is understood to be the null set. Because 
of the overlapping of the regions G,, and also of the regions ©*, one can show 
that Z,,C By ~B mi (m=0, +1, +2). It follows from the hypothesis (v) 
and the conditions (7.4) that we lose no generality by assuming that #_,¢ #_, 
and &2,C#,. The regions &,, cover Y and each point of Z is in at least two 
such regions. Our final assumptions on @ are: 


(vi) Corresponding to each region &,,, there is a point z,, such that each point 
z€B,, may be joined to z,, by an analytic curve I lying in B,, and such that on I’, 
HE 1s monotone (monotone increasing from z,, to z if m is odd, monotone decreasing 
from 2,, to z tf m 1s even). 


(vil) The integrals 


i Q(t, A) dt 
AO 

B 
are uniformly bounded for z€Q, and for |a|>N. 


The function Q(z, A) first appears in equation (5.2). The function b(t) is defined 
by the conditions (9.11). Of course, if 2 is bounded, hypothesis (vil) 1s auto- 
matically fulfilled. 

There are certain shadow zones which must be excluded from the regions &,, 
in the discussion of subdominant solutions of the given equation. Their images 
on the x-plane adjoin the segment —1<x*x<0. Their size and presence depends 
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upon arg A. They are defined as follows: 


a = 12|\2e,,, R(oE)>N, |E(—x(2))|>N, 
(10.3) and the paths joining z to z,, all pass through 
the neighborhood of &,, where |&(—x)|<N}. 


The shadow zones are excluded for the reason that the subdominant solutions 
of the given equation undergo a change in their asymptotic forms as one crosses 
into them. The presence of these shadow zones (see Fig. 2 for an example) is a 
phenomenon which has no counterpart in the discussion of asymptotic solutions 
of differential equations having just one 
turning point. 

As we have remarked, we wish to de- 
termine the asymptotic structure through- 
out Y of a pair of linearly independent 
solutions of the given equation. We first 
prove that, for each m, there exists a 
solution u,, which is asymptotically re- 
presented by w,, in the region Z,,— F&,. 
We then derive connection formulas 
among the solutions u,, from which the behavior of the solutions u,, may be 
derived for all z€Y. Lastly, we discuss certain characteristic solutions, which 
may well be of special interest in applications of the theory. 


arg B=-290 


Fig. 2. A shadow zone 


11. Subdominant solutions when |&(+a, v)|>N. For each region %,,— F,, 
there is a solution of the given equation which is subdominant, exponentially 
small, for z€ (&,,— F%,,) and R(o €)>0. This solution is unique up to a constant 
factor. We now establish its existence and give its asymptotic properties. We 
assume for the moment that 7€(Z,—%,)°Z,.,. In equation (10.1) let 
W,—W,, and W=wW,,,,, let w=w,,, and let the path of integration be a curve I’ 
which originates at z,,. For brevity, we adopt the notation t=&(x(d),v) and 
& =&(x/(z),v). If the curve /’ should intersect one or both of the regions where 
|z(+x)| SN, then R(o&) <0. Thus, if we replace the subarcs of I’ on which 
|x(+)| SN by arcs on which |t(-- x)| =N, andif we note that on the remaining 
portion of J) R(oé)<R(oT), then we see that on the (new) path of integration 
exp [20(€—T)] is bounded. 

The choices just made determine a solution u,, of equations (10.1) and (2.5). 


If we adopt the abbreviations 
Un (2) = Ain (2) ty, (2, A) 
W,, (2) =A 6 (2) ©, (2, A), 


where A,,(z, A) and b(z) are defined by (9.10c) and (9.11), we can rewrite the 
integral equation (10.1) in the form 


(11.1) 


Bt, OMAN 
(14.2) Un (2) =m @) +f {om 2) Woes 2) Wt) Wavea 2) PEM} a 


In writing (11.2), we have used the evaluations (9.8) and (9.5) of W(w,,, @m41) 2): 
The relations (9.10), (9.12) and (41.1) imply that for |&(+-%,»)|>N and 
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2€ (By — Sm) 2m 41, both W,, and W,,41 are bounded. We have already observed 
that exp[2o0(—T)] is bounded on I. All functions in equation (11.2) are con- 
tinuous. In virtue of these facts and hypothesis (vi) of §10, the hypotheses of 
the Lemma are fulfilled with respect to equation (11.2). Thus for z€ (%,,—F%,) > 
Barr, \E(Ex(2),»|>N, and |A|>N, 


(11.3) U,,(2) = W,,(2) + A-* 7 0(A). 


We can show by analysis analogous to that above that this relation also holds. 
for 2€ (By— Sy) \Am—1. In terms of u,, and w,,, the result (11.3) takes the form 


mn A) O 
(11.4) ty 2, A) = ty (2, 2) + EDEN) (|| >). 


Since os SL) CBy~ Amer, (V4) holds for all 2€(4,—7,) such) that 
|E( &(4 a L x ( v)| > NG 
We can use the estimates (11.4) and (9.12) to obtain the the form of u,,(z, A) 


directly from the equation found by differentiation of (10.1). It is 


Am (z, 4) O (1) 
bye 


(11.5) Um (z, A) = w,,(z, A) 4 


where the prime denotes differentiation with respect to 2. 


12. The solutions u,,, when |§(+ a, v)|<N. There are two cases to be con- 
sidered: |&(x,y)| SN and |&(—x,v)|<N. For the moment we suppose that 
|E| SN. We recall that the solution u,, of the given equation and equation 
(10.1) is that solution determined by the choices w=w,, and zy, =z,, in (10.1), 
the path of integration being a curve [’.. We divide J’ in two parts Jj and J} 
on which |t|=N and |t|<N, respectively. As was done in § 11, we replace 
any part of Jj on which |7(—x)|<.N by an arc on which |t(—x)| =N. We now 
write the integral equation for w,, in the form 


ot) Up, (2) = Wy, (2) ane) K(z, t, 2) U,,(t) at el (2, t, A) u,, (t) dt, 


where 


K(z t i) _. Wm(2) Wty (t) — Wm (2) Wmt1 (2) : Qi, A) 
a W (Wm, Vm, t) tibia aie 


the upper or lower signs being used according as 1} CZ,,,, or Uae * 

For the sake of argument let ¢#%,,,,. On Ij, the behavior of w,, is then 
described by formula (41.4). Gene this result and the relations (9.9a), (9.10), 
and (9.12), we find that 


J B61) mul dt Ch Cuy | EWE OCN 26 NEO a 
é tb) Aetna 


? 


fhe, Tr 
Cm 


=f flow ree oun Mea 


Dy 


But (ot) 20 on J, except perhaps for an arc on which R (ot) = —N. Therefore 
by hypothesis (vii), 
(12.2) PRGA ae apne 


n+2 


1 
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From the estimate (9.9a) and the order relation Pag =0(A~*), we deduce that 
(12.3) WES AOA ae a) HO) ace Geekery 
The estimate ee enables us to fewit equation (12.1) a 


Um m ST za tn (t 
ae = Te ete OU) + [ Ke, t, A) tO ap 


The relation (12.3) shows that the hypotheses of the Lemma are fulfilled with 
respect to this equation; hence, 


(12.4) thy (2, A) = Wy (2,2) +Cy 2 (E(x, )| SN). 


m fet 2 


Since the order relation (9.9a) is independent of arg(z—/»/A), this result is valid 
in a full neighborhood of « =1, or of z=J/v/A, where || N and not just in that 
portion of Z,, where |é|<N. 

We now consider the case |é(—x)| <N. This situation is more complicated. 
If, for such z, Ji (o €)< —N, then the estimate (12.2) no longer holds; and therefore 
we cannot determine the behavior of u,, when |&(—x)|<N as we did above for 
|é| SN. We resolve this difficulty by showing that it is only an apparent one. 
Let us momentarily suppose that m= 2 and R(of)S—N. Then R(oé)=N 
in the neighborhood of &,,,4 where |€(—x)|<N. By the argument given in 
the next paragraph we are able to determine the behavior of w,,., in a full 
neighborhood of z = — |/v/A, or of x = —1, where |&(—x)|<.N. We are also able 
to determine the behavior of w,,,, in such a neighborhood. Similar considerations 
apply to values of m other than +2. Thus, we shall always be able to determine 
the behavior of two linearly independent solutions of the given equation in a 
full neighborhood of z= —J/v/A where |&(—x)| SN. 

As promised above, we now obtain the form of w,, if |&(—x)|<N and if, 
for such z, Wt (ao &)=—N. To do this we replace w,,(z) and w,,,,(z) in (12.1) by 
Wmig(e?™*z) and W,»~,(e*"*z), respectively, and similarly replace w,,(¢) and 
W,,+1(¢) for t such that |7(—x)|<N. In virtue of the identity (8.10), the analysis 


mx 


used in the derivation of (12.4) may be easily modified to yield the conclusion that 
(12.5) tg A) = ya (OF 2, A) men 8) 
when [= (=e, ”)| <N provided that for such z, R(o &) = —N. In (12.5) the upper 
signs are to be used if m<O, lower signs if m>0, and either if m0. As is the 
case when [El<N, we need not impose the restriction here that z le in &,, 
We need only suppose that the portion of the path of integration on ehicd 
|€(—x)|<N is rectifiable. Thus the result (12.5) is valid in a full neighborhood 
of z= —/A where |&(—x)|<N. 

We can use the estimates of w,, obtained thus far to find the form of w;,,(z, A) 
when |£|<WN or |&(—x)|<N directly from the equation obtained by differentia- 
tion of (10.1). The results are: 


(12.6a) Uy (2, A) (2, pea (Ee, v)| SN), 
(12.66) i, (2,A) = whsole***z, a) + Smee) (|£(— 2,9) SN), 


in neighborhoods of z=J/r/A and z=—Y/A, respectively. 
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We summarize our results in a theorem: 


Theorem 1. Under the hypotheses (i) through (vii) and corresponding to each 
region #,,, the given equation (2.5) has a solution u,, with the following asymptotic 
behavior when |A|>N- 

a) When 2€(Bn—F,) and |E(4%, v)|>N, u,, 1s exponentially large or small 
according as R(oé) <0 or R(a é)>0 and is described by the formulas (11.4) and 
(11.5). 

b) When |z—Yv]A|Se and |E(x,»)| SN or |z+)v]A| Se and |E(—x, »)| SN, 
u,, is oscillatory in character and is described by the relations (12.4), (12.6a) and 
(12.5), (12.6b), respectively, provided |/v/A and —vjA lie in Z,, and R(of) = —N. 

c) If y is real, then u,, 1s bounded and oscillatory for —15%S1,m=0, +1, +2. 


The regions &,, and %,, the variables x and &, and the parameters y and 
are defined by (10.2), (10.3), (9.4), (7.5), (9.2), and (6.1), respectively. The 
functions and constants appearing in the relations (11.4), (41.5), (12.4), (12.5) 
and (12.6) are defined in Parts I and II. The bounds for the error terms in all 
of the estimates for solutions of the given equation which have been derived 
are uniform in arg 4 for |argA| <a. This is so for the same reasons as those 
cited in § 8 for the uniform boundedness of the error terms in the expansions 
(8.8a) and (8.8b) with respect to arg ». 


13. Connection formulas. It remains to give a description of the solutions 
u,,(z, A) when |&(--x)|>N and when z is not confined to the specific region 
#,, as assumed in Theorem 1. The dependence relation among three solutions u,, 
u,, and u,,,, of the given equation may be written in the form 


(13.1) Uy, = Ay mUm + By mn Umt1 


The coefficients are, of course, given by the formulas 


W (Uy, Um, 2) W (Um, U4, 2) 
13.2 A a Bing —— 
( ; ) att W (Um, Um+1> 2) ; % W (Um, Um+1) 2) ; 

Since the given equation lacks a first derivative term, the Wronskians in 
(13.2) are all constants and may be evaluated at x =]/v/A, where x =1 and & =0. 
If each of the solutions u,, u,,, and u,,,, fulfills the relation (12.4), then we may 
conclude that 


W (Uy, uj) = W (wy, wi) + Cy COA"), 


where k,7 =l,m, or m+1. Thus by relation (9.8), 


W »°m Lomtv == 
Ana = Fie ey + Cmts CON a, 
(13.3) m> “m+1 


_ _W(Um, r%) toni —n—3/ 
PUNS SH TRA + Cy, Cre OW). 
Since v,=v, if k= j (mod 4), the relations (9.5) always suffice to give the values. 
of the Wronskians in (13.3). The constants C; are defined by (8.9a). 

The indices / and m for which the formulas (13.1) and (13.3) hold depend 
upon arg 9; hence, the pairs of linearly independent solutions of the given equa-- 
tion whose behavior can be determined over all of QZ, depend upon arg y. The 
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selection of these fundamental pairs of solutions may be made using Theorem 1 
and the connection formulas (13.1). We omit the details of the selection and 
state only the results. 


Theorem 2. For each value of A such that |arg 4| <a, the given equation (2.5) 
has a pair of linearly independent solutions whose asymptotic behavior over all of 


QZ, may be found from Theorem’ and the connection formulas (13.1); namely, 


if OSargyv<a U_g, Uy 
mSargyvSat+e ae U_g, Us 

the fundamental pair is lee 

= aig y= 0 U_1, Uy 
—nw—eSargvs—a Ug, Us. 


14. Characteristic solutions. If $(¢»—1) =k, a non-negative integer, 
Ym (0) =e HOM HY, (VER FD), 


where y,,(*, ¥) is a solution (6.1) of the transformed Weber equation (8.1) and 
H,,(x) is the Hermite polynomial of order k. The Hermite polynomials are even 
or odd according to the parity of k so that when 3(o»—1) =k, 


Vin (x, v) a cS 1)* Vinee (%, v) < 
These facts are reflected in-the behavior of the solutions of the given equation. 
If |arga|<e, 
W (tio, U2) = 4 (Yor Vas) +9? OO e¥ e#™* O(a") OC"). 


where the upper or lower signs are to be used according as arg y<0 or arg v>0. 
When v=2k+1, W(¥o, V2) =0. Thus, as shown by Lancer [4], for certain 
values of A the Wronskian W(u,, uw.) vanishes. These characteristic values are 
countable in number and are given by the formula 


A, = 2k +1+0(A,"%). 
At the characteristic values, 


ty (2, Ay) = (—1)" Wao (2, An) [1 +0," )], 


where the upper or lower signs are to be used according as arg (v(A,)) is negative 
or positive. We call uo(z, A,) a characteristic solution of the given equation. 
We can observe the unique behavior of a characteristic solution “)(z, A;,) for 
real z by referring to Theorem 1. 

Since | 9 (A, ®)|<N for values of z such that —15 RzS1 and|3z|<N|A,|4, 
a characteristic solution is oscillatory for z in a strip of width O(A;") centered 
on the interval [—1, 1]. In certain regions outside of such a strip, a characteristic 
solution is subdominant and decreases exponentially to zero as |z|—> oo. These 
regions are sectors of width }2—e with vertices at z= +/A,/v(A,) and which 
are centered on the real axis. In particular, a characteristic solution is uniformly 
bounded for real z. Only characteristic solutions and their multiples have this 


property. 
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On S ingularities of Solutions 
of Partial Differential Equations in Three Variables 


ERWIN KREYSZIG 


Communicated by S. BERGMAN 


1. Introduction 


In the case of linear partial differential equations in two independent variables 
there exist various operators which transform analytic functions of a complex 
variable into complex or real solutions of these differential equations. Such 
mappings of the algebra of analytic functions into the linear space of the afore- 
mentioned solutions can be considered as a natural generalization of the well- 
known relation between analytic functions and harmonic functions in two 
variables. Some of these operators preserve many properties of the analytic 
functions involved. BERGMAN operators [/|* possess this important feature. 
In this way theorems on analytic functions can be “‘translated”’ into theorems on 
solutions of partial differential equations in two variables, 7.e. the theory of these 
solutions can be based on the theory of analytic functions in a systematic manner. 

In the case of elliptic partial differential equations in three variables with 
analytic coefficients there exist similar possibilities. For example, solutions of 
LAPLACE’S equation 


Av Cae ale t Ux, xp t On 4 =O 


can be obtained from functions f(w,¢) of two complex variables ¢ =e" and 


U =X, 4 a (Gs 1 Hy) € = % (%g — tq) 64) 


by means of the WHITTAKER representation [2] 


22 
V(X, X2, Xs) = se i f(u,e") dt = say fieoer ac. 
0 ¢|=1 
Using this representation, BERGMAN [3—6] investigated harmonic functions v 
corresponding to different classes of “associated functions” f, and solutions of 
equation (1.1), below. 
The present paper concerns partial differential equations in three variables 


of the type 


(1.1) Ux f Vy 2,1 Ux x A(X, X%3) VU = 0. 
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We mention that our investigation can be extended to certain classes of more 
general partial differential equations in three variables with analytic coefficients 
whose solutions exhibit a similar behavior. The basic tools to show this are 
Bergman operators. 

Before stating in Section 3 the problem considered in this paper, we introduce 
suitable notation. 

2. Representation of solutions of (1.1) 


We assume that the coefficient «(x,, %3) in (4.1) is an analytic function. If 
we introduce new variables 


X=kXy. Zee ele 2 Se es 
and use 
a peek? O20 ee 
OZ Ge, ige OLE On. Onn 
equation (1.1) takes the form 
(2.1) Wx x —Wzzx+F(Z,2*)w =0 


where w(X, Z, Z*) =v(%,, %, %3) and F(Z, Z*) =a (Xp, Xs). 
We assume that F(Z, Z*) is an entire function. 


Theorem 1 (BERGMAN [7]). Let f,(Z) be an analytic function which is regular 
in a simply-connected domain DY containing the origin. Let the functions 


Gaep\Z,2"), SW — 2), w =O, 1... aly Deh aed eee eee 


: Gn n0 = O, 
be solutions of the system 


e NS ) n,S,P+ 
(2.2) (4 a Os) (s = 2) (s 1) Qn, s+2.p ret Test - (Insp 5 Ovo) a =O 


with the property 
(2.3) Grp (Ss 0) =. 


Then the functions 


(2-4) [bn] co ZZ, Zy1 
== ane (dor hi (Z) ts ps Tn nae, p\fae) A J uy J ips (Z,) aZ, “ee aZ,) 
re p=y 00 0 


are solutions of (2.1) which are regular for ZED, | Xx 
symbols. 


<oo. 6,, etc. are Kronecker 


Solutions which are real for real values of the original variables 
=X, %=—1(Z4+2Z*), x=Z—Z* 
will be termed real solutions and denoted by capital letters. Obviously, 
12:5) WAX, Z,2*; ft.) = (w,(X, 2,2"; 2), Oe) 


is a real solution corresponding to the complex solution (Xi, ZZ Ne Ch (2A 
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Remark. We shall consider solutions of (2.1). In this connection it is important 
to study what properties of those solutions are completely independent of the 
special form of the coefficient F or depend only on certain properties of F. Our 
aim is to work out tools for investigating those properties. In some cases of 
differential equations (2.1) our problem can be reduced to a similar one concerning 
certain differential equations in two variables. Then we may derive the desired 
relations by means of previous results [10]. 


3. The problem considered in this paper 

The representation (2.4) is helpful in considering coefficient problems, 7.¢., in 
determining properties of solutions of (2.1) in terms of the coefficients of the 
corresponding power series development. Applications of this type have been 
discussed by BERGMAN [7]. 

The properties of the operator P,(7,) are similar to those of the so-called 
Bergman operator of the first kind introduced by BERGMAN [8] for generating 
solutions of partial differential equations in two variables. Indeed, there exist 
simple relations between the solution w, (X, Z, Z*; f,,) and its “associated function” 
f,(Z). From (2.3) and (2.4) we obtain 


(3.1) @,(X,Z,0; f,) =X" f,(Z), 

(3.1 b) Wy, (X,0,2*; fn) =X” f, (0). 

In the case of the corresponding real solution (2.5) we thus have 
(3.2a) W,(X,Z, 03 fy) =3 X" (7,(Z) +74(0)), 
(3.2b) W(X, 0,2*5 fn) = 3X" (fn (0) + fn(2Z*)) 


We assume that /,,(0) is real. 


The relations (3.1) and (3.2) are basic in connection with coefficient problems. 
First of all, instead of w, (X,Z,Z*) we may consider solutions of (2.1) of the form 


N 
(3-3) w(X,Z,Z*\ = > w,(X,2,2*;7,) 
n=0 
where the associated functions /)(Z), ..., fy(Z) are regular in a common domain 
of the above type. Let 
N oo lo) 
(3.4) w(X,Z,Z*) =) > Ds Peale Le. 
n=0 m=0 p=0 
Then from (3.1) it follows that 
(3.5) PZ) ys oe Oe oer 
m=0 


Hence the domain of regularity, the nature and location of singularities and other 
properties of w(X, Z, Z*) can be determined by means of the NV +1 subsequences 
Carine =O) ty WN of the cociicients—im the corresponding 
development (3.4). Cf. BERGMAN [7]. 

It may be expected that properties of w(X, Z, Z*) can be obtained in terms 
of other subsequences (Ry 5), 7% =0,1,-.-, 7=0,1,.--, N, #>0 and fixed, of 
the coefficients of w(X,Z,Z*). The present paper is devoted to this problem. 


Ae 
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4. Definition of the functions a, (X, Z) 


In order to investigate the problem just stated it suffices to establish relations 
between the N +1 subsequences (imo), 77 =0,1,-.-, #=0,1,...,N, and the 
nm +1 subsequences (Ry mp), 7% =0, 1,..., 47 =0, 4,..., N, where # is an arbitrary 
natural number. Relations of this type may then form the basis for a more 
general investigation of coefficient problems (cf. Sections 7 and 8). For this 
purpose we shall now introduce auxiliary functions a,(X,Z) in such a manner 
that relations between the above subsequences of coefficients are equivalent to 
relations between the functions a,(X, Z). 

Lemma 2. Let w(X,Z,Z*) be an arbitrary solution of (2.1) which can be 
represented in the form (3.4). Let the coefficient F(Z, Z*) be an entire function which 
has the representation 


(4.1) KZ) > Cpe oleae 
m=0 p=0 
and let 
(4.2) ¢p(Z) = 2 mp 2”, p=0,A, 


Then the functions 


N fee) 
(4.3) ap(X, Z) — De Raspes p = 0,1, OTN) 


n=0 m=0 
are solutions of the following system of partial differential equations: 


0a Cra Pact 
p p— 
(4.4) Eee ee >) p-s-1 4s = 9, p=1, 2, -.<. 
s=0 


Proof. From (4.3) and (3.4) it follows that w(X, Z, Z*) can be represented 
in the form 


(4.5) w(XyZ,2*) =) apy eZ ye 
p=0 

From (4.1) and (4.2) we obtain 

(4.6) VOTE) = 2¢5(2) Le, 


co 


0 ap gk —1 , Oa eP-1 1 , £ 
gee 7) 2a es eee 
=: p=0 


Dax 
p=0 

In order that this equation be satisfied for every value of Z* the coefficient of 
each power of Z* must vanish. This yields (4.4). 

The coefficients in the development (4.3) of each function a,(X,Z) form just 
one of the above-mentioned sets of N +1 subsequences of the coefficients in the 
representation (3.4) of w(X,Z,Z*). Relations between the functions a (xa) 
are thus equivalent to the desired relations between those sets of uieaaneeee’ 
These relations will be obtained from (4.4). The first 7 equations (4.4) involve 
only the functions a (X, Z), a,(X,Z),..., a,(X, Z). 


s=0 
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5. Singularities of a,(X,Z), r>0, and their relations to those of ag(X, Z) 


We consider any one of the functions (4.3), say a,(X,Z), 7>0 and fixed, 
corresponding to a solution w(X,Z,Z*) of (2.1). Since we can determine the 
domain of regularity of w(X,Z,Z*) from the subsequences (Ry mo), 7 =0, 1, ---; 
=O, Ute. Neelef. (3.4)), ¢.e. from the function a)(X, Z), it is natural to inquire 
how the singularities of the above function a,(X,Z) are related to those of the 
function a)(X, Z). 

First of all, the distribution of possible singularities of a,(X, Z) is relatively 
simple. If a,(X, Z) is singular at a point, say (Xy, Z,), then it is singular at every 
point of the plane Z=Z,, as follows from (4.3). That is, while in general the 
singularities of an analytic function of two complex variables lie on arbitrary 
analytic surfaces, the functions under consideration may be singular only on 
certain planes Z =const. 


Lemma 3. I} the function F(Z, Z*) in (2.1) is entire and if a,(X,Z), r>0, 
is singular at the points of a plane Z =Z, lying in a simply connected domain D 
containing the origin, then ay(X, Z) is singular at every point (X,Z,) in Z. 

Proof. Since the coefficient /(Z, Z*) in (2.1) is assumed to be an entire func- 
tion, all the functions c,(Z) in the differential equations (4.4) are entire. Con- 
sequently, if a,(X, Z) is singular on the plane Z =Zj, at least one of the functions 
a,(X, Z), y—1=s=O, appearing in (4.4) must be singular on this plane. By 
repeating this argument the statement follows. 

Combining Theorem 1 and Lemma 3, we obtain the following result. 


Theorem 4. Let the coefficient F(Z, Z*) of the partial differential equation (2.1) 
be an entire function. If a solution w(X,Z, Z*) (cf. (3.4)) of (2.4) ts regular for 
all values of Zin a simply connected domain G containing Z =0, then the correspond- 
ing functions a,(X,Z), p=1,2,..., defined by (4.3), are regular for all values 
ZEDand |X|<oco. 


6. The domain of regularity of solutions of (2.1) 
in terms of arbitrary subsequences (k,.,,), Y > 0 and fixed 


The properties of the functions a,(X,Z), a,(X,Z),... stated in Theorem 4 
result from the regularity of the corresponding solution w(X, Z, Z*) of (2.1) ina 
certain domain. Unfortunately the condition that these functions have the above 
property is only necessary but not sufficient for the regularity of w(X, Z, Z*) 
in that domain. Indeed, let us consider the simplest of the equations (4.4), 7.e. 

0a, 0 Ay 
(6.1) mo oer ro) Ao- 
From (3.14), (3.5), and (4.3) we obtain the representation 


(6.2) a9(X, Z) = X" fy (Z). 


Hence if a)(X,Z) has a plane of poles resulting from /,(Z) or /,(Z) and if this 
plane coincides with a plane of zeros of cy(Z) of sufficiently high order, then 
a,(X, Z) is regular at the points of this plane. By considering the other equa- 
tions (4.4) in a similar manner we thus obtain the following intermediate result. 


156 Erwin KREYSZIG: 


Lemma 5. Suppose that the coefficient F(Z,Z*) of the partial differential equa- 
tion (2.1) is an entire function. Let Z=Z, be a plane of poles of order s of a)(X, Z) 
resulting from the associated functions f,(Z) or f,(Z). Let Z=Z, be a zero of order 
m,=s of ¢;(Z),7=0,1,...,7—1. Then the functions a,(X,Z),...,a,(X,Z) are 
regular on the plane Z =Z,. 

The coincidence of poles of a,(X,Z) with zeros of c;(Z) is the only reason 
that a plane of poles of a)(X,Z) might correspond to a plane of regular points 
of a,(X,Z), #>0. This can be seen from the various equations (4.4). If a)(X, Z) 
has a plane of poles of the order s, say the plane Z=Z , and if c (Zo) +=0 then 
a,(X,Z) has at Z =Z, a plane of poles of the order s — 1 (or a plane of logarithmic 
singularities if s=1), as follows from (6.1). In a similar way we find from the 
equation (4.4), f=2, that if a)(X, Z) has a plane of poles Z =Z, and cy (Zo) +0, 
c,(Z))=- 0 then a,(X,Z) is singular at Z=Z,. Indeed, the functions involved 
in that equation have singularities of different orders which cannot cancel each 
other. Similar arguments hold with respect to essential singularities. The other 
functions a,(X,Z) can be treated in a similar way. Using Theorem 1 we thus obtain 

Theorem 6. Suppose that the coefficient F(Z, Z*) of (2.1) is an entire function. 
Let w(X,Z, Z*) be a solution of (2.1) which can be represented in the form (3.4). 
Suppose that in this representation the coefficients k,,, %=0,1,...,N, m=O, 
1,...,7>0 and fixed, are such that a,(X,Z) (cf. (4.3)) is regular in a simply-con- 
nected domain D which includes the origin but no zeros of th functions c;(Z), 
7=0,1,...,.7—1. Then w(X,Z,Z*) is regular in Bx(|X|<oe). 


7. Singularities of solutions of (2.1) in terms of arbitrary subsequences of 
coefficients in the corresponding power series development 

By means of Theorem 1 various results on analytic functions can be “‘trans- 
lated’’ into theorems on solutions of (2.1). These theorems involve only the 
subsequences (R,, 9), 7 =0,1,...,N, m=O, 1, ..., of the coefficients in the series 
development (3.4). Cf. BERGMAN [7]. This fact is a consequence of the relation 
(3.2a) between the solutions w(X,Z,Z*) of (2.4) and their associated functions 
fo(Z), ..-, fy(Z). We shall now show how to apply the results in the previous 
sections to obtain more general theorems involving arbitrary subsequences (h,,»,,); 
w==0,41,...,N, m=0,1,...3 7>0 and fixed. of the coefficients in (3.4). It 
suffices to outline this procedure by means of some illustrative examples; the 
basic idea of deriving such theorems on solutions of (2.1) from theorems on 
analytic functions is similar in each particular case. 

Let us generalize a theorem of Borer on rational analytic functions. 


Theorem 7, Let w(X,Z,Z*) be a solution of (2.1) which can be represented 
am the form (3.4). Suppose that the function F(Z, Z*) in (2.1) is entire, and that 
the functions c,(Z), 7 =0,1,...,7—14, (cf. (4.1), (4.2)) ave different from zero in a 
simply connected domain B of the Z-plane which includes the origin. Then w (Xi L525) 
has in B x(|X|<co) at most finitely many planes of singularities (which are 
planes of poles of order not exceeding s +r and (for all Z*==0) logarithmic branch 
planes) tf and only if all but at most a finite number of the determinants 


(7.1) Dain =| tive Goes ta ae 
M = 0,41, i015 S10; eee ee Omanaaexed. 
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with column vectors 


Cpr 
IBM {Ry, mA-7—1,7 Ry, m+i,4r°°* Sa tare gh 


vanish. 

Proof. The coefficients k,,,, appearing in the determinants (7.1) are just 
those which figure in the representation (4.3) of the function a,(X,Z), 7>0 and 
fixed. We may represent a,(X,Z) in the form 


N 

(7.2) AOR YS be AVA 
n=0 

where 

(7.3) hy (Z) = Pam 2”, 
m=0 


Wi ON eeeeN » R= Oand. fixed. 


For each fixed value n =v, O0Xv<N, the condition involving the determinant 
(7.1) is necessary and sufficient in order that h,(Z) be a rational function which 
has poles of order not exceeding s. Cf. E. Borer [11]. Hence the only singularities 
of a,(X,Z) are the planes of poles of order not exceeding s which correspond to 


the poles of the functions h,,(Z), 7 =0,1,...,.N. In consequence of Lemma 3 
the function a)(X, Z) is singular on these planes. Since the domain # does not 
include zeros of the functions c;(Z),7 =0, 1, ...,7—1, these planes of singularities 


are the only ones of this function in #, which follows from Theorem 6. The 
nature of the singularities of a)(X,Z) can be determined by means of (4.4). 
One finds that the planes of singularities of a)(X, Z) are planes of poles of orders 
less than or equal to s +7. In consequence of (3.2a) these singularities correspond 
tommpoles Or thesassociated functions j(Z),.2—0,A, ....V, ofthe ‘solution 
w(X,Z, Z*) under consideration. Since the coefficient F(Z, Z*) of (2.1) is assumed 
to be an entire function, all the functions q,,,,(Z7, 2*) occurring in the representa- 
tion (2.4) are entire. The integrals in (2.4) can be transformed into single integrals. 
By means of repeated integration by parts we obtain 


st Beet oe, pl = ae 
(7.4) SSS tn(Zp)4Zp_1...dZ, = 21d; 2? Jef, (f) at 
00 0 j=0 0 
where 0;,,7=0,1,...,—1, are constants; cf. BERGMAN [7]. If we substitute 


(7.4) into (2.4) and use the above result on the singularities of a)(X, Z), the state- 
ment on the type of singularities of w(X, Z, Z*) in # follows. 


8. Singularities of real solutions on the boundary of the domain of convergence 
of their power series development 

Let us consider another application of the results in the preceding sections 
of this paper and generalize a theorem by BERGMAN [7]. 


Theorem 8. Let 


N N co co 
W(X, Z, Z*) ==> UX 2 3'1,,) =) a yas IG ep LZ: Ge 
n=0 n=0 m=0 p=0 
8.1 i 
( ees =Kipm 
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be a veal solution of (2.1) whose associated functions },(Z), 7=0, 1,--- N, are 
veal at Z =O; cf. (3.3) and (2.5). Let F(Z,Z*) be entire; ct. (2.1). Suppose that 


it 
lim 1 sup Kye =O = const; «OA aoe, 


for an arbitrary ie value r>0. Suppose further that c,(Z)=-0 for all values of 
|Z| So and 4=0,1,...,.7—41. Then W(X, Z,2*) 18 sinedler On LL oe a), 
and only 1f 


. : Ss 7 | \se 2 
(8.2) lim sup (| > (5) Boor \P=2 
Proof. We introduce the functions 
N co 
(8.3) A,(X,Z) = % Pt Gg Ds A OES aks 
n=0 m=0 
Then 


to 


Wik; Z727) ae PO 4 CH 
By substituting this representation into (2.1) we see that the functions 4,(X, Z) 
are related by a system of differential equations similar to (4.4). From this 
system it follows that the results of the preceding sections are also valid for real 
solutions of (2.1) for which each /,,(0), 7 =0, 1, ..., N, is real. Using these results 
the above statement is a consequence of the following lemma (cf. BlEBERBACH [72], 
Paso) eLet Be 
7 (2) = pa R, 2° 

s=0 

be an analytic function regular at the origin. ao the point 2=2, =e" 0 


the circle of convergence |z|=0, 0 = lim n sup | Rk, \s is a singular point of f(z) i 
and only if 


where 
vooy Leo 
o=0 


Indeed, from this lemma it follows that the above assumption on the sequences 
(Apo t =0), 45.53, Ny mm =04, .. 70 and fixed: isrthe necessary and suf- 
ficient condition for the function A,(X,Z) to be regular for |Z|<e, |X|<oo 
and singular on Z=Z,, n=0,1,...,.N, |X|<oo. Hence, by Lemma 3, the 
function Ay(X,Z) is singular at those points. Because of the assumptions on 
c;(Z), 7=0,1,...,~—1, it follows from Theorem 6 that these are the only 
singularities of Ay(X,Z) in |Z|<o, |X|<co. From this, (3.3), and (3.2a) the 
statement follows. 


9. Final remark 


So far we have assumed that the coefficient F(Z, Z*) of the partial differential 
equation (2.1) is an entire function. This is essential to the proofs of all the 
preceding results. Indeed, the situation changes if we permit F (Z, Z*) to have 
singularities. We confine ourselves to a brief remark in this direction. 
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Suppose that F(Z, 0) =cy(Z) (cf. (4.2)) has a pole at a point Z=Z,. Let 
w(X,Z,Z*) be a solution of (2.1) which can be represented in the form (3.4) 
and has the property that w(X, Z), 0) =a)(X, Z))=-0. Then the corresponding 
function a,(X, Z) is singular on Z =Z,. However, if all the associated functions 
fo(Z), ..-, fy(Z) of w(X, Z,Z*) have zeros of sufficiently high order at Z=Z,, 
then a,(X, Z) is regular in the plane Z =Z,. This and similar statements regarding 
the functions a,(X,Z), $=1, 2,..., follow from the differential equations (4.4). 

Lemma 3 will no longer be true if /(Z, Z*) is not entire. A plane of singularities 
of a,(X,Z) may now result from either a)(X, Z) or F(Z, Z*). The situation with 
respect to other results in this paper is similar. 


Theorem 9. Suppose that the coefficient F in (2.1) ts a function of Z only 
and has a singularity at a point Z=Z,. Then for every solution w(X,Z, Z*) of 
(2.1) which can be represented in the form (3.4) and which has the property that 
w(X,Z,,0)=0, the corresponding functions a,(X, Z), p=1, 2,..., are singular 
at each point of the plane Z=Z, in the domain of existence of the solution. 

This result follows from (4.4). 
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Sonic Limit Singularities. Part I: General Theory 
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Communicated by HILDA GEIRINGER 


The sonic limit line was first identified by H1LpA GEIRINGER [1]. In the 
same paper she also gave a classification of limit singularities for M=1 which 
was later adopted by the present authors [2]. In analogy with the treatment of 
limit singularities in the supersonic region [1], this classification was based on 
the vanishing of y; and p, at points in the hodograph, where y is the stream- 
function and & and y are the characteristic variables. 

Further investigation has shown that this classification is not adequate since 
the &- and y-directions become indistinct at the sonic circle and consequently 
the vanishing of y; and y, at a sonic point may depend on the path of approach 
in the hodograph plane. In the present paper we suggest a new classification 
(Sections 3—6) which is based on the vanishing of y, and ws, where g and # 
are the hodograph variables, and we discuss the geometry in each case without, 
however, reproducing the results derived in [2] and [5]. To this end we first 
show (Section 2), under the assumption that y is a regular function of the hodo- 
graph variables, that a sonic limit point cannot be an isolated limit point, 7.e. it 
must either lie on a sonic limit line or ona supersonic limit line. Examples of 
flows exhibiting these singularities will be presented in Part II*. 


1. Introduction 
First we recall briefly some of the basic equations stated in [2]. 
In the hodograph variables g, # the streamfunction y satisfies Chaplygin’s 
equation 


(1) VYaq a Yoo Closter 


where 


a yt 1—@° 
For any non-constant solution y(g, #) of this equation the physical coordinates 
x(q,8) and y(q, 8) are determined by 


dx = {[—sin 8 yy + a 


24 


cos Ppp| dq + [gcosd y, — sin d po] ao} 


f[2—4 


dy = 97 {C088 Ye + 2 sin 9 y,| dq ain [qsin } yp, + cos B pe | a9\ 


* GEIRINGER [6] has pointed out the occurrence of corresponding singularities in 
plasticity theory, where the present results should have counterparts. 
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where 9 is the density, 


1 
Gag i 
The Jacobian of the transformation from the hodograph to the physical 
plane is a(x, ) 4 
aa a(q, 8) F oq? Lope, = (1) 5 | 


In the supersonic region (M>1) it is possible to introduce the characteristic 
coordinates &, 7 of (1) by 


ire The Ati dds 
2¢= f ee ee 2n = f dg 6 


and write the Jacobian in the form 


etal eh ay 
Breakdown of the transformation will occur at points in the hodograph plane 
where J =0, and their images are called limit points. Curves consisting of such 
image points are called limit lines. 

Clearly when M>1 the vanishing of J implies that either y;, or y,, or both 
vanish. A curve in the physical plane along which y;=0, p,=- 0 is called a limit 
line #,, and one along which wp; +0, y,=0 is called a limit line Z.* The point 
of intersection of two such curves is called a double limit point. 


2. Non-isolation of the sonic limit point 


A question which immediately suggests itself is whether limit singularities 
at sonic points can be isolated or must lie on a line of singularities. It is well 
known that limit points in the subsonic region must be isolated, whereas ordinary 
and double limit points occurring in the supersonic region must lie on supersonic 
limit lines. Under the assumption that w is regular in the hodograph plane, the 
main result (Theorems 1—4) proved in this section is that every sonic limit 
point must lie on some kind of limit line. This line either consists exclusively 
of sonic limit points — and is then called a sonic limit line — or it contains 
supersonic limit points, in which case the point under consideration is a sonic 
point on a supersonic linut line. An investigation of the geometry in the neighbor- 
hood of these points (which differs in several respects from that for supersonic 
limit points) will be the subject of Sections 3—6. 

Throughout the present discussion we shall assume that p is a regular func- 
tion of g and ®. It then follows that the Jacobian J = — [q? y7 —(M?—1) y5 J/e2@? 
vanishes for M=1 (¢=gq,) if and only if y,=0. 

A point at which M=1, y,=0, ps0 will be called an ordinary sonic limit 
point, and one at which M=1, y,=0, ys=0 will be called a higher-order sonic 
limit point, It follows immediately from Chaplygin’s equation that at a sonic 
limit point y,,= 0. 

A line along which M=1, y,=0 is called a sonic limit line, #,. Along such 
a line Y= Y_o= Yq00= Ygo09= °° = 0, aNd Yog=Yyq9= Yqq99—= Yaqoo0— "=O. 


*x Their images in the hodograph plane will be denoted by /, and /,, respectively. 
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Theorem 1. Every ordinary sonic limit point les on a sonic or supersonic 
limit line. 
Proof. At an ordinary sonic limit pomt M=1, y,=0, po= 0. Hence with. 


F=—o%qI: 1 1 2 
F=y; (ce — zz) ¥3 = 0 
{ 1 peadearel 43 Ute lyy Bee 
E= 20,0, 20Poa (a ra Yo aa ta i a @ eo? 


1 1 
Fy = 2 Vo0 — 2V0 Yoo ( = =a =U; 


and by the implicit function theorem there exists a curve g=/(#) through the 
point * on which # = 0 with the direction 


EE Mig ttle af = 
Zak ied ES 0,  @ = COnse. 

For higher-order sonic limit points the non-isolation cannot be deduced from 
the implicit function theorem since F, = k= 0 (and if F, ,=0 then Fo — FE, F59=0)- 
To handle this case we consider the sign of F*= y,+ ()M2— 1 /q) pp on certain 
families of curves in the hodograph. 


Lemma. If at a sonic point all derivatives of p with respect to g and & vanish. 
up to and including the n order, then all derivatives of order n+-1 vanish, except 
possibly Yo9...0 and Poo... 

Proof. This follows by successive differentiation of Chaplygin’s equation (1). 

Thus the discussion of the general higher-order sonic limit point where ale 
derivatives of y up to and including the 2 order vanish can be broken down 
into the following two cases: 

(a) The first (in the order of decreasing g-subscripts) non-vanishing derivative 
of w of (n+ 1)* order is pp. 9. 


1+n 
(b) The first non-vanishing derivative of y is p, 9... 9. 
l+n 


In the second case, we have to consider » odd and m even separately. 


Theorem 2. In case (a) the sonic limit point is not isolated. 


CALE Let — Yoo...9=R+0. Then expanding in a neighborhood of the 
point *%, a 


Py = 0 (5g, 68)"™, 7 = 0 (dq, 60)?" 
Yo — h(O0)" + O09, OUT Ye =F (08) 010g, 09) 7% 


M2— 4 y+1 
eae Su é O 2 


ate rs Strictly speaking we should write “image point in the hodograph’’. The abbrev- 
lation is used whenever the meaning is clear. 


plage 69)” denotes a regular function of dg and 6% with terms (0q)* (68)! where- 
Pl en. 
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so that M2 — * 
Bm yh — Sip vb = — 71 Bbq) (68) + 0 (59, 08)8"* 
t 
Thus on the straight line 6g=m 60 the principal part of F behaves as follows 
li ytig 
(oyrrH = wi Rem. 


Thus for sufficiently small 69 +0 we have that F = 0 according as 6g =m 60 S.0 
This means that F takes on positive and negative values for various choices of 
the parameter m as shown in Figure 1. 


F<0 


F<0 


/4=1 


Fig. 1. The sign of F in case (a) Fig. 2. Vanishing of F in case (b) with n odd 


Now Ff is a continuous function on any small circle C about the point, and 
hence there exist points on C such that /=0. Thus F must vanish on a curve 
through the point, and this curve cannot lie in the subsonic region, since limit 
points are isolated there. 


Theorem 3. In case (b) with n odd the sonic limit point 1s a sonic point on a 
supersonic limut line. 


Proof. The fact that the first non-vanishing derivative of y, with respect 
to # is of odd order implies by Taylor’s theorem that yw, changes sign on g=4q;. 
This implies that the continuous function Ft= y,-+ (//M?—1/¢) ps tends to values 
with opposite signs at the ends of any sufficiently small “‘semicircle’’ C about 
the point on the supersonic side. (See Figure 2.) Hence F*=0 for some super- 
sonic point on this semicircle. The same result holds for F-. 


Lemma. Consider case (b) with n even and the derivatives of py of order 1+n+k 
at the point. Then for k=1, 2,...,3” 


0, pb<n — 2k, 
: Ce p=n—2k, 


n 


(Yaq...4 aes 

n+k—p p 

and for k=1,2,...,4n—1, 
0, p<n—2k—1, 


Oo a me! rs ~ 
one pa K(P) Poo...0 > Di aN 
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where, apart from k, K(k) depends only on the first derivative of (1—M?)/q? at 
M=1. That is, the first nonvanishing derivative of order 1+-n-+-k 1m the order 
of decreasing g-subscripts contains 1+-3k differentiations with respect to q and 1s 
a multiple of Wqg.,.9 only. 

Proof. Successive differentiation of Chaplygin’s equation (1) shows that all 
derivatives of y of order 1-++ 2+ k which involve more than 1+ 3 differentiations 
with respect to g can be expressed in terms of derivatives of order m and lower, 
and derivatives of order 1+ other than y,5,..5 and Ws9,..9. By the previous 


lemma, then, these derivatives of order 1+”-+k must vanish at the point. 

Similarly, we see that the derivative of order 1-+-”-+ which involves pre- 
cisely 113 differentiations with respect to g can be expressed in terms of 
Yqo...9 28 indicated. 

It is important to note that the same results (including the function kK (A), 
are obtained if in Chaplygin’s equation the coefficient (1— M?)/q? is replaced by 
the first term in its Taylor expansion at M = 1, and the coefficient of , is replaced 
by zero. Thus instead of Chaplygin’s equation we may consider Tricomi’s equation 


ae. y a tS 
qi 


whenever we are concerned solely with the derivatives indicated in the lemma. 


(3) Pyy —V Yee = 9, oq, 


Theorem 4. In case (b) with n even the sonic limit point is a sonic point on 
a supersonic limit line. 


Proof. From the lemma, the sum of the leading terms (one for each k) in the 
Taylor expansion of F*=y,-+(\/M?—1/¢) ys on the semi-cubical parabola 
(Viy+ 1)/q? 6g)! = 3p 60 is a multiple of W,0...9-* Now, from the remark above 


——— 


ee n 


* From the lemma, the principal parts of the Taylor expansions of yy and wp 
about the point reduce to 


a n—2k 
Tce Mary 2 manne ) K() Yq0...0(59)°* (60)"—>*, 
r iF n 
1 
vo= a7 (1) Yooo...0(84) (69)"—1 + 
Ley. fi +h 
VW c . ‘ ai vid 
+ 2) Geear (se 41) BA) Yee...0(89)* (30-27, 


k=1 


n 
Hence the principal part of F* on (j/(y + 1)/q?6q)!=3p 68 is 


k=1n 
1 5 n-k er ae 
jets — JN ee NR ae. ie = = 
1E: Seay 3k )tae Me 4)... 8-5-2) (5-2 | + 


k=l ( 
k=4hn-1 
We ahi XI 1 n+k 3 .\2k41 
eames GAT (gh-41) 034-1) ---8-5-23 (59) |} vue. 00) 


However, it is difficult to discuss changes in sign of F* from this representation. 
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it is clear that the same leading terms result if Tricomi’s equation (3) is used 
and G*=y,+y'y, is expanded on yi= 3x. Moreover, these leading terms 
remain unaltered if we construct any other solution of Tricomi’s equation which 


gives the same Py 20.3 , at the point considered and has the same earlier derivatives 
Fe 
equal to zero there. 


Now it is well known [3] that the general solution of Tricomi’s equation is 
given by 


= a is aaa 


+ fils flee tor a]eve morte 


where t(*) == (x; 0), v@)= sl, es 


and in view of the above results we select the solution determined by the Cauchy 
data 


namely 


3) 2 4 ( 1 = 
y= ss fis Bes, y(2t— 4) 8 — 0) dt. 
Tae 3 


GO)sar iY) Hairs 2h 


0 
«(4 +o(1 + +n) Bi) +inp i 
we now show that G*(f) changes sign for certain values of #. Assume that 
=—-y ~>> 0. Then for #=0 we have 
ea ERO 


5(14 — a) 8 dt*. 


ea (N= "ie few -bdt=ax">0, 


and for the characteristic =1 of the Tricomi equation 
1 
ae n— | n—i pars 
2 nt|"-* (1 —t)- eat 
0 
P(g) P(e +3) 

* This integral is of course identical to the corresponding polynomial appearing 

in the previous footnote. 


G* (1) = ax” 


1(2+3n)>0, 
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and re) 1 : 3 
e wn V's) 5n-1 gle: LIN pis (hed ya 
G(t) = a2 a2 f(2e—30@ —a]eta—9 
0 
ice n Pe rin —*) qn —2 2 n) <0 
TEP ee te 
Similarly we have for the other characteristic f= — 1 
1 
n P(3) ney 3 ios 33 | —4 SP \ tae gt 
G*(—1) = a3" Farge 2 | [a —(1 +5 0) (2 1) Font —yrt 
0 
Now letting ‘=1—s 
1 
n I'(3) n— 3 — De ec \ ee ae 
Gao ah 2 ee +(1+3n)@s—1) #4 n (1 —s)~hs"-8ds 


0 
so that G*(— 1) =G* (1) and hence G-(— 1) > 0 and G*(—1) <0. 
iG. 
G=0 ho G=0 


a 
ey 8 =y/ 
(Characteristic) 


(Characteristic) — 
+ 


— a 
+ at 
+ oP 
a 
D=0 20 


Fig. 3. The sign of Gt(F+); case (b) with m even 


Thus on yi= 2x the quantity G* takes on positive and negative values in 
the neighborhood of the point as shown in Figure 3 and there must be at least 
two supersonic curves on which Gt=0 and two supersonic curves on which 
G~=0 entering the point. 


3. Ordinary sonic limit point on a sonic limit line (wy, + 0) 


(Examples of such points are found in the radial flow and in the flow discussed 
in [2].) 
Theorem 5. Neither an YZ nor an LZ, can come into the point. 


Proof. The proof of this theorem is presented in [5] and depends on the fact 
that, due to the presence of the sonic limit line (y,=0 on g=q), we have from 
Taylor’s formula applied to y, that y,=O(6q). Then ¢y,//M?—1=0(6q)}, so 
that ye = (qy,//M2—1) —py tends to the value —yy at g=q, irrespective of 
the path of approach. But if the sonic point of an Y, (on which y;=0) is ap- 
proached along the 4, wy; clearly tends to zero. Hence ys=0 if such a point 
belongs to a sonic limit line. A similar argument holds for the 4. 


The geometry of the flow in the neighborhood of such a point was derived in 
[2]; see also [5]. 
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Theorem 6. If w5g+-0, the streamline in the hodograph inflects at the image 
point on passing from the subsonic to the supersonic region. 


Proof. First we note that there is a streamline through the image point. 
This follows from the implicit function theorem since ws-— 0. 

Differentiating Chaplygin’s equation with respect to gq and applying the 
conditions of an ordinary sonic limit point we have that p,,,==0 if (and only if) 
Wo9:0. Hence using g as parameter on the streamline we find that 

Lee ee a ao 


wees pe ccae = a de Ee 
dq ee dq Yo 25 


Theorem 7. Any curve in the physical plane without inflection can be a sonic 
limit line, the two flows which meet at it lying on the convex side. 


Proof. Let the curve be denoted by y =/(«) and define the stream direction # 
at each point by 


(4) cot # = — f(x). 
Considering # as a parameter on the curve: x=g(#), y=/(8), we now set 
(5) F(#) = —cosec® g'(#) = sec Hh’ (9), 

and construct the hodograph solution w(g, 0) which satisfies the conditions 


Os Vo = 2% F (8) on ¢= 4%. 


Then the corresponding flow in the physical plane clearly contains a sonic limit 
line, whose coordinates (see (2)) are given by 


x= — f F(d)sind dd = g (0), 
y= f F(9) cosddd=hi(d), 


when the integration constants are suitably chosen. It is easily checked that 
the flow is unique, except for reversal of direction, since (4) must hold at a sonic 
limit line (see [5]). 

Consider now the map of a radial element at the sonic line in the hodograph 
plane. By differentiating (2), we find dx/dq—dy/dq=0 and 


ee 2 (p= 1)cose 6 (pv —1) Che 1) h! 
dq? oe eee pee 
d? eagh hh ee eee bes 
ae at = sind F(9) = — “7 g'(8), 


where (5) has been used (cf. Eq. (8) on p. 232 of [2]). Hence whether the element 
points into the subsonic or supersonic region, it maps into an element in the 
direction (h'(8), —g’(8)) perpendicular to the sonic limit line. But according to 
(4): (x) g’(8) =cosec? 8> 0, so that g’ has the same sign as f’’. Hence the image 
element points towards the convex side of the curve. 
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4. Higher-order sonic limit point on a sonic limit line (y, = 0) 

Any point at which an Y, or %, meets a sonic limit line is such a point (see 
Theorem 5), and any higher-order sonic limit point is necessarily isolated from 
others of the same kind (since otherwise p(q, #) = constant). In particular, it 
follows that neither an Y, nor an #, can coincide with the sonic line, or in 
other words a sonic limit line cannot be an Y, or #2. 

(i) Higher-order sonic limit point such that py9==0. (An example of a flow 
with this type of singularity was constructed in [2]. Other examples will appear 
in. Part I) Section. 57) 


Theorem 8. An odd number of 4, and an odd number of #, enter the pornt. 


Proof. As in the proof of Theorem 5 we have, since y,=0 along the sonic 
circle, that y= (q p,{VM2— 1) —ws tends to the value —y, at g=g,. Now the 
9= 9, function y; (g,@) is continuous on any sufficiently small 
supersonic ‘‘semicircle’’ C about the point and tends 


to values with opposite signs at the ends (see Figure 4). 


yo 20\ fe Hence there exist an odd number of points on C at 
EM va which y;=0. Similarly for py, =0. 
Yoo. #0 oan, Theorem 9. There is yust. one LZ, and one Lom 
me euce entering the point and their images 1,, 1, together form 
ve SO Gate ao a cusp of the first kind* in the radial direction. 
= 


Proof. Making use of the conditions y,=y»s9= 
Voq=Y79= 9, Yoo=O (note that the sonic limit line 


M=1 se : ‘ ; 
Fig. 4. Higher-order sonic limit point CONdition is embodied in y,,=0), and the fact that 
(v9.9 ©) Yeo O only if y,,,4+0 (cf. Theorem 6), we have 
I L; fy Fo q Lis Ly9 Loa Loge Fy90=0, 
vond 


Foo = — 2y5o( qG )+0, Figg = ©Viqq# O- 
Thus the leading terms in the Taylor expansion of the curve F(q, 0) =0 are 
1 1 
F = 57 B yoo (04) (00) + gy gag (09) +] + =0. 


Hence any branch enters the point either in the g-direction or in the 9-direction. 
For any branch in the g-direction we may write 6g=x, 60—=ux, where u 
tends to zero with x, and in the x, w-plane this branch is determined by 


1 1 1 
G(s) Pte a: [3 Foo!) + a, + Regg too] +::=0. 


However, since 0G/0x = F,,,,/4! == 0, there is a unique curve 
F 
a = f(u) = —12 1" yet... 
99949 


through the origin. Its image (cusp of first kind) 


(59)8 = — 12-22% (gg) 4... 


9999 


* See [4], where the method of our proof is given. 
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in the hodograph plane therefore provides all branches (viz. two) of F=0 which 
enter in the radial direction. 

Similarly we see that there is only one branch entering in each of the positive 
and negative #-directions. These together form the sonic limit line. Thus the 
cusp must be formed by a unique J, and a unique /, (cf. Theorem 8). 


Theorem 10. The sonic limit line in the physical plane is cusped at the point. 


Proof. From the transformation equations (2) we have along a sonic limit 
line (using ? as a parameter) 


dx 1 ; dy 1 
HRY — 0g (— sin ? ws) 6 ae = 0g (cos # wa) 
and 


d®x 1 ay 1 : 
ae og (— cosP yy — sin B Wo»), clare (— sin dps + cosd woo). 


Hence we see that at the point both first derivatives, but not both second deri- 
vatives vanish. 


Theorem 11. The streamline w(q,0)=c which goes through the image point 
in the hodograph plane has a cusp of the first kind there. The cusp tangent has the 
radial direction. 


RO meet p—p—c¢. Lhen 


P= Yq = Vo = Pog = Yq = 9, = Yow =O and hence Wqqq 7 O 


(cf. Theorem 6). The Taylor expansion of ® is 


B = 5 [Boo(9)"] + 57 Pyaa(O9)@ +--+ --- 


Thus according to [4] the curve #(qg, #) forms a single cusp of the first kind with 
cusp tangent in the radial direction. 


(ii) Higher-order sonic limit point such that wgg—0. There exist cases where 
no supersonic limit line enters the point. A necessary condition for this (cf. 
Theorem 8) is that ws does not change sign on g=gq,. An example is presented 
iin Pewee INL, Secor 7 

The sonic limit line need not be cusped. For from the proof of Theorem 10 
we see that dx/d0 =dy/d0 = a?x/d0*= d*y/d#*=0, but if yo99=0, then not both 
third derivatives vanish. 


Theorem 12. If a supersonic limit line enters the point, then its tmage 1s erther 
perpendicular or tangent to the sonic circle. 


Proof. Since the point lies on a sonic limit line, the first nonvanishing deri- 
vative must be of the type wps,..9=-0 (see Section 2). Then the expansion of 


——2 


Theorem 2 holds: 1+n 


F = — "°F (6q) (69)°" + 0 (6g, 9)"**, 


12% 
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and the curve on which F(g, 0) =0 is either tangent (6g—=0) or perpendicular 
(6% =0) to the sonic circle. 

Theorem 13. If ygo-:0, then the streamline p(q,0)=c through the image 
point has a triple point with three branches entering from the subsonic side, and 
one from the supersonic side, all in the radial direction. Two of the branches on 
the subsonic side form a cusp of the first kind. 


Proof. As in Theorem 11, let p=y—c. Then 


P= Yq = Yo = Yaq = Ya0 = Poo Poaq = Vago = Yaoe = 0, 
= a <1 ee od ae 
YoooF 9,  Paqqg=% Pago = — qe e08 0, Pyqqqq= 9, 


and the Taylor expansion about the point is 
1 is 
B(q,9) = 5 [Poo (09)")] + 47 LQ) Pyoa0 (59)? (88)] +--+ + 
+ <7 (0) Perego (9)! (88) +] ++. 


Thus any branch of the curve w(¢g, #) =0 coming into the point has the direc- 
tion 69=0. Now let b=3! Pjqq0/Poo0= —3! (vy +1)/G@<0 and b'=3Y,o970/Pooo 
and put 60=ux, 6g=x. Then 


3!¥(q, 9) 1 
Yoo9 x 


G(%, u) = 7 = w+ x[bu+---]+27[b'u+.---]+---=0, 


and the curve G(x, u) =O enters along bux =O, t.e. x=0 or u=O. 


By repeating the argument, we find as in the proof of Theorem 9 that x =0 
corresponds to a unique curve, whose image in the hodograph 


(09) = — (dq)? + 


forms a cusp of the first kind. Since )<0, this falls on the subsonic side of the 
sonic circle. 


By a similar analysis, “0 may be shown to correspond to a unique curve, 
whose image in the hodograph passes through the point in the radial direction. 


5. Sonic point of a supersonic limit line: Ordinary point (7, = 0) 
(Examples of such a point are found in the flows of RINGLEB and of TEMPLE 
& YARWOOD. See ¢.g. [5].) 


Theorem 14. The supersonic limit line is unique and is tangent to the sonic 
circle at the image point. 


Proof. As in Theorem 1, we have F = F,=0, F,=-0. Hence there is a unique 
line g=/(#) through the point on which F(g, ®)=0, and it has the direction 


d FE 
na = =0, 94g¢=const. 
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Theorem 15. The geometry in both physical and hodograph planes is the 
same as for an ordinary point of a sonic limit line. 


Proof. The directions of the isovel, equipotential, characteristics, streamline 
and isocline in the hodograph plane depend only on the property Y,=0, Yo 0 
at the point. Moreover the effect on these curves of the mapping from the hodo- 
graph into the physical plane depends only on y,=0, y»==0 at the point. (Cf. 
Section 3 and [2].) 


6. Sonic point of a supersonic limit line: Higher-order point (y, = 0) 


(i) First non-vanishing derivative is of the type wp... 9==0. (The simplest example 
Ita 

of this type of flow is one where n=1, 1.e. at the point: y=Wp=Yq=Yq0= 9, 

Ws90. Such a flow can be constructed by superposition of Chaplygin solutions 


as follows: 


Nice Bea Bf aq L 7 (ind + cos) +CA(q) sin 29 


where F(q) = @F(3, —3,3,9°), A= — th G) +A@WileB (a), B=41lea, 
and C = 1/29; & (q)-) 


Theorem 16. The images of the branches of the supersonic limit line come in 
tangentially or perpendicularly to the sonic circle. 


Proof. As in Theorem 2 the Taylor expansion of F about the point is 


F = — "7° #9 (6q) (69)"" +0 (5g, 08)"*?; 


hence the curve on which F=0 has the direction 6g=0 and/or 6?0=0. 
(ii) First non-vanishing derivative (q decreasing) is of the type Wao... 9-0. 


itn 
(The simplest example of this type of flow is one where m = 1, 1.¢. y= Y= Vqq= 9; 
Y,9==0. Such a point is g=g,, ®=O0 in the flow p=k# + (1/q) sin #, where 
k=—4/q.) 


Theorem 17. The images of the branches of the supersonic lint line come in 
perpendicularly to the sonic circle. 


Proof. Let <> yo...9= #0. Then as in Theorem 2 we find the Taylor 


series +n 


F = k'?(60)?" + O(6q, 60)?**1. 
Hence F =0 has the direction 60 =0. 
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Sonic Limit Singularities. Part I: Examples 


G. S. S. LUDFoRD & S. H. ScHOT 


Communicated by HILDA GEIRINGER 


1. Introduction 


It was shown by HILDA GEIRINGER [/] that in addition to the well-known 
supersonic limit lines in the two-dimensional inviscid potential flow of a com- 
pressible fluid there may occur limit lines (exhibiting a different geometry) which 
coincide entirely or in part with the sonic line. This corrected an earlier statement 
by RINGLEB [2], that the only flow exhibiting such a sonic limit line is the well- 
known radial or source flow. 


In the first part of the present paper we gave the general theory of sonic 
limit singularities. Here we discuss special types of sonic limit line flow, providing 
tangible examples of this theory. In particular, three such flows are explicitly 
calculated, the resulting flow patterns being of some interest in themselves. 

In a previous paper [3] the present authors showed how flows containing 
sonic limit lines may be constructed by forming suitable product solutions 
[see (1) in next section] of the linear hodograph equations, for any value of the 
separation constant »20. As an example of a sonic limit line flow with no 
source-like character, we explicitly constructed the flow corresponding to m =1. 
(The case  =0 leads to the source flow.) 

In this paper we investigate the case of general m (which may be supposed = 0) 
and show that the corresponding class of flows has certain well-defined geometrical 
features. For m less than a certain value m, the flow in the physical plane re- 
presents a compressible analogue of the incompressible flow about a corner* 
(with a certain corner angle) and must be envisaged as taking place on a qua- 
druply-sheeted surface. The sheets are joined at a supersonic limit line and at 
the sonic limit line which is a portion of a hypocycloid (n> 1) or epicycloid (<1). 
In certain cases, portions of the same sheet overlap in the large. 

In addition we discuss the shape of the sonic limit line when solutions cor- 
responding to different are superposed and show how then the supersonic 
limit line can be eliminated so that an isolated sonic limit line is obtained. 


Throughout the discussion we will make use of the general results established 
in Part I (this journal, Vol. 2, pp. 160—172), which will be referred to by means 
of the prefix I. 


x For a discussion of the corner flows corresponding to Chaplygin’s product 
solutions, see [4] and references cited there. 
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2. Sonic limit line flows for any n 
It was shown in [3] that flows with a sonic limit line may be constructed for 
any positive real number m, by taking the stream-function to be* 


(1) p(t, 0) = P,(r) sinnd = i" H, (1) sinnd, t= 4, 
where H,,(r) is, in general, a combination of hypergeometric functions, 
(2) H, (1) = % E(t) +7" Fe (r), 

satisfying the equation 

3) t(t—) A+ [(@—1)-@+1—P)r] H, +an(n + 1)8 HR, =0, 


the constant x, in the combination being chosen such that 


dP, a 
(Cae clone ee 


1” H,(t)]=0 atthe sonic value t=1,= a (3 = — : 
Any linear combination of such solutions (1) will also be a solution containing 
a sonic limit line. The discussion of these composite flows will be deferred until 
Section 6. 

Corresponding to this solution in the 1, #-plane (hodograph), the coordinate 
functions, see (I.2), are given parametrically in terms of t and #@ by 


Laos =i een - 
poet (4= 7) aera n P,) cos(n +4) 8 


: (27 sae) + nP,) cos (” — 4) 6 


nm—1 


(5) 
AE eae! —8 / Ge? = 
Y= Zt (4 — 1) | . 3 (27 re nP,)sin(m +1) 8 + 
1 dk, | 7 
es (27 =a n P,) sin (1 1) 9| 


if m=-1. For n=1, the coordinate functions were given in [3] and will not be 
repeated here. The case n =0, leading to the well-known radial or source flow 
will also be omitted in the discussion. 

We now investigate the shape of the sonic limit line in the physical plane 
corresponding to any n>0. From (4) and (5) it is given by** 


= n Fy (74) z i 1 “i 
Xn art(1—z)i lw +1 cos(m + 1) 8 — — cos (n 1) 9), 

o nF, (74) 1 : : 
Vn 2? (1 — 1,)8 aia sin (n +1) 8 + 5 sin(w — 1) 3). 


These are the parametric equations of curves (roulettes) which can be generated 
by rolling a circle on the outside or inside of a fixed circle. For »>1 the ratio 


ie In general p(t, 0) = P,(t) sin(nd + e,) = 13” H,,(t) sin(w@+ «,). Since the in- 
clusion of the phase angle e, merely amounts to a rotation of axes in the physical 


plane, there is no loss in generality in taking «, to be zero throughout the following 
discussion. 


*x* The same result follows from Theorem Me, 
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of the radius of the rolling circle to that of the fixed circle is (7 —1)/2n and the 
curves are hypocycloids. For n<1 this ratio is (1—)/2n and the curves are 
epicycloids. For n=1 (the case discussed in [3]) the curve is a cycloid and for 
n =0 (the radial flow) the roulette degenerates into a circle. When 2n/(n—1) 
is rational, the sonic limit line will close up after a certain number of revolutions. 
Thus for »=4,%, 2,3 the curves are particularly simple, namely a cardioid, 
nephroid, astroid, and tricuspid, respectively. 

In Part I we called a point where y,=ys=0 for t=1; a higher-order sonic 
limit point. Thus for any 7, the sonic limit line will have higher-order points 
gue Yo = 13" H,,(t,) ncosnd =0, 


that is, at the points corresponding to 


(7) palm st on M0, hte se ok 


nN 


It is precisely at these points that the roulettes have their cusps. A sonic limit 
line does not necessarily have cusps at its higher-order points, however. This 
depends on the vanishing of wy, (cf. Section I.4).) 


3. The associated supersonic limit line and streamlines in the hodograph 
In addition to sonic limit lines the flows given by (1) also possess other 
(supersonic) limit singularities (cf. Theorem 1.8). The Jacobian of the trans- 
formation from the hodograph to the physical plane shows that these occur 
where 
(8) Cy Dis = ($4) a kop cr. ty tage ($4) =0. 
GAO. Oil ie 27(1—7)Pt1 \ ad 


d(t, 3) (1 — 7)? 


For the solutions (1) these supersonic limit lines are the curves 


me sel (ae 2t 1 AaB, 
(9) cotnd == | Gp) r—1 8 Be dr 


To investigate the shape of these curves we note from (3) that the function 
P(t) satisfies the equation 


U 


d i ae n* 1— (287) 2 p 
: = QO. 
(10) dt|(4—71)' at 4 ri — 1h 
Hence if we let 
elena 
(11) Si n ‘ P. : dt ’ 


then S,, will satisfy the Riccati equation 


aSy B | 2% 2 RAD cet) tah ee 
(12) dt ie ie EERO Sn v= Bs 
At the sonic circle t=1;= 5 ae we have the conditions 
Sy 
(13) So 0), a 0 


which follow from (4) and (12). 
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Now the expression on the right of (9) becomes indeterminate at t=7,, but 
application of L’Hospital’s rule leads to 


1 
(14) cotn? =0, pa 2M tt on m=0,; 441, £2,..2. 


Hence the supersonic limit line intersects the sonic limit line at the higher-order 
points (7). This was seen to be a general property of sonic limit lines in Section I.4. 
Moreover it was shown in Theorem I.9 that, when wy ,==0, the supersonic limit 
line enters the higher-order point in a cusp. At the maximum circle tr =1, again 
(14) holds. For (2m/n) -42<0<[(2m + 2)/n] - 3a the supersonic limit line must 
extend from t=t, to t=1, in two symmetric branches (corresponding to the 
plus and minus signs in (9)) about the line of symmetry # = [(2m +1)/n] - 2. 

In order to determine the number of sheets required to represent the flow in the 
physical plane it is important to know whether the branches of the supersonic 
limit line in the hodograph plane intersect again on the line #=[(2m+1)/n]-32 
in the interval t,<1t<1. Such supersonic double limit points certainly 
occur for some combinations of F(t) and F®)(r) (see for example CRAGGs’ 
doublet [5]). We shall show, however, that for any m less than a certain value 7, 
which we do not compute explicitly, the flow represented by (2) does not possess 
such points. 

From (9) and (41) such a double limit point would have to occur at a zero 
of S,. Now at the initial point 1, of the interval t,<t1tSX1 we have, on dif- 
ferentiating (12) and using (13), that 

GES n(2B + 1) 
at PAG) Sa) 


== 10) 


so that from (13), S,,(t) is negative in the neighborhood of t,. Thus if the curve 
S,=S,,(t) is continuous in the interval t,<1t<1, it cannot cross or touch the 
t-axis in this interval; for it would have to do so with non-negative slope, whereas 
from (12) we see that at a zero of S, in the supersonic region dS,,/dt would have 
to be negative. 


We now investigate the possibility that S,—=S,,(t) has an infinity in the 
interval t,<t<1. Consider the zero isoclines of (12). The ordinates of these 
in the t, S,-plane are given by 


eee Br+Vp tT? — n?(1— 71) [((2B+1)t — 1] 


‘ n(1—T) 


and are real when 
A = f2 — *(= — )( “ly = 
pr—n are 2B +1 = =) 


Now the maximum of (1/7—1)(26 +1—4/t) occurs at t=1/8 +1 and is Ba 
Hence for any <1, we have 4=0 in the whole interval, and the zero isoclines 
extend from t=0 to t=1 with non-positive slope. From the sketch (Fig. 1) 
it is then clear that the integral curve S,,(t) which satisfies the conditions (13) 
at t=Tt, cannot have an infinity in the range t,<1t<1, for it must lie below 
the t-axis and cannot cross the zero isoclines, which therefore bound it from 
below. Thus S,<0 for all .<1<1. 
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Consider now the case n>1. For »>1, we have 1<0 for some values of 1, 
and the zero isoclines lie as shown in Fig. 2. Thus the integral curve S,, (zt) which 
satisfies (13) is no longer bounded from below by the isoclines, and may have a 
singularity at some point 7) in the interval t,<1<1 such that lim S,,(t) =— oo 


TT. 


and iim S,(z ) =-++ oo. As shown in the figure, the branch on which lim SA) 


tT Tt 


= -+ co will then cross the t-axis with negative slope, as required by (12). 

It can be shown from the asymptotic behavior of (12) that such infinities 
do in fact occur for large enough values of m; moreover, when there is an infinity 
(and hence a zero) for one value of 7 there is also for all larger values. The smallest 


tt 


Fig. 1*: 1< 1 Fig. 2*. 1<m<n 


value of m for which the integral curve has such a singularity will be denoted by 
MN, (just proved greater than unity). This value 1) for which S,,(t) first has zeroes 
will not be calculated explicitly. However, the example worked out in Section 5 
for the case m =2 exhibits no supersonic double limit points, so that we may 
conclude > 2. 

From these results it is clear then, that for any »<m, the supersonic limit 
line in the hodograph plane in a flow given by (41) consists, for (2m/n) -31<0< 
[(2m +2)/n]-42, of a drop-shaped portion with the line #=[(2m-+1)/n]-32 
as axis of symmetry (see for example Fig. 7). We will not discuss the shape of 
the supersonic limit line when n= 7. 

We now turn to a discussion of the streamlines in the hodograph plane for 
the case »<m,. The straight streamlines are given by dy=y,dt+y,d0 =0, 
a0 =0 which implies y,=0, or from (1): 


It 
G=m-—, m=O, Ap De here 
nN 


By symmetry we need only consider the circular sector 0<0<a/n. If the 
straight streamlines are interpreted as boundaries in the physical plane it is 
evident that (1) represents flow about a corner of angle 7—(z/n). The total 
angie through which the stream turns is 7/n. 


* The plus and minus signs indicate regions in which the integral curves have 
positive and negative slope, respectively. 
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Let us now consider the non-straight streamlines in the hodograph. Dif- 
ferentiating (1) along a streamline 


dy = Gr inn Odt +uP,cosndd? =0, 
dt 

so that the inclination « of the streamline to the radius vector at any point of 

the hodograph is given by 


bee ee 
OO aa ea tae 2T oH 
and on substituting (1) we find 
cot nm 
(15) COU Ore ree oe 


Now we have seen above that 
SO @at, Cae: 5S, <0-throughout” a, 215 
By a similar argument to that in the proof above, it can be shown that 


S,<0 throughout 0<t< 1; 
so that forn<mp, 


sign(cot «) =sign(cotn#) throughout 0<1T<1. 


Since all streamlines pass through the origin, it is clear then, that as # increase 
from 0 to z/n the streamlines form closed loops, which are symmetric about the 
line } = 2/2 (see e.g. Fig. 7). The curvature of a streamline changes continuously, 
and becomes zero at the sonic circle t = 1; (in accord with the theory of the sonic 
limit line, see [3] and Theorem 1.6). The only exception to this general behavior 
of the non-straight streamlines is the streamline, say y,, which passes through 
the higher-order sonic limit point at t=1,, ®=a2/2n. This is given by 


Ye = P,(r) sinn (=) = P, (x) 


and is known to have a cusp at that point (Theorem I.11). It divides the stream- 
lines which do not intersect a limit line — sonic, /;, or supersonic, / — from those 
that do, and will be called the critical streamline y,. Another streamline of 
interest is the one passing through the point t=1, ®=a/2n. Given by 


Wa — i (1) sin n (=) = Ine (1) ? 


it will be called the asymptotic streamline, for reasons which will become clear 
in the next section. 


The streamlines in the hodograph plane may then be divided into three classes: 
(a) those which intersect neither the J, nor the J: y,<p< , 


(b) those which intersect both the J, and the 7: Pa<py; 
(c) those which intersect only the J,: 0<p<y,- 


The behaviors of these three types of streamlines in the physical plane are radi- 
cally different. 
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4. The physical plane 


Consider now the flow in the physical plane corresponding to the hodograph 
solution in the typical circular sector 0<1t<1, 0<#<a/n. From the geometry 
in the neighborhood of limit lines (sonic and supersonic) we know that the stream- 
lines are cusped there and consequently the regions to one side are multiply 
covered. Thus the limit lines in the physical plane delineate regions in which 
the flow changes its multi-valued character. 
This can be envisaged as a folding of the 
physical plane into a multiply-sheeted con- 
figuration with folds at the limit lines, so 
that each sheet is simply covered by the 
appropriate portion of the total flow, and 
at the limit line folds the flow passes from 
one sheet to the appropriate next sheet. 
Due to our choice of the representative 
circular sector in the hodograph the entire 
multi-sheeted configuration representing the Va 
physical plane is of course bounded by the 
map of the sides of the sector, 7.e. the 
straight streamlines ? =0, #=a/n. 


tT=0 T= 6=0 


Fig. 3. Circular sector in the hodograph plane 


Accordingly, for 7<,, the physical plane must be envisaged as a quadruply- 
sheeted surface, the sheets corresponding to the regions in the hodograph bounded 
by (see Fig. 3): 


(I) O<t<t, ®B=0; t=, Oa; Qe ieee 


2 — h Te a S= i ae 
i eeteerr ete OS Oa Lb; t=%, Oar 


Milesy = 7,: 


Ope wy = eee eed OS PT, 
2n n nN 2n n 


Examples (Section 5, Figs. 6, 9) show that portions of sheet I may overlap. 


We have already discussed the shape of the sonic limit line in the physical 
plane. The supersonic limit line maps into a wedge-shaped curve with a cusp 
at the higher-order sonic limit point, the cusp direction coinciding with that of 
the sonic limit line. 

The streamlines in the physical plane behave as follows*. Initially all 
streamlines have the #=0 direction. They then bend around in a clockwise 
direction and finally all assume the #=<z/n direction. Between these two final 
states their behavior depends on the group to which they belong in our classifi- 
cation of streamlines in the hodograph plane (see end of Section 3): 


(a) These streamlines bend smoothly about the cusp of the YF, and never 
leave sheet I. 


* Examples are given by Figs. 4, 5 andy]. o: 
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(b) These start on sheet I, are cusped at the Y, and pass onto sheet IT, con- 
tinue on sheet II until they meet the &, reflect onto sheet III, continue on 
sheet III until they reach the &, pass there to sheet IV, continue on sheet IV 
up to the &, and reflect back onto sheet I. 

(c) These start on sheet I, reflect at one branch of the % onto sheet II and 
continue on this sheet to infinity (the circle t=1 maps into infinity in the 
physical plane); they then return on sheet IV, and reflect at the other branch 
of the & back onto sheet I. 


Note that the dividing streamline between group (b) and group (c), namely y,,. 
asymptotically approaches the two branches of the supersonic limit line. 

An example (1=1) having the features described above was worked out 
in [3]. Other such flows will be constructed explicitly in the next section. 


5. The sonic limit line flows for n= 2 and n= 4 


As a representative example of the above results consider n =2, y=. The 


straight streamlines are #=0 and ? =4$z2, corresponding to a compressible flow 
about a 90° corner. This particular choice of is exceptional, in that one of the 
solutions of the hypergeometric equation (3) is a polynomial, and the other is. 
an infinite series involving a logarithmic term. 


The streamfunction in this case is taken to be 


(146) y=A(r) sin 20 = |x, ct POE, —3,3; 7) + rt POs, —5, —1; t)| sin 20, 


where 


Bae Se tad 6 7] 


2 
(= Ce eal ee rs 120 wy (vy +4)F(v—S5) 
re ee F(2,— 3: l Os 2 ? 

i Tee ed oe Oe TQ) 2 To age 
eel d = d 
Mg |e ROVER) Se 123.42). 


The coordinate functions are, from (5), 


H_ = tT 4(1 — 1) 8 (3 (t B — B) cos 38 + (cB + B) cosd] 


(17) 1 i 
=T (1 — 7) 4[§ (eB — P) sin 38 — (xB + B) sin 9] 


MY 


tw 


and the sonic limit line in the physical plane is consequently the astroid. 
%_,= — 3[t, *(1 — %)§ B(z,)] (cos 39 — 3 cos 3), 
Yeo=—3([t, *(1—7,) 73 P,(z,) | (sin 3.9 + 3 sin 8), 


or, on eliminating #, 


(%_ + y»)§ 4 (2 Yo)3 =([4 V2 77 3(4 = 1%)? B, (a) ]8. 
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a Se a 
—wimit line 


y=4.2256 


some limit line y=0 
Fig. 4. Case nm = 2. Hodograph plane 


y-0 p=4.2256 p=ls y=22.586 p=28 
U 


y=33 


y=0 


Y=4.2256 


limit ine____. 


y=/5, 


Y=22.586 


y=28 


W=33 Y=38 
Fig. 5. Case m = 2. Physical plane 


The supersonic limit line in the hodograph is given by 


1 Tete te 
= Ze 


(18) taiw2@) == 


t 


and its image in the physical plane is obtained from (17). 


The critical streamline is 
the (Tg 215 86), 
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and the asymptotic streamline is* 
5 SP ak Seay at 
Ya = Plt) =m FOZ, 35) 1) + FO 5 ht) 


P'(3) L'( a) 


= ee Ee 
mel Ie al a Rey (ETE == ADS O- 


) 1"(6) (4) 
The streamline corresponding to any y=Yp is obtained from (16) and (17). 

Graphs of the limit lines and 
selected streamlines in the hodograph 
and in the physical plane, respectively, 
are presented in Figs. 4 and 5**. The 
manner in which the sheets corre- 
sponding to the physical plane are 
connected is indicated diagrammati- 
cally in Fig. 6. 

As another example of the fore- 
going, this time exhibiting a sonic limit 
line of the epicycloidal type (namely 
a nephroid), we select the case n =. 
Here both solutions of the hyper- 
Fig. 6. Case »=2. Sketch of multi-sheeted surface corre- eooimcute equation ee unexceptional 

sponding to physical plane and are represented by infinite series: 


Minit line 


(19) p(t, 8) =P (x) sing 0=[x, tt F® (a, b, 3; t) +7 -t FO (a—h, b—$, 3; 1)| sin $0 


where 
a,b=4(—4+)31), 
Bota b, 2. ae T(2) SS P(v+a)I(v +b), 
2 LI’ (a) I'(b) =r I'(vy + 3) yp! : 
Pog 2 poet ee OO Gee aie a eee 
(a , ait] @-pre—H 2 I'(v+)y! (Gis 
d = 9 d 1 
aes lag FO) | (al Be _, = 342054. 


The mapping to the physical plane is given by 


pout 741 — 1)-4[— 3 (27 BY — $B) cos$0 + (27 B’ + $B) cosh 9], 
1 
B 


yy = 4(1 — 1) 4 [—4 (27 Bi — $B) sin $0 + (27 Ph’ +3) sin3 9), 


* The expansion 


Dn 2 6) ic a==5 
F(a, b,c; ean pero F(a,b,1+a+b—c¢c;1—1)+ 


I'(a+b—c)TI(c) =p 
++ — T'(a) I'(b) ~(1 — 1)* @ F(¢—a,c—b,1+¢—a—b;1—1) 
is used near t = 1. 
** In the actual computation of these flows use was made of the tables in [6] 
and [7], taking due account of the errors in [7] which were pointed out in [8]. 
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and the sonic limit line is the nephroid 


X= 3, *(1— 7%) 2 A(z) (cos $B + cos39) 


Ve=$t, (1 — 7%) *B(z,) ($sin $0 + sin$ 9). 


sonic limit line 


Pe 
/ 


an 


~ a 
linitme = ¥ 9° 


Fig. 7. Case m = 4. Hodograph plane 


y-44s 7. 

ed ; 

Minit line Innit line 
y=410 y=309 PIAS w=34S W=369 W=#I0 


Fig. 8. Case »=%4. Physical plane 


The supersonic limit line in the physical plane is obtained from (20) and 


(en eel) Ore os 
(21) tan 0 = / 


ete 1—T Pe 


The critical streamline is 


Ve = P. (t) = 3.685 76, 
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and the asymptotic streamline is 


po 
ES tra 2 3 37076. 


Pa = P(t) = my F(a, 0, 35 1) Fa — 3,0 2, 33 4) 
aS 


=y Cm al 
VPA 
The streamlines in the physical plane are obtained from (19) and (20). 

Graphs of the limit lines and selected streamlines in the hodograph and 
physical plane for this case are shown in Figs. 7 and 8. Fig. 9 is a diagrammatic 
representation of the sheeted surface corresponding to the physical plane. 

An interesting feature of both these flows is that sheet I overlaps itself (cf. [2] 
and [3]). 


- sore lint line 


init line 
a 


~ 


Fig. 9. Case n = 4. Sketch of multi-sheeted surface corresponding to physical plane 


6. Superposition of sonic limit line flows for different n. 
The isolated sonic limit line 


Superposition of a finite number of sonic limit line solutions (4) for different n 
will again yield a flow with a sonic limit line. The coordinate functions may 
be found by adding the coordinate functions (5) of the individual flows. Natu- 
rally, the simple geometric character of the sonic limit line in the physical plane 
for the separate solutions is thereby destroyed. 


However, a simple rule for the geometrical construction of this line may still 
be given. Draw the sonic limit lines which correspond to the individual solutions. 
Then add the coordinates of the points on the individual sonic limit lines at 
which the slopes of the latter are the same. The curve constructed in this manner 
will be the sonic limit line of the combined solution. This follows from the fact 
that a sonic limit line is always perpendicular to the streamline, so that one is 
in fact adding for the same value of #, 7.e. for the same point in the hodograph. 

In Section 3 we noted that, for any n, the flow constructed from (1) necessarily 
contained a supersonic limit line in addition to the sonic limit line. The question 
arises whether, by superposition, a flow can be constructed which is devoid of 
any supersonic limit singularities and contains only a sonic limit line. 
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Now given the solution (1) for any ”, we can, by adding to it a suitable radial 
flow solution (7 =0), obtain a flow such that (i) there is still a sonic limit line 
and (ii) ~» has the same sign along the sonic line. Hence the sonic line consists 
entirely of ordinary limit points (y,=+-0) and it must therefore be an isolated 
sonic limit line by Theorem I.5. In fact we shall see that even in the limiting 
case of the above construction, when ws is allowed to take zero values (without 
changing sign: py y=0, cf. Theorem I.8), the flow still has an isolated sonic limit 
line (see Section I.4 (ii); an example will be presented in Section 7). 


For any , we have for the solution (1) 
yo (t:, 8) = P,(t)) ncosnd, 
whereas for the radial flow, whose streamfunction is 


(22) WIT). = RO, 


we have 


We (ze ?) = k. 


If we add (1) and (22) and require that ws does not change sign anywhere on 
T =T,, we must choose & such that 


(23) |k|=+nP,(t,) accordingas P,(t,) SO. 


Clearly the addition of the radial flow does not affect the vanishing of y, on 
the sonic line; hence the resulting flow will also possess a sonic limit line. Moreover 
if the strict inequality holds in (23) we have yw =-0 for t =t, and the sonic limit 
line is isolated. 

Consider now what happens to the supersonic limit line, which appears in 
the flow (1), when the radial solution (22), with | k| increasing from zero, is added. 
From the streamfunction of the combined flow 


w(t, 0) =kd + P(t) sinnd 


we see that the supersonic limit line is given by 


(24) k+Pincosn? + \PQOP, sinnd? =0, 
see (8), where 

ee ar(1— ret oT epee 
(25) (2B+1)7™—1’ Q (1 —7)? 


Equation (24) may be written 


Bs k 
cos [nF 6(t)] = ings eran 


hence for given k, only those values of t for which 
(26) A=n?P?+ POP?—k>0 and y4St=1 


will give real values of #. 
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With » fixed, we investigate the region in the k, t-plane for which (26) is 
satisfied. Certainly the as is bounded by the lines t=7, and t=1. Dif- 
ferentiating k2 =n? P? + PQ P,” with respect to t and making use of Chaplygin’s 
equation (I.1) in the form PO Pi’ + PQ’ P, +n? P,=0, we have 


a) _ 4 (2p? + PQ Pp’) =Pi*(P'Q — PQ’) =— BEGPEUT pir <o, 


jhe) ~~) Ghee [(2B + 1) —1}? 


so that the boundary curve k?=n?P?+PQP,* of the ee (26) is mono- 
tonically decreasing with increasing t in the range tStS1. A sketch of 
the region (26) for k>0 is presented in 
Fig. 10. 

From this graph we may deduce the 
extent of the supersonic limit line in 
the hodograph. Clearly for k=O the 
supersonic limit line extends from t = 7, 
to t=1. This remains true for values 
oh kup to k=|nP; sh However as k 
increases from \nP, (1 )| to |nP,(z,)| the 
supersonic limit line detaches cou the 
maximum circle t=1 (while remaining 
7 AAG eee attached to the sonic circle t=t,) and 
Fig. 10. Extent of the supersonic limit line in the hodo- 2 : 4 

graph plane for each k : shrinks towards the sonic circle. For values 

of & such that |2P,(1)|<k<|nP,(z,)| 

the supersonic limit line thus extends only up to a circle of radius t,,, say, and 
Tm—>t, aS k-+|nP,(t,)|. For |nP,(z,)|<k no real branches of the supersonic 
limit line exist. 


kt 


pao 


24) (== 


7. A flow with isolated sonic limit line 
As a simple example of a flow with an isolated sonic limit line (see Section 6) 
we superpose the sonic limit line flow for »=1 (see [3]) and the radial flow 


neice p=kO+ cle + (1 — 7+] sind, 


where 
geen c/a CaN ae eer, 
x 2st ( 1.267876 for rere 
since e (t= ‘ar Gas a k must be chosen such that 
Sep 4 (B+1 Ba 4) 2p. \P 


(see (23)). In view of our earlier study of higher-order sonic limit singularities 
(Section I.4 (ii)) and the preceding section, we choose k to be exactly 


a Apreaoh (= 1.811629 for y=), 


then 


Yo =k+ et ¥[x + (4 — 1)'*1] cos@, 


Yoo = —T *[x + (1 — 2)?*1] sin@, 


and the point t =1,, 8 =0 is a higher-order point on a sonic limit line at which 
Yoo= 0. 
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The coordinates in the physical plane are given by (cf. [3]) 
2% =[xrt(4—t) P+ RB4 t*] cos 20 + 
+ xB lo — 1)? + dz + (28 +1) logt + 2kt~4(1 — 1) FP cos®, 
2y =|[xt4(1— 1) 8 +8 +74] sin 20 — (6 +1) (29) + 2kc-4(1 — 1)* sin ®, 


Ha 
yp=2 Y=O0.5 
T=05 

y=07 

‘ y=0.07 

ee Retr ica 
sonie limit line ae) 
oF y=—0.0/ 
Ne vo 
ered Die 

W=—O7. 


=-05 
ee Mb : 


Fig. 11. Isolated sonic limit line. Hodograph plane 


where the integral in the first equation can be evaluated in closed form for 
y =. The sonic limit line in the physical plane is then given by 


x = (8 +1) cos#[2 — cos #] + constant, 
y = (8 +1) sind [2 — cos d] — (6B +1) 0. 
This curve is the locus of a point moving harmoni- y=O0rEe 


cally on the extended diameter of a circle, while 
the circle in turn rolls on a straight line. It is a 


trochoid-like curve and does not possess any cusps Cp om 


(cf. Section I.4 (ii)). 

The flow in the hodograph plane, and hence oe y=0-€ 
also in the physical plane, is of quite a different Shue eee 
character from that of the sonic limit line flows streamline y=0 
discussed previously. 

In the hodograph plane the streamlines no longer form lemniscate-like loops 
inside a circular sector but behave as shown in Fig. 11. For a given value of 
y=Yo> 0, the streamline in the hodograph consists of two portions: (7) one 
above the horizontal axis which meets the sonic limit line and ranges from t =0 
through all subsonic and all supersonic values of t to the maximum value t =1, 


188 G. S. S. Luprorp & S. H. ScHor: 


and (ii) one below the horizontal axis which starts at T=0 and loops around 
back to tr =0 without attaining supersonic speeds. The two portions which make 
up the streamline y = —wy)<0 form the symmetric image of y=Yo in the axis 
=0. This explains the symmetry of the figure. The streamline y =0 falls into 


aN 
Ha < : | | 
\ | / 
\ \ | ; | 
‘ ‘ | | I 
: \ | | I 
‘ \ | 
‘i | | 
\ | | 
Wy | | / 
ee SS I j 
Pe Me | / 4 ] 
i. \ | / if / 
~ / / 
= MD, 3 OS: fi 
Ve Ne Q y 
NF 
7 
Sake. ; 
a a eS / 
Raa SSS ae 
Sai 
| 
| 
i 
Sa a eos =ured 
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Fig. 13a and b, Isolated sonic limit line. Physical plane (schematic) 


the above pattern in a limiting sense (see Fig. 12); it degenerates into the straight 
segment #=0 plus a symmetric cusped curve which passes through the higher- 
order limit point t=t,, ®=0 (cf. Theorem I.13). 

The corresponding streamlines in the physical plane therefore consist of 
subsonic-supersonic portions which are cusped at the sonic limit line, and smooth, 
purely subsonic portions. 

To simplify the discussion of the streamlines in the physical plane we restrict 
ourselves first to that portion of the plane which corresponds to the upper half 
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(007) of the hodograph plane. In this range lie the cusped portions of some 
streamlines (y20) and the smooth portions of the rest. Thus we expect the 
streamlines in the physical plane corresponding to this half of the hodograph 
to lie on two sheets (a subsonic and a supersonic sheet), some entirely on the 
subsonic sheet and the others partially on both sheets. From this and the varia- 
tions of # along the streamlines in the hodograph we infer Fig. 13a. Note that 


Ta 


Fig. 14. Isolated sonic limit line. Sketch of multi-sheeted surface corresponding to physical plane 


parts of the subsonic sheet overlap in a manner similar to that in the examples 
of Section 5. By symmetry the portion of the sheeted surface corresponding to 
the lower half of the hodograph is as shown in Fig. 13b. 


To obtain the complete picture we must join the two parts depicted in Figs. 13a 
and 13b along the line of symmetry PQR; the result is Fig. 14. Note that the 
entire surface consists of only two sheets — a subsonic and a supersonic sheet. 
However, the individual sheets overlap in parts, so that there are regions of 
single-, triple-, and quintuple-valuedness in parts of the x, y-plane. 
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On Exterior Boundary Value Problems 
in Linear Elasticity 


R. J. DUFFIN & WALTER NOLL 


Communicated by E. STERNBERG 


1. Introduction 


The questions considered in this note were suggested by a paper by FINN & Nott 
on the uniqueness and non-existence of Stokes flow [/]. The equations of linear 
isotropic elasticity are only slightly more complicated than the equations govern- 
ing Stokes flow, and it proves possible to extend the treatment of Finn & NoLi 
to analogous elastic problems. 

The equations for the displacement field q of a linearly elastic body at rest 
and under no body forces may be written (cf. [2]) in the form 


(1.1) Aq=egradp, divq=—cpb. 


Here c is the positive constant c =1— 20, where a is Porsson’s ratio. 

If we set c=O in (1.1) and interpret q as the velocity field, then (4.1) are 
the equations of Stokes flow. The results of [7] may be phrased this way: Let 
the region of flow be bounded internally by closed surfaces (curves) 8. We seek 
solutions of (1.1) (with c=0) which, in three dimensions, tend uniformly to zero 
at infinity or, in two dimensions, are bounded at infinity. There is at most 
one solution q of this type which assumes a given constant value q=qp at the 
internal boundary 8. We note that this theorem remains valid also if the pres- 
cribed boundary value gq =q)j is not constant, although a simple physical inter- 
pretation for q as the flow past a rigid body can be given only when q is constant. 

In this paper the corresponding theorem is proved for the case when c>0 
in (1.1). In this way, we get three different uniqueness theorems for boundary 
value problems in elasticity: 

(a) Consider an infinite elastic medium bounded internally by closed surfaces 
%. There is at most one displacement q of the medium which assumes prescribed 
values q =q, at the internal boundaries 8 and which tends uniformly to zero 
at infinity. 

(b) Consider an infinite elastic medium bounded internally by an infinite 
cylinder % (not necessarily circular). There is at most one plane displacement q 
of the medium in the plane perpendicular to the cylinder which assumes pre- 
scribed values g=q, at 8 and which is bounded. 
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(c) Consider a thin infinite elastic plate bounded internally by closed con- 
tours 8. There is at most one displacement q in the plane of the plate which 
assumes prescribed values q = at 8 and which is bounded. All three results 
are special cases of Theorem 1, proved in Section 2. The proposition (b) cor- 
responds to the case of plane strain. (c) corresponds to a state of generalized 
plane stress. In this case, g must be interpreted as the average displacement, 
the average being taken over the thickness of the plate. As is well known [3, 
p. 208], q then satisfies the two-dimensional form of the equations (1.1) when 


Porsson’s ratio o is replaced by o’ = = Thus, in this case, the constant c must 
be given the positive value c=1—20'= Ss. We also prove a uniqueness 
theorem for the transverse deflection w of an infinite thin plate bounded 
internally by closed curves 8. The function w is biharmonic: A Aw =0. We show 
that, if the plate is clamped in a prescribed way at 8, so that w and dw/én are pre- 
scribed at 8, then there is at most one deflection w which is of the order O (7) asy>ov. 

The essential difference between the results obtained here and the uniqueness 
theorems previously known lies in the conditions at infinity. Previously, com- 
plicated regularity conditions at infinity had to be imposed on q or w. We 
merely require simple limit conditions. It should not be difficult to obtain 
similar results when the surface tractions, rather than the displacements, are 
prescribed at the interior boundaries. 

As in the case of Stokes flow so also in elasticity we observe the striking 
difference between the situations in two and three dimensions. In three dimensions 
we get a unique solution if we require that q(oo) =0, 7.e. that the “boundary 
at infinity” is fixed. If we impose the same condition in two dimensions there 
is, in general, no solution at all. For a physical interpretation consider, for 
example, a very large thin plate with a fixed external boundary. At first glance 
our theorem seems to imply that it is not possible to impose a rigid displacement 
GY =Q=const. on the internal boundaries, a result that contradicts physical 
experience. However, the correct interpretation is the following: The gradient 
of the displacement q, and hence the strain, is very small everywhere in the 
plate. In the limit as the external boundary expands to infinity the gradient 
of the displacement tends to zero. The displacement q converges to its internal 
boundary value q y. But this convergence is not uniform, and the condition 
that q=0 on the external boundary is not preserved in the limit. A similar 
explanation must be given in the case of the transverse deflection of a large 
plate. If the external boundary is clamped (w =@w/én =0) we get a unique 
solution. In the limit, however, the conditions w =é@w/dn=0 at the “boundary 
at infinity” get lost. 

If we impose a displacement q = qp) on the boundary of an internal cavity 
in a large elastic body with fixed external boundary, we shall get a displacement q 
which approximates the solution for the infinite medium determined uniquely 
by the condition q(co) =0. This solution will be given explicitly in Section 4 
for q)= const. and a spherical cavity. The displacement gradient is not small, 
and the strains and stresses will remain finite in the limit as the body becomes 
infinite. In this case, the boundary condition q=0 on the external boundary 
is preserved in the limit. 
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2. The basic lemmas 


Let q be a vector field defined and three times continuously differentiable 
in a neighborhood of infinity. 


Lemma 1. In three dimensions, assume that q=q,+0(1). Then 


(2.1) Acunl G)s==0 © amplies5 curl ¢ = O(r-”),, 
and 
(222) Aide g) == 0 implies “div =O (7~") ; 


Lemma 2. In two dimensions, assume that q is bounded. Then 


(2:3) “(curl Giz 0,12 implesy ~-curl g = O(7 4) , 
and 
(2.4) A(divig) =Oe  implies., div q—=O0 (r+): 


The statements concerning the curl are Lemmas 4 and 5 in [7]. The propositions 
concerning the divergence can be proved in precisely the same manner as these 
lemmas. 

Lemma 3. Assume that w is a scalar field defined and biharmonic in a neigh- 
borhood of infinity. Then w=O(r) implies that grad w is bounded. 

Proof. We denote the mean value of a function w over the circumference 


of a circle of radius 7 with center at Q by m,(w,7, Q). If w is biharmonic in a 
region containing this circle we have the following mean value theorem (cf. [4], 


p.-316,°(14).*)3 

(2.5) 3w(Q) = 4m9(w, 37, Q) — mo(w,7, Q). 

We may replace w by grad w in (2.5) because grad w is biharmonic if w is. If 
we do so, multiply then by 7, integrate between y=0 and and finally apply 
GREEN’S theorem, we get 

(2.6) 37 grad w(Q) = 8m, (w, 27, 2) — my (w, 7, Q), 

where 


mM, (w,7, Q) =< f#@ +re)edo. 


We have é 
|, (w,7, Q)| S max | w (P) 


’ 


where P varies over the disc of radius y around Q. Hence, by (2.6) 
r| grad w(Q)| <3 max |w(P)|. 


Applying this inequality to circles of increasing radius with centers that tend 
to infinity, we obtain the assertion of the lemma. 


x It is proved there in the case of three dimensions, but both the proof and the 
result are applicable in any dimension. The lemma stated here is also valid in higher 


dimensions. 
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3. Uniqueness of the exterior boundary value problems 


Theorem 1. Let & be a region in space (in the plane) exterior to an internal 
boundary % consisting of a finite number of piecewise smooth closed surfaces 
(curves). Then there is at most one vector field q, defined and twice continuously 
differentiable in &, with the following properties: 


(a) q is a solution of the equations 
(3.1) Aq=gradp, divaq=—cp, 
where ¢ is a given positive constant. 

(b) q is continuous up to 8 and assumes prescribed values 
(3.2) = Gana x. 


(c) In three dimensions q converges uniformly to zero at infinity: q (co) =0 
(in two dimensions, q is bounded). 


Proof. We observe that (3.1) implies 
(3.3) curl curlg = — (1 +c) grad p. 
Operating on (3.3) with curl and div yields 
(3.4) Afoml¢) =O zha=0: 


Here we have made use of the fact that q has continuous third derivatives. 
It is shown in [2] that this is a consequence of (a). 


Since the equations (3.1) are linear and homogeneous, it is sufficient to show 
that q)=0 implies g=0. We do this in the three-dimensional case. 


By Lemma 1 and (3.4) we have 
(3.5) culg =O(r2), > =——divg=0(r). 
We use the GREEN’S identity 
(3.6) J (curl q)?dV re. (q- curl curl q) dV = curlg:(qxn)dS, 


valid for any smooth vector field defined in the finite region t with piecewise 
smooth boundary G&. 


By (3.3) we have 
q-curlcurl gq = — (1 +c) q-gradp 
=—(1 +0) [div(pq) — p diva]. 

By use of (3.1), and the divergence theorem, the identity (3.6), applied to the 
region &, between the internal boundary 8 and a sphere G, of radius 7 containing 
%, gives 

(3.7) [| (carla)? + + (aiv 9)" dV=—f [(1+e)pq-n+teurlg-(qxn)]d5. 

%, S, 

Here we have used the fact that q=q,=0 on %. It follows from (3.5) and 
q=0(1) that the right-hand side of (3.3) tends to zero as ro. 
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Hence es 
i (curl q)? +++" (divg)|av =0. 
B 
Since the integrand is non-negative and continuous, it must vanish. Hence 
curl g =0, div g =0 whence Aq =0. Since g =00n B and g =0(1) the maximum 
and minimum principle for harmonic functions implies q = 0, g.e.d. 

The proof for the two-dimensional case is the same as that of Theorem II in 
[1] when modified along the lines indicated above. 


Theorem 2. Let & be a region in the plane exterior to an internal boundary 
% consisting of a finite number of piecewise smooth curves. Then there is at 
most one scalar field w with the following properties: 


(a) w is defined and biharmonic in 8: AAw=0. 
(b) w and its first partial derivatives are continuous up to 8 and w and 
dw/On assume prescribed values 
ow a 
Oye = | al B. 
(c) w=O(7) as_7— 00. 
Proof. We may regard w as the stream function of a plane viscous flow q. 
Then 
q = (grad w)", 


where | indicates the operation of rotating a vector counterclockwise by a 
right angle. The condition (b) implies that grad w, and hence q assumes pre- 
scribed values qg =q, at 8. By Lemma 3 it follows from (c) that q is bounded. 
The results of [/] imply that q, and hence grad w, is uniquely determined. This 
implies the uniqueness of w because the boundary values w =wy, at ¥ are given. 


4. A representation formula 
The representation theorem (Theorem III) of [/] may be extended to the 
case of the three-dimensional equation (1.1) in the following way: Let q be a 
solution of (1.1) defined in a neighborhood of infinity € which is star-shaped 
relative to a point Q. Then there is a harmonic vector-field wu, defined in © and 
harmonic at infinity, such that 


(4.1) q=u + grady, 
where 
mes ea a + [ stdj 
(4.2) p=vQ0+re = pag" | s divu(Q + se)ds. 


0 
The proof is identical to the one given in [/]*. Using the representation (4.1), 
we may obtain the displacement q caused by a rigid translation gq = q) = const. 
of a spherical cavity in an infinite elastic medium. We try 


(4.3) u(Q +re) =a qo +B 5 [Go— 3 (Me) €], 


* The equations (4.3), (4.4) and (6.2) in [1] contain a mistake in sign. 
Arch, Rational Mech. Anal., Vol. 2 13b 
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oA. 


which is harmonic at infinity. Let Q be the center of the sphere and a its radius. 


With (4.1) and (4.2) we can determine the constants « and f so that the 
boundary condition q(Q +ae) =q is satisfied. The result of the calculation is 
(ri=re) 


(44) g= za tan fao[3(t +20) 2 +(2)] +3r(q0-”) [4 — SI 


For c =0 this reduces to the classical solution of Stokes for the flow past a sphere. 
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